skip to main content

Title: Exploring the Structural, Dynamic, and Functional Properties of Metal‐Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic‐level modeling in research on metal‐organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster‐based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high‐throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Advanced Functional Materials
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A host of important performance properties for metal–organic frameworks (MOFs) and other complex materials can be calculated by modeling statistical ensembles. The principle challenge is to develop accurate and computationally efficient interaction models for these simulations. Two major approaches are (i) ab initio molecular dynamics in which the interaction model is provided by an exchange–correlation theory ( e.g. , DFT + dispersion functional) and (ii) molecular mechanics in which the interaction model is a parameterized classical force field. The first approach requires further development to improve computational speed. The second approach requires further development to automate accurate forcefield parameterization. Because of the extreme chemical diversity across thousands of MOF structures, this problem is still mostly unsolved today. For example, here we show structures in the 2014 CoRE MOF database contain more than 8 thousand different atom types based on first and second neighbors. Our results showed that atom types based on both first and second neighbors adequately capture the chemical environment, but atom types based on only first neighbors do not. For 3056 MOFs, we used density functional theory (DFT) followed by DDEC6 atomic population analysis to extract a host of important forcefield precursors: partial atomic charges; atom-in-material (AIM) C 6 , C 8 , and C 10 dispersion coefficients; AIM dipole and quadrupole moments; various AIM polarizabilities; quantum Drude oscillator parameters; AIM electron cloud parameters; etc. Electrostatic parameters were validated through comparisons to the DFT-computed electrostatic potential. These forcefield precursors should find widespread applications to developing MOF force fields. 
    more » « less
  2. Abstract

    The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

    more » « less
  3. null (Ed.)

    Precipitation strengthening of alloys by the formation of secondary particles (precipitates) in the matrix is one of the techniques used for increasing the mechanical strength of metals. Understanding the precipitation kinetics such as nucleation, growth, and coarsening of these precipitates is critical for evaluating their hardening effects and improving the yield strength of the alloy during heat treatment. To optimize the heat treatment strategy and accelerate alloy design, predicting precipitate hardening effects via numerical methods is a promising complement to trial-and-error-based experiments and the physics-based phase-field method stands out with the significant potential to accurately predict the precipitate morphology and kinetics. In this study, we present a phase-field model that captures the nucleation, growth, and coarsening kinetics of precipitates during isothermal heat treatment conditions. Thermodynamic data, diffusion coefficients, and misfit strain data from experimental or lower length-scale calculations are used as input parameters for the phase-field model. Classical nucleation theory is implemented to capture the nucleation kinetics. As a case study, we apply the model to investigate γ″ precipitation kinetics in Inconel 625. The simulated mean particle length, aspect ratio, and volume fraction evolution are in agreement with experimental data for simulations at 600 °C and 650 °C during isothermal heat treatment. Utilizing the meso-scale results from the phase-field simulations as input parameters to a macro-scale coherency strengthening model, the evolution of the yield strength during heat treatment was predicted. In a broader context, we believe the current study can provide practical guidance for applying the phase-field approach as a link in the multiscale modeling of material properties.

    more » « less
  4. Abstract

    Additive manufacturing, no longer reserved exclusively for prototyping components, can create parts with complex geometries and locally tailored properties. For example, multiple homogenous material sources can be used in different regions of a print or be mixed during printing to define properties locally. Additionally, heterogeneous composites provide an opportunity for another level of tuning properties through processing. For example, within particulate-filled polymer matrix composites before curing, the presence of an applied electric and/or magnetic fields can reorient filler particles and form hierarchical structures depending on the fields applied. Control of particle organization is important because effective material properties are highly dependent on the distribution of filler material within composites once cured. While previous work in homogenization and effective medium theories have determined properties based upon ideal analytic distributions of particle orientations and spatial location, this work expands upon these methods generating discrete distributions from quasi-Monte Carlo simulations of the electromagnetic processing event. Results of simulations provide predicted microarchitectures from which effective properties are determined via computational homogenization.

    These particle dynamics simulations account for dielectric and magnetic forces and torques in addition to hydrodynamic forces and hard particle separation. As such, the distributions generated are processing field dependent. The effective properties for a composite represented by this distribution are determined via computational homogenization using finite element analysis (FEA). This provides a path from constituents, through processing parameters to effective material properties. In this work, we use these simulations in conjunction with a multi-objective optimization scheme to resolve the relationships between processing conditions and effective properties, to inform field-assisted additive manufacturing processes.

    The constituent set providing the largest range of properties can be found using optimization techniques applied to the aforementioned simulation framework. This key information provides a recipe for tailoring properties for additive manufacturing design and production. For example, our simulation results show that stiffness for a 10% filler volume fraction can increase by 34% when aligned by an electric field as compared to a randomly distributed composite. The stiffness of this aligned sample is also 29% higher in the direction of the alignment than perpendicular to it, which only differs by 5% from the random case [1]. Understanding this behavior and accurately predicting composite properties is key to producing field processed composites and prints. Material property predictions compare favorably to effective medium theory and experimentation with trends in elastic and magnetic effective properties demonstrating the same anisotropic behavior as a result of applied field processing. This work will address the high computational expense of physics simulation based objective functions by using efficient algorithms and data structures. We will present an optimization framework using nested gradient searches for micro barium hexaferrite particles in a PDMS matrix, optimizing on composite magnetization to determine the volume fraction of filler that will provide the largest range of properties by varying the applied electric and magnetic fields.

    more » « less
  5. Abstract

    Flexible metal–organic Frameworks (MOFs) are an interesting class of materials due to their diverse properties. One representative of this class is the layered‐pillar MOF DUT‐8(Ni). This MOF consists of Ni2paddle wheels interconnected by naphthalene dicarboxylate linkers and dabco pillars (Ni2(ndc)2(dabco), ndc = 2,6‐naphthalene–dicarboxylate, dabco = 1,4‐diazabicyclo‐[2.2.2]‐octane). DUT‐8(Ni) undergoes a volume change of over 140% upon adsorption of guest molecules. Herein, a ligand field molecular mechanics (LFMM) study of the CO2‐induced flexibility of DUT‐8(Ni) is presented. LFMM is able to reproduce experimental and DFT structural features as well as properties that require large simulation cells. It is shown that the transformation energy from a closed to open state of the MOF is overcompensated fivefold by the host–guest interactions. Structural characteristics of the MOF explain the shape of the energy profile at different loading states and provide useful insights to the interpretation of previous experimental results.

    more » « less