Abstract: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.

Quantum-accurate machine learning potentials for metal-organic frameworks using temperature driven active learning

Quantum-Accurate Machine Learning Potentials for Metal-Organic Frameworks using Temperature Driven Active Learning

Machine learning potentials for metal-organic frameworks using an incremental learning approach

Machine learned force-fields for an Ab-initio quality description of metal-organic frameworks

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

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

Machine Learned Potential for High-Throughput Phonon Calculations of Metal-Organic Frameworks

Modeling the structural and thermal properties of loaded metal-organic frameworks. An interplay of quantum and anharmonic fluctuations

Quantum Informed Machine-Learning Potentials for Molecular Dynamics Simulations of CO2's Chemisorption and Diffusion in Mg-MOF-74.

Parametrization of Non-Bonded Force Field Terms for Metal-Organic Frameworks Using Machine Learning Approach

Sub-micrometer phonon mean free paths in metal-organic frameworks revealed by machine-learning molecular dynamics simulations

DFT-Quality Adsorption Simulations in Metal-Organic Frameworks Enabled by Machine Learning Potentials

Simulating CO2 diffusivity in rigid and flexible Mg-MOF-74 with machine-learning force fields

Sequential Design of Adsorption Simulations in Metal-Organic Frameworks

Machine Learning Interatomic Potentials for Amorphous Zeolitic Imidazolate Frameworks

Thermomechanical Properties of Two-Dimensional Covalent Organic Frameworks Studied by Machine-Learning Molecular Dynamics

Using Machine Learning and Data Mining to Leverage Community Knowledge for the Engineering of Stable Metal–Organic Frameworks

Using Machine Learning and Data Mining to Leverage Community Knowledge for the Engineering of Stable Metal-Organic Frameworks

Structure and electronic structure of Metal-Organic Frameworks within the Density-Functional based Tight-Binding method

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