We, for the first time, systematically investigated the crystal structures, adsorption properties, and microscopic mechanism of CO2 capture with ethylenediamine (en)-appended isostructural M2(dobpdc) materials (M = Mg, Sc-Zn), using spin polarized density functional theory (DFT) calculations. The binding energies of en range from 142 to 210 kJ/mol. The weakest binding materials are en-Cr2(dobpdc) and en-Cu2(dobpdc). Two typical models, the pair model and the chain model, have been considered for CO2 adsorption. Generally, the chain model is more stable than the pair model. The CO2 adsorption energies of the chain model are in the range of 30-96 kJ/mol, with a strong metal dependence. Among these, the en-Sc2(dobpdc) and en-Cu2(dobpdc) have the highest and lowest CO2 adsorption energies, respectively. Moreover, the dynamic progress of CO2 adsorption has been unveiled via exploration of the full reaction pathway, including transition states and intermediates. First, the CO2 molecule interacts with en-MOFs to form a physisorbed complex with a shallow potential well. This is followed by overcoming a relatively large energy barrier to form a chemisorbed complex. Finally, ammonium carbamate is formed along the one-dimensional channels within the pore with a small energy barrier for configuration transformation. These results agree well with the experimental observations. Understanding the detailed microscopic mechanism of CO2 capture is quite crucial for improving our fundamental knowledge base and potential future applications. This work will improve our understanding of CO2 adsorption with amine functionalized MOFs. We expect our results to stimulate future experimental and theoretical research and advance the development of this field.
Bibliographical noteFunding Information:
H.Z., C.L.S., and L.-M.Y. gratefully acknowledge support from the National Natural Science Foundation of China (21673087, 21873032, 22073033, 21903032), startup fund (2006013118 and 3004013105) from Huazhong University of Science and Technology, and the Fundamental Research Funds for the Central Universities (2019kfyRCPY116). The authors thank the Minnesota Supercomputing Institute (MSI) at the University of Minnesota for supercomputing resources. The calculations were carried out on the Mesabi supercomputer. Additional work was carried out at the LvLiang Cloud Computing Center of China, and some calculations were performed on TianHe-2.