Ab Initio Derived Force Fields for Zeolitic Imidazolate Frameworks: MOF-FF for ZIFs

An automated protocol to construct flexibility parameters for classical forcefields: applications to metal-organic frameworks

In this work, forcefield flexibility parameters were constructed and validated for more than 100 metal-organic frameworks (MOFs). We used atom typing to identify bond types, angle types, and dihedral types associated with bond stretches, angle bends, dihedral torsions, and other flexibility interactions. Our work used Manz's angle-bending and dihedral-torsion model potentials. For a crystal structure containing N atoms in its unit cell, the number of independent flexibility interactions is 3(N atoms - 1). Because the number of bonds, angles, and dihedrals is normally much larger than 3(N atoms - 1), these internal coordinates are redundant. To reduce (but not eliminate) this redundancy, our protocol prunes dihedral types in a way that preserves symmetry equivalency. Next, each dihedral type is classified as non-rotatable, hindered, rotatable, or linear. We introduce a smart selection method that identifies which particular torsion modes are important for each rotatable dihedral type. Then, we computed the force constants for all flexibility interactions together via LASSO regression (i.e., regularized linear least-squares fitting) of the training dataset. LASSO automatically identifies and removes unimportant forcefield interactions. For each MOF, the reference dataset was quantum-mechanically-computed in VASP via DFT with dispersion and included: (i) finite-displacement calculations along every independent atom translation mode, (ii) geometries randomly sampled via ab initio molecular dynamics (AIMD), (iii) the optimized ground-state geometry using experimental lattice parameters, and (iv) rigid torsion scans for each rotatable dihedral type. After training, the flexibility model was validated across geometries that were not part of the training dataset. For each MOF, we computed the goodness of fit (R-squared value) and the root-mean-squared error (RMSE) separately for the training and validation datasets. We compared flexibility models with and without bond-bond cross terms. Even without cross terms, the model yielded R-squared values of 0.910 (avg across all MOFs) ± 0.018 (st. dev.) for atom-in-material forces in the validation datasets. Our SAVESTEPS protocol should find widespread applications to parameterize flexible forcefields for material datasets. We performed molecular dynamics simulations using these flexibility parameters to compute heat capacities and thermal expansion coefficients for two MOFs.

Coarse-grained force field for ZIF-8: A study on adsorption, diffusion, and structural properties

Metal-organic frameworks (MOFs) are revolutionizing a spectrum of industries, from groundbreaking gas storage solutions to transformative biological system applications. The intricate architecture of these materials necessitates the use of advanced computational techniques for a comprehensive understanding of their molecular structure and prediction of their physical properties. Coarse-grained (CG) simulations shine a spotlight on the often-neglected influences of defects, pressure effects, and spatial disorders on the performance of MOFs. These simulations are not just beneficial but indispensable for high-demand applications, such as mixed matrix membranes and intricate biological system interfaces. In this work, we propose an optimized CG force field tailored for ZIF-8. Our work provides a deep dive into sorption isotherms and diffusion coefficients of small molecules. We demonstrate the structural dynamics of ZIF-8, particularly how it responds to pressurization, which affects its crystal structure and leads to local changes in aperture size and area. Emphasizing the game-changing potential of CG simulations, we explore the characteristics of amorphization in ZIF-8. Through computational exploration, we aim to bridge the knowledge gap, enhancing the potential applications of nanoporous materials for various applications.

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.

Deciphering the controlling factors for phase transitions in zeolitic imidazolate frameworks

Zeolitic imidazolate frameworks (ZIFs) feature complex phase transitions, including polymorphism, melting, vitrification, and polyamorphism. Experimentally probing their structural evolution during transitions involving amorphous phases is a significant challenge, especially at the medium-range length scale. To overcome this challenge, here we first train a deep learning-based force field to identify the structural characteristics of both crystalline and non-crystalline ZIF phases. This allows us to reproduce the structural evolution trend during the melting of crystals and formation of ZIF glasses at various length scales with an accuracy comparable to that of ab initio molecular dynamics, yet at a much lower computational cost. Based on this approach, we propose a new structural descriptor, namely, the ring orientation index, to capture the propensity for crystallization of ZIF-4 (Zn(Im)2, Im = C3H3N2-) glasses, as well as for the formation of ZIF-zni (Zn(Im)2) out of the high-density amorphous phase. This crystal formation process is a result of the reorientation of imidazole rings by sacrificing the order of the structure around the zinc-centered tetrahedra. The outcomes of this work are useful for studying phase transitions in other metal-organic frameworks (MOFs) and may thus guide the development of MOF glasses.

Development of a Transferable Force Field between Metal-Organic Framework and Its Polymorph

Conventionally, force fields for specific metal-organic frameworks (MOFs) are derived from quantum chemical simulations, but this method can be computationally intensive, especially in cases for large MOF structures. In this work, we devise a methodology to reduce the force field derivation costs by replacing the original MOF with a smaller polymorphic structure, with the hypothesis that the force field parameters will be transferrable among chemically identical, polymorphic MOF structures. Specifically, we demonstrate this transferability in MOF-177 structure for H2O and NH3 gas molecules and show that the force field parameters derived from a smaller polymorphic MOF-177 can be used accurately to the original MOF-177 structure. This methodology can accelerate the development of force field parameters for large porous materials, in which computational costs for conventional methods are expensive.

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.

Coarse-grained versus fully atomistic machine learning for zeolitic imidazolate frameworks

Zeolitic imidazolate frameworks are widely thought of as being analogous to inorganic AB2 phases. We test the validity of this assumption by comparing simplified and fully atomistic machine-learning models for local environments in ZIFs. Our work addresses the central question to what extent chemical information can be "coarse-grained" in hybrid framework materials.

The influence of force fields on the structure and dynamics of water confined in ZIF-8 from atomistic simulations

The complexity of modeling flexible crystals, such as ZIF-8, mainly stems from the handling of intramolecular interactions. Numerous force fields have been proposed in the literature to describe the interactions between atoms in ZIF-8. We employ seven force fields to examine the structure and dynamic behavior of water molecules confined in ZIF-8, with the aim of investigating the impact of force fields on simulation results. Various structural characterization methods consistently indicate that the choice of different force fields has quantitative effects but no qualitative effects on the structural characteristics of confined water. Additionally, the force fields do not impact the qualitative description of the diffusion mechanism. Both mean-square displacement and van Hove autocorrelation function reveal two characteristic movements of water molecules diffusing in ZIF-8: a short-time intra-cavity hopping process and a long-time inter-cavity hopping process. However, the framework flexibility is found to play a crucial role in determining the order of spatial arrangement and local structure, self-diffusion coefficient and reorientational dynamics of confined water. Specifically, the DREIDING force field gives rise to an unrealistic stiff framework, enhancing the order of spatial arrangement and diminishing the local ordered structure of confined water. Meanwhile, it results in much slower translational and reorientational dynamics. Hence, the general DREIDING force field cannot be considered for providing a quantitative description of the water structure and dynamics.

Sub-Micrometer Phonon Mean Free Paths in Metal-Organic Frameworks Revealed by Machine Learning Molecular Dynamics Simulations

Metal-organic frameworks (MOFs) are a family of materials that have high porosity and structural tunability and hold great potential in various applications, many of which require a proper understanding of the thermal transport properties. Molecular dynamics (MD) simulations play an important role in characterizing the thermal transport properties of various materials. However, due to the complexity of the structures, it is difficult to construct accurate empirical interatomic potentials for reliable MD simulations of MOFs. To this end, we develop a set of accurate yet highly efficient machine-learned potentials for three typical MOFs, including MOF-5, HKUST-1, and ZIF-8, using the neuroevolution potential approach as implemented in the GPUMD package, and perform extensive MD simulations to study thermal transport in the three MOFs. Although the lattice thermal conductivity values of the three MOFs are all predicted to be smaller than 1 W/(m K) at room temperature, the phonon mean free paths (MFPs) are found to reach the sub-micrometer scale in the low-frequency region. As a consequence, the apparent thermal conductivity only converges to the diffusive limit for micrometer single crystals, which means that the thermal conductivity is heavily reduced in nanocrystalline MOFs. The sub-micrometer phonon MFPs are also found to be correlated with a moderate temperature dependence of thermal conductivity between those in typical crystalline and amorphous materials. Both the large phonon MFPs and the moderate temperature dependence of thermal conductivity fundamentally change our understanding of thermal transport in MOFs.

Ab initio methods for the computation of physical properties and performance parameters of electrochemical energy storage devices

With the rapid development of electric vehicles and mobile technologies, there is a high demand for electrochemical energy storage devices and electrochemical energy conversion devices. Devices meeting these needs include metal-ion batteries (MIBs), supercapacitors (SCs), electrochromic devices (ECDs), and multifunctional devices such as electrochromic batteries and supercapatteries. Currently, the goal has been the enhancement of operational parameters and physical properties that results in a higher performance of these devices. In the case of batteries, SCs, and supercapatteries, scientists seek to improve the equilibrium voltage, energy density, power, capacitance, and charge rate. In the case of ECDs, the focus is on improvement of the optical modulation and coloration efficiency. However, synthesis and characterization of new materials, or of materials with optimized properties, is time consuming and highly expensive. Computational simulation of materials can expedite the experimental endeavor by modelling novel atomic structures and predicting device performance. This is possible using ab initio theories and applying physical principles that allow us to understand the underlying mechanisms governing the behavior of materials in these devices. Taking as a point of departure density functional theory (DFT), in this review, we discuss the first principles methods used for the computation of physical properties and performance parameters of electrochemical energy storage devices. A wide coverage of DFT is given, dealing with the strengths and weaknesses of the most popular functionals used in the field of electrochemical energy storage. With these tools, ab initio methods for the computation of basic properties such as effective mass, mobility, optical band gap, transmissivity, conductivity (ionic and electronic), and criteria for structure stability (cohesive energy, formation energy, adsorption energy, and phonon frequency) are addressed. We also highlight the first principles techniques for the calculation of performance parameters in MIBs, SCs, and ECDs.

Computer simulation of the early stages of self-assembly and thermal decomposition of ZIF-8

We employ all-atom well-tempered metadynamics simulations to study the mechanistic details of both the early stages of nucleation and crystal decomposition for the benchmark metal–organic framework (MOF) ZIF-8. To do so, we developed and validated a force field that reliably models the modes of coordination bonds via a Morse potential functional form and employs cationic and anionic dummy atoms to capture coordination symmetry. We also explored a set of physically relevant collective variables and carefully selected an appropriate subset for our problem at hand. After a rapid increase of the Zn–N connectivity, we observe the evaporation of small clusters in favor of a few large clusters, which leads to the formation of an amorphous highly connected aggregate. [Formula: see text] and [Formula: see text] complexes are observed with lifetimes in the order of a few picoseconds, while larger structures, such as four-, five-, and six-membered rings, have substantially longer lifetimes of a few nanoseconds. The free ligands act as “templating agents” for the formation of sodalite cages. ZIF-8 crystal decomposition results in the formation of a vitreous phase. Our findings contribute to a fundamental understanding of MOF’s synthesis that paves the way to controlling synthesis products. Furthermore, our developed force field and methodology can be applied to model solution processes that require coordination bond reactivity for other ZIFs besides ZIF-8.

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.

Combined Deep Learning and Classical Potential Approach for Modeling Diffusion in UiO-66

Modeling of diffusion of adsorbates through porous materials with atomistic molecular dynamics (MD) can be a challenging task if the flexibility of the adsorbent needs to be included. This is because potentials need to be developed that accurately account for the motion of the adsorbent in response to the presence of adsorbate molecules. In this work, we show that it is possible to use accurate machine learning atomistic potentials for metal-organic frameworks in concert with classical potentials for adsorbates to accurately compute diffusivities though a hybrid potential approach. As a proof-of-concept, we have developed an accurate deep learning potential (DP) for UiO-66, a metal-organic framework, and used this DP to perform hybrid potential simulations, modeling diffusion of neon and xenon through the crystal. The adsorbate-adsorbate interactions were modeled with Lennard-Jones (LJ) potentials, the adsorbent-adsorbent interactions were described by the DP, and the adsorbent-adsorbate interactions used LJ cross-interactions. Thus, our hybrid potential allows for adsorbent-adsorbate interactions with classical potentials but models the response of the adsorbent to the presence of the adsorbate through near-DFT accuracy DPs. This hybrid approach does not require refitting the DP for new adsorbates. We calculated self-diffusion coefficients for Ne in UiO-66 from DFT-MD, our hybrid DP/LJ approach, and from two different classical potentials for UiO-66. Our DP/LJ results are in excellent agreement with DFT-MD. We modeled diffusion of Xe in UiO-66 with DP/LJ and a classical potential. Diffusion of Xe in UiO-66 is about a factor of 30 slower than that of Ne, so it is not computationally feasible to compute Xe diffusion with DFT-MD. Our hybrid DP-classical potential approach can be applied to other MOFs and other adsorbates, making it possible to use an accurate DP generated from DFT simulations of an empty adsorbent in concert with existing classical potentials for adsorbates to model adsorption and diffusion within the porous material, including adsorbate-induced changes to the framework.

Mixing ReaxFF parameters for transition metal oxides using force-matching method

We present the development of the method for the refitting the ReaxFF parameters for a system consisting of the mixed transition metal oxides. We have tested several methods allowing to calculate the differences between the vectors of the forces acting on atoms obtained from the reference DFT simulation and the parameters-dependent ReaxFF. We conclude that the footrule method yields the best parameters among the investigated options. We then validate the parameters using the system consisting of the hematite supported (TiO₂)_n clusters. The results indicate the refitted parameters allow to obtain acceptable geometries of the clusters upon MD simulation on the ReaxFF level, and despite the short timescale lead to the stable structures.

Structural complexity in Prussian blue analogues

We survey the most important kinds of structural complexity in Prussian blue analogues, their implications for materials function, and how they might be controlled through judicious choice of composition. We focus on six particular aspects: octahedral tilts, A-site 'slides', Jahn-Teller distortions, A-site species and occupancy, hexacyanometallate vacancies, and framework hydration. The promising K-ion cathode material K_xMn[Fe(CN)₆]_y serves as a recurrent example that illustrates many of these different types of complexity. Our article concludes with a discussion of how the interplay of various distortion mechanisms might be exploited to optimise the performance of this and other related systems, so as to aid in the design of next-generation PBA materials.

Advanced strategies in Metal-Organic Frameworks for CO2 Capture and Separation

The continuous carbon dioxide (CO₂ ) gas emissions associated with fossil fuel production, valorization, and utilization are serious challenges to the global environment. Therefore, several developments of CO₂ capture, separation, transportation, storage, and valorization have been explored. Consequently, we documented a comprehensive review of the most advanced strategies adopted in metal-organic frameworks (MOFs) for CO₂ capture and separation. The enhancements in CO₂ capture and separation are generally achieved due to the chemistry of MOFs by controlling pore window, pore size, open-metal sites, acidity, chemical doping, post or pre-synthetic modifications. The chemistry of defects engineering, breathing in MOFs, functionalization in MOFs, hydrophobicity, and topology are the salient advanced strategies, recently reported in MOFs for CO₂ capture and separation. Therefore, this review summarizes MOF materials' advancement explaining different strategies and their role in the CO₂ mitigations. The study also provided useful insights into key areas for further investigations.

ResponZIF Structures: Zeolitic Imidazolate Frameworks as Stimuli-Responsive Materials

Zeolitic imidazolate frameworks (ZIFs) have long been recognized as a prominent subset of the metal-organic framework (MOF) family, in part because of their ease of synthesis and good thermal and chemical stability, alongside attractive properties for diverse potential applications. Prototypical ZIFs like ZIF-8 have become embodiments of the significant promise held by porous coordination polymers as next-generation designer materials. At the same time, their intriguing property of experiencing significant structural changes upon the application of external stimuli such as temperature, mechanical pressure, guest adsorption, or electromagnetic fields, among others, has placed this family of MOFs squarely under the umbrella of stimuli-responsive materials. In this review, we provide an overview of the current understanding of the triggered structural and electronic responses observed in ZIFs (linker and bond dynamics, crystalline and amorphous phase changes, luminescence, etc.). We then describe the state-of-the-art experimental and computational methodology capable of shedding light on these complex phenomena, followed by a comprehensive summary of the stimuli-responsive nature of four prototypical ZIFs: ZIF-8, ZIF-7, ZIF-4, and ZIF-zni. We further expose the relevant challenges for the characterization and fundamental understanding of responsive ZIFs, including how to take advantage of their flexible properties for new application avenues.

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.

High-rate nanofluidic energy absorption in porous zeolitic frameworks

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.

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 Zn₂(bdc)₂(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.

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.

Density Functional Theory Studies on Zeolitic Imidazolate Framework-8 and Ionic Liquid-Based Composite Materials

The identification of suitable density functional methods for predicting the properties of nanoporous composite materials is highly significant in the field of chemical and material sciences. The stability of the composite materials depends on the nature of bonding and dispersive interaction at the interface. Thus, we have studied the effect of dispersion correction in the incorporation of hydrophobic and hydrophilic ionic liquids (ILs) into zeolitic imidazolate framework-8 (ZIF-8) nanostructures using the density functional theory (DFT)-based approaches. These structures were analyzed employing selected methods (Becke-Lee-Yang-Parr and Perdew-Burke-Ernzerhof) with dispersion correction (D2 or D3) and different basis sets (such as double-zeta valence polarized (DZVP), triple-zeta valence polarized (TZVP), and triple-zeta valence doubly polarized (TZV2P)) for the understanding of microscopic features of IL@ZIF-8 nanopores. It is found that the result obtained from DFT-D2/TZVP is more reliable for the prediction of the experimental crystal structure as well as stability and spectral information of the complexes. Furthermore, the microscopic analysis of geometries reveals that ILs are highly dispersed and stabilized at the nanopores of ZIF-8, particularly the ZIF-8 structure is highly preferable for the hydrophobic group in ILs. It is found that fluorine-containing anions are highly dispersed on the ZIF-8 surface compared to the nonfluorinated anion (i.e., [BMIM]⁺[Cl]^-). This is confirmed from the adsorption energies (E _ads), charge transfer, electron density analyses, and IR spectral analysis. These findings can provide more insights into the stability of composite materials, which are suitable for applications of catalytic conversion at the confined state, gas storage, and separation techniques.

Tuning the Electric Field Response of MOFs by Rotatable Dipolar Linkers

Recently the possibility of using electric fields as a further stimulus to trigger structural changes in metal-organic frameworks (MOFs) has been investigated. In general, rotatable groups or other types of mechanical motion can be driven by electric fields. In this study we demonstrate how the electric response of MOFs can be tuned by adding rotatable dipolar linkers, generating a material that exhibits paraelectric behavior in two dimensions and dielectric behavior in one dimension. The suitability of four different methods to compute the relative permittivity κ by means of molecular dynamics simulations was validated. The dependency of the permittivity on temperature T and dipole strength μ was determined. It was found that the herein investigated systems exhibit a high degree of tunability and substantially larger dielectric constants as expected for MOFs in general. The temperature dependency of κ obeys the Curie-Weiss law. In addition, the influence of dipolar linkers on the electric field induced breathing behavior was investigated. With increasing dipole moment, lower field strengths are required to trigger the contraction. These investigations set the stage for an application of such systems as dielectric sensors, order-disorder ferroelectrics, or any scenario where movable dipolar fragments respond to external electric fields.