Investigating the gas sorption mechanism in an rht -metal-organic framework through computational studies
AuthorsPham, Tony T.
Forrest, Katherine A.
Georgiev, Peter A L
Mullen, Ashley L.
Eubank, Jarrod F.
KAUST DepartmentAdvanced Membranes and Porous Materials Research Center
Physical Sciences and Engineering (PSE) Division
Chemical Science Program
Functional Materials Design, Discovery and Development (FMD3)
MetadataShow full item record
AbstractGrand canonical Monte Carlo (GCMC) simulations were performed to investigate CO2 and H2 sorption in an rht-metal-organic framework (MOF) that was synthesized with a ligand having a nitrogen-rich trigonal core through trisubstituted triazine groups and amine functional groups. This MOF was synthesized by two different groups, each reporting their own distinct gas sorption measurements and crystal structure. Electronic structure calculations demonstrated that the small differences in the atomic positions between each group's crystal structure resulted in different electrostatic parameters about the Cu2+ ions for the respective unit cells. Simulations of CO2 sorption were performed with and without many-body polarization effects and using our recently developed CO2 potentials, in addition to a well-known bulk CO2 model, in both crystallographic unit cells. Simulated CO2 sorption isotherms and calculated isosteric heats of adsorption, Qst, values were in excellent agreement with the results reported previously by Eddaoudi et al. for both structures using the polarizable CO2 potential. For both crystal structures, the initial site for CO2 sorption were the Cu 2+ ions that had the higher positive charge in the unit cell, although the identity of this electropositive Cu2+ ion was different in each case. Simulations of H2 sorption were performed with three different hydrogen potentials of increasing anisotropy in both crystal structures and the results, especially with the highest fidelity model, agreed well with Eddaoudi et al.'s experimental data. The preferred site of H 2 sorption at low loading was between two Cu2+ ions of neighboring paddlewheels. The calculation of the normalized hydrogen dipole distribution for the polarizable model in both crystal structures aided in the identification of four distinct sorption sites in the MOF, which is consistent to what was observed in the experimental inelastic neutron scattering (INS) spectra. Lastly, while the experimental results for the two groups are quantitatively different, the sorption mechanisms (for both crystal structures and sorbates) are broadly similar and not inconsistent with either set of experimental data; the theoretical sorption isotherms themselves resemble those by Eddaoudi et al. © 2013 American Chemical Society.
SponsorsThis work was supported by the National Science Foundation (Award No. CHE-1152362). Computations were performed under a XSEDE Grant (No. TG-DMR090028) to B.S. This publication is also based on work supported by Award No. FIC/2010/06, made by King Abdullah University of Science and Technology (KAUST). The authors also thank the Space Foundation (Basic and Applied Research) for partial support. B.S. would like to acknowledge the use of the services provided by Research Computing at the University of South Florida. In addition, this work is based in part on experiments performed on the TOFTOF instrument operated by FRM-II at the Heinz Maier-Leibnitz Zentrum (MLZ), Garching, Germany. P.A.G acknowledges support from the Project "Beyond Everest." under EU Programme REGPOT-2011-1. Lastly, the authors thank Professor Randy W. Larsen and his research group for interactive discussions on this project.
PublisherAmerican Chemical Society (ACS)