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.

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

Machine learned Force-Fields for an ab-initio Quality Description of Metal-Organic Frameworks

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

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 Learning Interatomic Potentials for Amorphous Zeolitic Imidazolate Frameworks

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

Towards a Generalizable Machine-Learned Potential for Metal-Organic Frameworks

Transferable Force Field for Metal–Organic Frameworks from First-Principles: BTW-FF

Transferable Force Field for Metal-Organic Frameworks from First-Principles: BTW-FF

A general forcefield for accurate phonon properties of metal-organic frameworks

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

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

A Refined Set of Universal Force Field Parameters for Some Metal Nodes in Metal-Organic Frameworks

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

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

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

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

A foundation model for atomistic materials chemistry

Efficient Machine Learning Force Field for Large-Scale Molecular Simulations of Organic Systems