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.

Functionality Optimization for Singlet Fission Rate Screening in the Full-Dimensional Molecular and Inter- molecular Coordinate Space

Functionality Optimization for Singlet Fission Rate Screening in the Full-Dimensional Molecular and Intermolecular Coordinate Space

Optimizing Conical Intersections Without Explicit Use of Non-Adiabatic Couplings

Efficient Bayesian Inference of Atomistic Structure in Complex Functional Materials

A new configurational bias scheme for sampling supramolecular structures

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

Dynamic Cooperativity of Chromophores in Promoting Singlet Fission in Perylenediimide Crystal

Efficient Enumeration-Selection Computational Strategy for Adaptive Chemistry

Molecular level interaction, electrons excitation properties, solvent effect using IEFPCM investigation, topological surface, and docking analysis of 4-pyrrolidin-2-yl-pyridine

Molecule superstructures for computer-aided molecular and process design

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

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs

Quantum mechanical and molecular mechanical simulation approaches bridging length and time scales for simulation of interface reactions in realistic environments

Interplay Between Non-adiabatic Dynamics and Frenkel Exciton Transfer in Molecular Aggregates: Formulation and Application to a Perylene Bismide Model

Diffusion Generative Models for Designing Efficient Singlet Fission Dimers

Computational Investigation on Structural and Reactive Sites (HOMO-LUMO, MEP, NBO, NPA, ELF, LOL, RDG) Identification, Pharmacokinetic (ADME) Properties and Molecular Docking Investigation of (E)-4-((4-chlorobenzylidene) amino) Benzene Sulfonamide Compound

Large-Scale Computational Screening of Molecular Organic Semiconductors Using Crystal Structure Prediction

Maximizing Singlet Fission by Intermolecular Packing.

Accurate interaction energies by spin component scaled Möller-Plesset second order perturbation theory calculations with optimized basis sets (SCS-MP2mod): Development and application to aromatic heterocycles

Computational Screening of Organic Semiconductors: Exploring Side-Group Functionalisation and Assembly to Optimise Charge Transport in Chiral Molecules

Optimisation of ultrafast singlet fission in 1D rings towards unit efficiency