Structure and Synthesizability of Iron-Sulfur Metal-Organic Frameworks

Andrew L Ferguson,Jianming Mao,Ningxin Jiang,Andrea Darù,Alexander S Filatov,Jessica E Burch,Jan Hofmann,Simon M Vornholt,Karena W Chapman,John S Anderson

DOI: https://doi.org/10.26434/chemrxiv-2024-swb7z

2024-11-18

Abstract:Sulfur-based metal-organic frameworks (MOFs) and coordination polymers (CPs) are an emerging class of hybrid materials that have received growing attention due to their potential magnetic, conductive, and catalytic properties with applications in electrocatalysis and energy storage. In this work, we report a high-throughput virtual screening protocol to predict the synthesizability of candidate metal-sulfur MOFs/CPs by computing the thermodynamically stable structures resulting from a particular combination of metal cluster, linker, cation, and synthetic conditions. Free energies are computed using all-atom classical mechanical thermodynamic integration. Low-free energy structures are refined using ab initio density functional theory, and pair distribution functions and powder X-ray diffraction patterns calculated to complement and guide experimental structure determination. We validate the computational approach by retrospective predictions of the stable structure produced by experimental syntheses, and a subsequent screen predicts Fe4S4-BDT--TPP as a new thermodynamically stable 1D CP comprising a redox-active Fe4S4 cluster, 1,4-benzenedithiolate (BDT) linker, and tetraphenylphosphonium (TPP) countercation. This material is experimentally synthesized and the 1D chain structure of the crystal confirmed using microcrystal electron diffraction. The computational screening pipeline is generically transferable to neutral and ionic MOFs/CPs comprising arbitrary metal clusters, linkers, cations, and synthetic conditions, and we make it freely available as an open source tool to guide and accelerate the discovery and engineering of novel porous materials.

Chemistry

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a computational workflow for predicting the synthetic feasibility and structural characteristics of iron - sulfur - based metal - organic frameworks (MOFs) and coordination polymers (CPs). Specifically, the authors proposed a high - throughput virtual screening protocol, aiming to predict the synthetic feasibility of candidate iron - sulfur - based MOFs/CPs by calculating the thermodynamically stable structures under specific combinations of metal clusters, linkers, cations, and synthesis conditions. The core assumption of this protocol is that, under given environmental conditions, the structure with the highest thermodynamic stability (i.e., the lowest free energy) is the most likely to be synthesized.

To achieve this goal, the research team adopted a series of methods:

1. **Structure Generation**: A library of structural prototypes containing 4 - coordinated Fe₄S₄ clusters and 2 - coordinated 1,4 - benzenedithiolate (BDT) linkers was constructed using the ToBaCCo - 3.0 code, generating 390 3D and 22 2D structures. Together with the previously reported 1D CP, there were a total of 390 putative MOF structures.

2. **Classical Mechanics Free - Energy Calculation**: Partial charge assignment, neutralizing cation addition, structural relaxation were carried out for each generated structure, and its relative free energy was evaluated. These calculations used molecular mechanics optimization and molecular dynamics simulations to determine the most stable structure under specific conditions.

3. **Density Functional Theory (DFT) Optimization**: The low - free - energy structures screened by classical mechanics were subjected to DFT optimization and energy calculation to further refine the structures and verify their stability.

4. **Experimental Verification and Characterization**: Based on computational predictions, a new Fe₄S₄ - BDT - TTP material was synthesized, and its 1D chain - like crystal structure was confirmed using microcrystal electron diffraction (MicroED).

5. **Surrogate Model Development**: To improve the efficiency of high - throughput virtual screening, the research team also explored low - cost surrogate models based on geometric or energy descriptors to predict the free energy of MOFs, thereby reducing the computational cost.

Through these steps, the research team not only verified the effectiveness of the computational method but also successfully predicted and synthesized a new iron - sulfur - based MOF/CP material, demonstrating the potential of this computational workflow in accelerating the discovery and engineering design of new porous materials.