Abstract:Various kinds of flexibility have been observed in metal-organic frameworks, which may originate from the topology of the material or the presence of flexible ligands. The construction of free energy profiles describing the full dynamical behavior along the phase transition path is challenging since it is not trivial to identify collective variables able to identify all metastable states along the reaction path. In this work, a systematic three-step protocol to uniquely identify the dominant order parameters for structural transformations in flexible metal-organic frameworks and subsequently construct accurate free energy profiles is presented. Methodologically, this protocol is rooted in the time-structure based independent component analysis (tICA), a well-established statistical modeling technique embedded in the Markov state model methodology and often employed to study protein folding, that allows for the identification of the slowest order parameters characterizing the structural transformation. To ensure an unbiased and systematic identification of these order parameters, the tICA decomposition is performed based on information from a prior replica exchange (RE) simulation, as this technique enhances the sampling along all degrees of freedom of the system simultaneously. From this simulation, the tICA procedure extracts the order parameters-often structural parameters-that characterize the slowest transformations in the material. Subsequently, these order parameters are adopted in traditional enhanced sampling methods such as umbrella sampling, thermodynamic integration, and variationally enhanced sampling to construct accurate free energy profiles capturing the flexibility in these nanoporous materials. In this work, the applicability of this tICA-RE protocol is demonstrated by determining the slowest order parameters in both MIL-53(Al) and CAU-13, which exhibit a strongly different type of flexibility. The obtained free energy profiles as a function of this extracted order parameter are furthermore compared to the profiles obtained when adopting less-suited collective variables, indicating the importance of systematically selecting the relevant order parameters to construct accurate free energy profiles for flexible metal-organic frameworks, which is in correspondence with experimental findings. The method succeeds in mapping the full free energy surface in terms of appropriate collective variables for MOFs exhibiting linker flexibility. For CAU-13, we show the decreased stability of the closed pore phase by systematically adding adsorbed xylene molecules in the framework.

Predicting kinetics of polymorphic transformations from structure mapping and coordination analysis

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Polymorphic phase transitions: Macroscopic theory and molecular simulation

The rules of metastability: detailed transformation mechanisms in chemical elements by means of molecular dynamics techniques

Functional materials discovery using energy–structure–function maps

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Free-energy landscape of polymer-crystal polymorphism

Crystallization and Vitrification Kinetics by Design: The Role of Chemical Bonding

Collective Variables for Conformational Polymorphism in Molecular Crystals

The computational prediction of pharmaceutical crystal structures and polymorphism

A Hybrid Machine Learning Approach for Structure Stability Prediction in Molecular Co-crystal Screenings

A Robust Crystal Structure Prediction Method to Support Small Molecule Drug Development with Large Scale Validation and Prospective Studies

Structural transformations driven by local disorder at interfaces

Rationalizing polymorphism with local correlation-based methods: A case study of pnictogen molecular crystals

Selective formation of metastable polymorphs in solid-state synthesis

Finite temperature string by K-means clustering sampling with order parameter as collective variables for molecular crystal: application to polymorphic transformation between β-CL-20 and ε-CL-20

Deciphering Cryptic Behavior in Bimetallic Transition Metal Complexes with Machine Learning

Nudged-Elastic Band Calculations of Polymorph Transitions and Solid-State Reactions in Molecular Crystals

Protocol for Identifying Accurate Collective Variables in Enhanced Molecular Dynamics Simulations for the Description of Structural Transformations in Flexible Metal–Organic Frameworks