Why Porous Materials Have Selective Adsorptions: A Rational Aspect

Feb 15, 2017 - Gas storage/separation is a typical application of porous materials such as metal organic frameworks (MOFs). The adsorption/separation ...
0 downloads 7 Views 3MB Size
Article pubs.acs.org/IC

Why Porous Materials Have Selective Adsorptions: A Rational Aspect from Electrodynamics Qiang Chen,†,‡,§ Yuguang Ma,∥ Wei-Chao Song,‡ Ze Chang,†,§ Jian-Rong Li,⊥ Jun Zhang,‡ Hong-Wei Sun,‡ Perla B Balbuena,∥ and Xian-He Bu*,†,‡,§ †

School of Materials Science and Engineering, National Institute for Advanced Materialsand, TKL of Metal and Molecule-Based Material Chemistry, Nankai University, Tianjin 300350, China ‡ College of Chemistry and §Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Nankai University, Tianjin 300071, China ∥ Department of Chemical Engineering, Texas A&M University, College Station, Texas 77843, United States ⊥ College of Environmental and Energy Engineering, Beijing University of Technology, Beijing, 100124, China S Supporting Information *

ABSTRACT: Gas storage/separation is a typical application of porous materials such as metal organic frameworks (MOFs). The adsorption/ separation behavior results from the host−guest and/or guest−guest interaction and equilibration (host, porous material; guest, adsorbates). Although the driving forces for gas adsorption have been investigated, a detailed picture of interactions between gas molecules and MOFs has not clearly emerged. Herein, a new cobalt microporous MOF [Co(tipb)(adc)](DMF)3(H2O)1.5, which possesses a rare self-interpenetrated gra topology, has been prepared with both tipb and H2adc ligands (tipb = 1,3,5-tris(p-imidazolylphenyl)benzene, adc = 9,10anthracenedicarboxylate). This MOF shows high stability and exceptional selective adsorption of CO2 over N2, O2, and CH4. In particular, a theoretical assumption of a “regional dynamic electric field effect” is proposed to clarify the selective adsorption. Moreover, we suggest that the proposed effect may be one of the most important factors impacting gas separation and storage in porous materials.



INTRODUCTION In view of energetic, biological, and environmental factors, gas separations and storage play significant roles in industry and daily life.1,2 In particular, as one kind of greenhouse gas and a renewable carbon source, CO2 has attracted much attention. The separation of CO2 over N2/O2 for its enrichment and recovery as well as that over CH4 for biogas upgrade have motivated strong interest in the development of novel porous materials3−6 such as microporous metal−organic frameworks (MOFs). Their potential is due not only to their fascinating structural architectures and topologies but also to their tailorable structures, properties, and associated potential applications.7−11 As good candidates for gas separations and storage, the driving force for gas adsorption in MOFs has been investigated and various strategies have been developed to improve their working capacity and selectivity.12−16 The adsorption/separation behavior results from the host−guest and/or guest−guest interaction and equilibration (host, porous material; guest, adsorbates). However, a clear picture of the interactions between gas molecules and the framework has not emerged, and some explanations of selective gas adsorption in MOFs are somewhat ambiguous, except for the pore size exclusion effect17,18 or quadrupole moment19−22 of gas © XXXX American Chemical Society

molecules (see the Supporting Information for more details), particularly in cases where special adsorption behaviors were observed. Herein, a new CoII-based MOF [Co(tipb)(adc)](DMF)3(H2O)1.5 (1, tipb = 1,3,5-tris(p-imidazolylphenyl)benzene; adc = 9,10-anthracenedicarboxylate) is presented where a high adsorptive selectivity of CO2 over N2, O2, and CH4 was observed at different temperatures. Most importantly, the selective adsorption could be reasonably well explained by the presence of the “regional dynamic electric field effect” in the framework. Furthermore, as electric fields commonly exist in all porous materials, this effect can be considered as a universal/ common influence factor for gas adsorption in porous materials.



RESULTS AND DISCUSSION Syntheses and Structural Characterization of 1. The solvothermal reaction of Co(NO3)2·6H2O, tipb, and H2adc in N,N′-dimethylformamide (DMF) resulted in the formation of red prismatic crystals of 1. X-ray crystallographic analysis revealed that only one crystallographically independent CoII Received: November 20, 2016

A

DOI: 10.1021/acs.inorgchem.6b02746 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

(Figure S4).24 Additionally, a stability test confirmed that 1 has high thermal and chemical stability (see the Supporting Information for more details). Selective Gas Adsorptions. In order to characterize the nature of the pores and gas adsorption properties of 1, isotherms are measured for N2, H2, O2, CH4, and CO2. The N2 adsorption isotherm of activated 1 is measured at 77 K (Figure S10). It shows a classic type I isotherm with a Brunauer− Emmett−Teller (BET) surface area of 514 m2 g−1 and Langmuir surface area of 680 m2 g−1. The fit of the adsorption data to the Horvath−Kawazoe method demonstrated the pore width distribution from about 0.5 to 0.9 nm (Figure S11). At 195 K, 1 exhibits a type I isotherm toward CO2 gas and gives an uptake of 28.51 wt % (145 cm3 g−1 at STP) at 620 mmHg. Most interestingly, 1 shows high adsorption selectivity for CO2 over N2 and O2 as shown in Figure 2. The adsorption selectivities are calculated via ideal adsorption solution theory (IAST).25 The CO2/N2 (in a 15:85 molar ratio) selectivity is estimated as 2094 at 195 K (80 KPa) and 27.8 at 273 K (100 KPa). The predicted CO2/CH4 selectivity is 113 at 195 K (80 KPa), 7.4 at 273 K (100 KPa), and 5.4 at 303 K (100 KPa) from equimolar gas-phase mixtures (Figures S14−16). At 303 K, almost no N2 and O2 is adsorbed (Figure S13). It is insignificant to calculate the selectivity. Assumption of a Regional Dynamic Electric Field Effect. Compound 1 indeed has shown remarkable adsorptive selectivity for CO2, and the selective adsorption orders are CO2 > CH4 > N2/O2 at different temperatures. However, the origin of the selectivity could not be attributed to a size/shape-based effect since the kinetic diameters of CO2, CH4, N2, and O2 are all much smaller than the diameter of the channel of 1. By checking carefully the structure of 1, the selective adsorption is assumed to originate from the interactions between the gas molecules and the local electric fields of the framework: Commonly, electronegative atoms or groups tend to distribute around electropositive metal centers symmetrically due to the electrostatic repulsion of the homocharges, forming symmetrical structures such as the classical paddlewheel26,27 secondary building unit (SBU). However, in 1 the steric hindrance of coordinated imidazole rings may compel the electronegative carboxyl groups to locate on one side of the CoII atoms (Figure S5). This results in a separation of positive and negative charge centers in this region and accordingly leads to regional electric fields around the CoII atoms (Figure 3a). This phenomenon can be defined as “unequal coordination”. As shown in Figure 3b, the CoII centers are quite close to the wall of the 1D channel, and the distance between the nearest

center is involved in the fundamental asymmetric unit. CoII locates in the center of a distorted octahedral coordination geometry, defined by three nitrogen atoms from three independent tipb ligands and three oxygen atoms from two different adc2− ligands (Figure S1). Each tipb ligand connects three CoII atoms, affording a 63 layer parallel to the ab plane (Figure S2). The 63 layers packed in an ABCD fashion were interconnected by adc2− ligands to give an overall threedimensional (3D) framework with one-dimensional (1D) channels (0.96 × 0.73 nm in size) along the c axis (Figure 1). After removing the DMF and H2O molecules from the

Figure 1. (a) Schematic diagram of the proposed interactions between regional electric fields and gas molecules. 1D channel of 1 viewed along the (b) c and (c) a directions. Adjacent 63 layers are presented in different colors to show the ABCD packing fashion of single networks.

channels, PLATON23 calculations show that the porous structure of 1 composes a void of 3706.0 Å3 that represents 35% per unit cell volume [10586.0 Å3]. Taking into account the vdW surface, the free volume is 5352.4 Å3, occupying ∼51% volume of a unit cell. Topology analysis (see the Supporting Information for more details) revealed that the framework of 1 could be simplified into a rare self-interpenetrated gra net

Figure 2. Adsorption isotherms of 1 for CO2, CH4, N2, and O2 at (a) 195 and (b) 273 K (filled symbols, adsorption; open symbols, desorption). B

DOI: 10.1021/acs.inorgchem.6b02746 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

Figure 3. (a) Schematic diagram of the electric field formed by the permanent dipole near CoII center. (b) Approximate scope of the regional electric fields in the channel of 1.

Figure 4. (a) Calculated electrostatic potential on an isodensity surface (ρ = 0.005) near the CoII center. (b) Calculated electrostatic potential contours on the (001) plane.

For gas molecules diffusing inside the regional electric field, the interaction between a gas molecule and the electric field can be expressed by the multipole expansion of the energy of a charge distribution28 (Figure 5)

CoII is 8.4 Å. Therefore, a cutoff of the regional electric field was set to be ∼8.4 Å (the minimum diameter, Figure 3b). In this region, the electric field should be dominated by the permanent dipole located around the metal center. The strength of the permanent dipole moment is evaluated using a cluster model (see the Supporting Information for more details). The first-principles calculation indicates a large dipole moment (Figure S19) of ∼12.7 D. This large dipole moment suggests a strong electric field around the CoII centers. To confirm the existence of the electric field, the electrostatic potential of the cluster model was calculated by density functional theory (DFT) (Figure 4a). The distribution of the electrostatic potential is not even. The carboxyl O atoms exhibit strong electronegativity, while the aromatic H atoms are slightly electrically positive. As illustrated in Figure 4b, the contours of electrostatic potential indicate the direction and strength of the electric field, which varies locally according to the position in the plane. These results suggest that that there exists a nonuniform electric field around the CoII center defined as a regional electric field. The term “regional” refers to the local environment of the porous material and to the subsequent effect on the electric field characteristics. It should be pointed out that the permanent electrical dipole near CoII is the most obvious source of the regional electric field, but it is not the only one. From the calculated results (Figure 4), the imidazole rings and anthracene nucleus are also important for the formation of the regional electric fields in 1. This indicates that the electric fields may be originated from a variety of effects. The electric field and the potential gradient play important roles in gas adsorption in 1.

W = qΦ(0) − p ·E(0) −

1 6

∑ ∑ Q ij i

j

∂Ej ∂Xi

(0) + ···

In the equation, q is the total charge while p and Qij are the electric dipole moment and the quadrupole moment tensor components of a gas molecule. In our study, the charge of gas molecules is zero, and the effect of other higher order multipoles can be neglected compared with the dipole/ quadrupole moments. Therefore, the major contribution comes from the interaction of the dipole/quadrupole of the gas molecules and the regional electric field. Moreover, for nonpolar gas molecules, their permanent electric dipole moments are equal to zero. However, in an external electric field, nonpolar gas molecules can be polarized to give induced dipole moments, as well as instantaneous dipole moments. These two types of dipoles are related to the polarizability of the gas molecules. They are also of importance in the MOF− gas interactions. As both the polarizability and the quadrupole moment of CO2 are larger than those of CH4, N2, and O2 (Table S1), CO2 should have stronger interaction with the regional electric field in the 1D channels of 1. This can explain the order of adsorption selectivity for CO2 over CH4/N2/O2. On the other hand, although the quadrupole moment of CH4 is zero, the polarizability of CH4 is far larger than that of N2 and O2. As a C

DOI: 10.1021/acs.inorgchem.6b02746 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

adsorbed adsorbates will interact with adsorbates which are going to be adsorbed (van der Waals force). Consequently, electric fields may exist extensively in MOFs, and the fields are dynamic. The current research of multipoint interactions also substantiated the origins of electric fields are various.41,42 The interactions between the framework and the gas molecules can be described as a regional dynamic electric field effect. The term “dynamic” refers to the electric fields of effective adsorption space that are variable for the polarity of the guests and/or other influencing factors. The assumption can also be extensively applied to explain the adsorption of other polar/nonpolar molecules. For polar molecules, the permanent dipole will play a leading role in the interaction between electric field and polar molecules. Some reported results of MOF materials in HKUST-1 and Ni/ DOBDC corroborate that H2O is more easily adsorbed than CO2.43 This may be attributed to the strong permanent dipole moment of H2O (about 1.85 × 1018 esu·cm). While, for gas molecules with a weak permanent dipole such as CO, the effect of instantaneous dipole, induced dipole, and quadrupole moment shall be the crucial factors. For example, though the permanent dipole of CO2 is zero, both its polarizability and its quadrupole moment are higher than that of CO. Considering all factors, CO2 shall have a stronger interaction with MOFs than that of CO, which has been proved by previous studies.44,45 Besides, there are no significant correlations between the selective adsorption order and the ionization energy of the adsorbates (Table S1). Thus, the effect of induced dipole moments is far greater than instantaneous dipole moments in our research system. Besides MOFs materials, the assumption can be used to explain the adsorption behavior in other porous materials. For instance, partially charged single-walled carbon nanotubes by methylene blue encapsulation can enhance the CO2 absorptivity;46 Zhou et al. grafted a porous polymer network with polar sulfonic acid (PPN-6-SO3H) and its lithium salt (PPN-6SO3Li), which exhibit significant increases in isosteric heats of CO2 adsorption and CO2-uptake capacities;47 Cooper et al. reported a series of conjugated microporous polymer (CMP) networks.48 The carboxylic acid-functionalized network shows the highest isosteric heat of sorption for CO2. All these results could be rationally attributed to the presence of a regional dynamic electric field in these materials and explained with the assumption.

Figure 5. Outline of the theoretical assumption “regional dynamic electric field effect”.

result, the combined interaction of the gas molecules and the framework has an order of CO2 > CH4 > N2/O2. This is verified by periodic DFT calculations (see the Supporting Information for more details): The calculated interaction energy is −4.2 kcal/mol for N2, −6.2 kcal/mol for CH4, and −7.8 kcal/mol for CO2. The calculated values are in good agreement with the experimental results and can be attributed to the effect of the regional electric field. Notably, the dipole moments play more important roles than the quadrupole moments in this case. In summary, the presence of a regional electric field can reasonably explain the selective adsorption order of CO2 > CH4 > N2/O2 in 1. It should be pointed out that the explanation above does not involve specific porous materials. It can be expressed as a universal hypothesis, which may apply to (Figure 5) other MOFs and porous materials as discussed next. As shown in 1, the electric fields originated from a variety of effects. It mainly involves two aspects (Figure 5). One aspect is the host framework of porous materials: as we discussed, the regional electric field can be originated from the unequal coordination. Also, it may be generated by unsaturated metal centers29−31 in some MOFs. Moreover, it can be produced by polar ligands and/or groups32−35 in porous materials: some investigations have confirmed that polar groups such as acylamide and carboxyl groups could enhance the CO2 uptake and the corresponding heat of adsorption. In addition, there are some other triggers such as ionic frameworks,36−38 delocalization of π electrons, instantaneous dipole moment of the groups in the framework, etc. Another aspect is guest molecules or/and adsorbates. The adsorbed/encapsulated polar molecules or ions in MOFs are important sources of electric fields: for example, recent studies have confirmed that small amounts of water can help to enhance the adsorption of CO2 in HKUST-1.39,40 The



CONCLUSIONS In conclusion, we successfully synthesized a new MOF which presents an interesting self-interpenetrated gra topology and exceptional adsorption selectivity for CO2. It is found that the selectivity is closely related to the effect of a regional dynamic electric field present in the framework. The factors that influence physical adsorption were discussed in detail from the point of electrodynamics: The regional dynamic electric fields exist extensively in MOFs and other porous materials, and the intrinsic properties of adsorbates (such as polarizability, quadruple moment) are the crucial factors that influence gas adsorption. According to the regional dynamic electric field effect, the adsorption capacity and selectivity could be enhanced by heightening the electric field intensity and potential gradient of the effective adsorption space of porous materials through structure modification, which is the central advantage of MOF materials. D

DOI: 10.1021/acs.inorgchem.6b02746 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry



EXPERIMENTAL SECTION

Xian-He Bu: 0000-0002-2646-7974

Synthesis of [Co(tipb)(adc)](DMF)3(H2O)1.5 (1). A mixture of Co(NO3)·6H2O (29.1 mg, 0.1 mmol), tipb (50.4 mg, 0.1 mmol), H2adc (27.0 mg, 0.1 mmol), and 2 drops of pyridine in 5 mL of DMF was sealed in a 16 mL Teflon-lined stainless steel container and heated at 160 °C for 3 days. The container was then cooled to room temperature at a rate of 0.1 °C/min, and the red crystalline product was separated by filtration, washed with DMF, and dried in air. Yield: 62% (based on tipb). Despite extensive efforts, elemental analyses are not well reproduced, presumably because the solvent molecules in the frameworks are variable. IR (KBr pellets, cm−1): 3134, 1664, 1598, 1521, 1400, 1309, 1061, 963, 833, 656. Calculation of Electrical Dipole Moment. The cluster model of 1 was extracted from the periodic framework for evaluation of the strength of the dipole moment. Taking into account the entire conjugation effect, 1 includes all the nearly coplanar molecular fragments. A single-point DFT calculation was performed using the B3LYP functional49 (Becke three-parameter exchange hybrid functional combined with the Lee−Yang−Parr correlation functional) with the 6-31g(d) basis set for C, H, O, and N atoms and the Stuttgart RSC 1997 effective core pseudopotential50 for Co atoms. The result shows that the centers of positive electrical charge and negative charge are near CoII. The direction of the dipole moment is from the oxygen atoms to the cobalt atom (Figure S19). Notably, the magnitude of the dipole moment is as large as ∼12.7 D, which is rarely reported in other MOF systems. Calculation of Binding Energy. Periodic density functional theory (DFT) calculations were performed using the Vienna ab initio Simulation packages,51,52 in which the Kohn−Sham equations are solved by self-consistent algorithms. The exchange-correlation functional used in calculations was in the framework of the generalized gradient approximation (GGA) proposed by Perdew and Wang (PW91).53,54 Valence electrons were described by a plane wave basis set with a cutoff energy of 400 eV, while inner core electrons were described by using Vanderbilt’s ultrasoft pseudopotentials for computational efficiency.55 A 2 × 2 × 2 Monkhorst−Pack k-point mesh was used for Brillouin zone sampling in the reciprocal space.56 A conjugation gradient algorithm was applied to ionic relaxation and until the energy difference was less than 10−4 eV. An semiempirical dispersion potential57 was added to the conventional Kohm−Sham energy for evaluation of van der Waals (vdW) interactions in the DFT calculations. The framework atoms were fixed to their positions in the crystal structure, whereas the gas molecules were relaxed to reach a local minimum energy configuration. The binding energy is calculated by the following equation: Ebind = Efrwk/gas − Efrwk − Egas, in which Efwk/gas is the total energy of the framework and the adsorbed gas molecule(s), Efrwk is the energy of the framework, and Egas is the energy of the gas molecule(s).



Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the 973 Program of China (2014CB845601), National Science Foundation of China (21531005, 21421001, and 21202088), MOE Innovation Team of China (IRT13022), and China Postdoctoral Science Foundation (2015M571260). We thank Prof. Michael O’Keeffe for the advice on the topology analysis and Prof. Susumu Kitagawa and Prof. Michael J. Zaworotko for the helpful discussion.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.inorgchem.6b02746. Expermental details, crystallographic data, supplementary figures, IR, TGA, PXRD, adsorption data, calculation and modeling (PDF) Crystallographic data for 1 (CIF)



REFERENCES

(1) Li, J. R.; Sculley, J.; Zhou, H. C. Metal-organic frameworks for separations. Chem. Rev. 2012, 112, 869−932. (2) Van de Voorde, B.; Bueken, B.; Denayer, J.; De Vos, D. Adsorptive separation on metal-organic frameworks in the liquid phase. Chem. Soc. Rev. 2014, 43, 5766−5788. (3) Datta, S. J.; Khumnoon, C.; Lee, Z. H.; Moon, W. K.; Docao, S.; Nguyen, T. H.; Hwang, I. C.; Moon, D.; Oleynikov, P.; Terasaki, O.; Yoon, K. B. CO2 capture from humid flue gases and humid atmosphere using a microporous coppersilicate. Science 2015, 350, 302−306. (4) Foo, M. L.; Matsuda, R.; Hijikata, Y.; Krishna, R.; Sato, H.; Horike, S.; Hori, A.; Duan, J.; Sato, Y.; Kubota, Y.; Takata, M.; Kitagawa, S. An Adsorbate Discriminatory Gate Effect in a Flexible Porous Coordination Polymer for Selective Adsorption of CO2 over C2H2. J. Am. Chem. Soc. 2016, 138, 3022−3030. (5) Keith, D. W. Why Capture CO2 from the atmosphere? Science 2009, 325, 1654−1655. (6) Yang, Q.; Vaesen, S.; Ragon, F.; Wiersum, A. D.; Wu, D.; Lago, A.; Devic, T.; Martineau, C.; Taulelle, F.; Llewellyn, P. L.; Jobic, H.; Zhong, C.; Serre, C.; De Weireld, G.; Maurin, G. A Water Stable Metal-Organic Framework with Optimal Features for CO2 Capture. Angew. Chem., Int. Ed. 2013, 52, 10316−10320. (7) Gandara, F.; Furukawa, H.; Lee, S.; Yaghi, O. M. High Methane Storage Capacity in Aluminum Metal-Organic Frameworks. J. Am. Chem. Soc. 2014, 136, 5271−5274. (8) Luo, F.; Yan, C.; Dang, L.; Krishna, R.; Zhou, W.; Wu, H.; Dong, X.; Han, Y.; Hu, T.-L.; O’Keeffe, M.; Wang, L.; Luo, M.; Lin, R.-B.; Chen, B. UTSA-74: A MOF-74 Isomer with Two Accessible Binding Sites per Metal Center for Highly Selective Gas Separation. J. Am. Chem. Soc. 2016, 138, 5678−5684. (9) Makiura, R.; Motoyama, S.; Umemura, Y.; Yamanaka, H.; Sakata, O.; Kitagawa, H. Surface nano-architecture of a metal-organic framework. Nat. Mater. 2010, 9, 565−571. (10) Sumida, K.; Stueck, D.; Mino, L.; Chai, J.-D.; Bloch, E. D.; Zavorotynska, O.; Murray, L. J.; Dinca, M.; Chavan, S.; Bordiga, S.; Head-Gordon, M.; Long, J. R. Impact of Metal and Anion Substitutions on the Hydrogen Storage Properties of M-BTT MetalOrganic Frameworks. J. Am. Chem. Soc. 2013, 135, 1083−1091. (11) Zhang, H.; Deria, P.; Farha, O. K.; Hupp, J. T.; Snurr, R. Q. A thermodynamic tank model for studying the effect of higher hydrocarbons on natural gas storage in metal-organic frameworks. Energy Environ. Sci. 2015, 8, 1501−1510. (12) Deshmukh, M. M.; Ohba, M.; Kitagawa, S.; Sakaki, S. Absorption of CO2 and CS2 into the Hofmann-Type Porous Coordination Polymer: Electrostatic versus Dispersion Interactions. J. Am. Chem. Soc. 2013, 135, 4840−4849. (13) Li, J. R.; Kuppler, R. J.; Zhou, H. C. Selective gas adsorption and separation in metal-organic frameworks. Chem. Soc. Rev. 2009, 38, 1477−504. (14) Lin, J.-B.; Zhang, J.-P.; Chen, X.-M. Nonclassical Active Site for Enhanced Gas Sorption in Porous Coordination Polymer. J. Am. Chem. Soc. 2010, 132, 6654−6656.

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Ze Chang: 0000-0002-3865-2165 Jian-Rong Li: 0000-0002-8101-8493 Perla B Balbuena: 0000-0002-2358-3910 E

DOI: 10.1021/acs.inorgchem.6b02746 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry (15) Nugent, P.; Belmabkhout, Y.; Burd, S. D.; Cairns, A. J.; Luebke, R.; Forrest, K.; Pham, T.; Ma, S.; Space, B.; Wojtas, L.; Eddaoudi, M.; Zaworotko, M. J. Porous materials with optimal adsorption thermodynamics and kinetics for CO2 separation. Nature 2013, 495, 80−84. (16) Pham, T.; Forrest, K. A.; McLaughlin, K.; Tudor, B.; Nugent, P.; Hogan, A.; Mullen, A.; Cioce, C. R.; Zaworotko, M. J.; Space, B. Theoretical Investigations of CO2 and H2 Sorption in an Interpenetrated Square-Pillared Metal−Organic Material. J. Phys. Chem. C 2013, 117, 9970−9982. (17) Thallapally, P. K.; Tian, J.; Radha Kishan, M.; Fernandez, C. A.; Dalgarno, S. J.; McGrail, P. B.; Warren, J. E.; Atwood, J. L. Flexible (Breathing) Interpenetrated Metal−Organic Frameworks for CO2 Separation Applications. J. Am. Chem. Soc. 2008, 130, 16842−16843. (18) Xue, D.-X.; Belmabkhout, Y.; Shekhah, O.; Jiang, H.; Adil, K.; Cairns, A. J.; Eddaoudi, M. Tunable Rare Earth fcu-MOF Platform: Access to Adsorption Kinetics Driven Gas/Vapor Separations via Pore Size Contraction. J. Am. Chem. Soc. 2015, 137, 5034−5040. (19) Herm, Z. R.; Swisher, J. A.; Smit, B.; Krishna, R.; Long, J. R. Metal−Organic Frameworks as Adsorbents for Hydrogen Purification and Precombustion Carbon Dioxide Capture. J. Am. Chem. Soc. 2011, 133, 5664−5667. (20) Tsivion, E.; Long, J. R.; Head-Gordon, M. Hydrogen Physisorption on Metal-Organic Framework Linkers and Metalated Linkers: A Computational Study of the Factors That Control Binding Strength. J. Am. Chem. Soc. 2014, 136, 17827−17835. (21) Vaidhyanathan, R.; Iremonger, S. S.; Shimizu, G. K. H.; Boyd, P. G.; Alavi, S.; Woo, T. K. Direct Observation and Quantification of CO2 Binding Within an Amine-Functionalized Nanoporous Solid. Science 2010, 330, 650−653. (22) Wang, Q.; Bai, J.; Lu, Z.; Pan, Y.; You, X. Finely tuning MOFs towards high-performance post-combustion CO2 capture materials. Chem. Commun. 2016, 52, 443−452. (23) Spek, A. L. Single-crystal structure validation with the program PLATON. J. Appl. Crystallogr. 2003, 36, 7−13. (24) O’Keeffe, M.; Peskov, M. A.; Ramsden, S. J.; Yaghi, O. M. The Reticular Chemistry Structure Resource (RCSR) Database of, and Symbols for, Crystal Nets. Acc. Chem. Res. 2008, 41, 1782−1789. (25) Gao, S.; Morris, C. G.; Lu, Z.; Yan, Y.; Godfrey, H. G. W.; Murray, C.; Tang, C. C.; Thomas, K. M.; Yang, S.; Schroder, M. Selective Hysteretic Sorption of Light Hydrocarbons in a Flexible Metal-Organic Framework Material. Chem. Mater. 2016, 28, 2331− 2340. (26) Tranchemontagne, D. J.; Mendoza-Cortes, J. L.; O’Keeffe, M.; Yaghi, O. M. Secondary building units, nets and bonding in the chemistry of metal-organic frameworks. Chem. Soc. Rev. 2009, 38, 1257−83. (27) Yan, Y.; Juricek, M.; Coudert, F.-X.; Vermeulen, N. A.; Grunder, S.; Dailly, A.; Lewis, W.; Blake, A. J.; Stoddart, J. F.; Schroeder, M. Non-Interpenetrated Metal-Organic Frameworks Based on Copper(II) Paddlewheel and Oligoparaxylene-lsophthalate Linkers: Synthesis, Structure, and Gas Adsorption. J. Am. Chem. Soc. 2016, 138, 3371− 3381. (28) Jackson, J. D. Classical Electrodynamics; Wiley: New York, 1962. (29) Bae, Y.-S.; Spokoyny, A. M.; Farha, O. K.; Snurr, R. Q.; Hupp, J. T.; Mirkin, C. A. Separation of gas mixtures using Co(ii) carboranebased porous coordination polymers. Chem. Commun. 2010, 46, 3478−3480. (30) Forrest, K. A.; Pham, T.; McLaughlin, K.; Belof, J. L.; Stern, A. C.; Zaworotko, M. J.; Space, B. Simulation of the Mechanism of Gas Sorption in a Metal-Organic Framework with Open Metal Sites: Molecular Hydrogen in PCN-61. J. Phys. Chem. C 2012, 116, 15538− 15549. (31) Sumida, K.; Horike, S.; Kaye, S. S.; Herm, Z. R.; Queen, W. L.; Brown, C. M.; Grandjean, F.; Long, G. J.; Dailly, A.; Long, J. R. Hydrogen storage and carbon dioxide capture in an iron-based sodalite-type metal-organic framework (Fe-BTT) discovered via highthroughput methods. Chem. Sci. 2010, 1, 184−191.

(32) Yang, Q.; Wiersum, A. D.; Llewellyn, P. L.; Guillerm, V.; Serre, C.; Maurin, G. Functionalizing porous zirconium terephthalate UiO66(Zr) for natural gas upgrading: a computational exploration. Chem. Commun. 2011, 47, 9603−9605. (33) Zheng, B.; Bai, J.; Duan, J.; Wojtas, L.; Zaworotko, M. J. Enhanced CO2 Binding Affinity of a High-Uptake rht-Type Metal− Organic Framework Decorated with Acylamide Groups. J. Am. Chem. Soc. 2011, 133, 748−751. (34) Noro, S.-i.; Tanaka, D.; Sakamoto, H.; Shimomura, S.; Kitagawa, S.; Takeda, S.; Uemura, K.; Kita, H.; Akutagawa, T.; Nakamura, T. Selective Gas Adsorption in One-Dimensional, Flexible CuII Coordination Polymers with Polar Units. Chem. Mater. 2009, 21, 3346−3355. (35) Sun, J.-K.; Ji, M.; Chen, C.; Wang, W.-G.; Wang, P.; Chen, R.-P.; Zhang, J. A charge-polarized porous metal-organic framework for gas chromatographic separation of alcohols from water. Chem. Commun. 2013, 49, 1624−1626. (36) Nouar, F.; Eckert, J.; Eubank, J. F.; Forster, P.; Eddaoudi, M. Zeolite-like Metal-Organic Frameworks (ZMOFs) as Hydrogen Storage Platform: Lithium and Magnesium Ion-Exchange and H2(rho-ZMOF) Interaction Studies. J. Am. Chem. Soc. 2009, 131, 2864− 2870. (37) Yang, S.; Lin, X.; Blake, A. J.; Walker, G. S.; Hubberstey, P.; Champness, N. R.; Schroeder, M. Cation-induced kinetic trapping and enhanced hydrogen adsorption in a modulated anionic metal-organic framework. Nat. Chem. 2009, 1, 487−493. (38) Wang, S.; Yang, Q.; Zhang, J.; Zhang, X.; Zhao, C.; Jiang, L.; Su, C.-Y. Two-Dimensional Charge-Separated Metal−Organic Framework for Hysteretic and Modulated Sorption. Inorg. Chem. 2013, 52, 4198− 4204. (39) Chen, Y. F.; Babarao, R.; Sandler, S. I.; Jiang, J. W. Metal− Organic Framework MIL-101 for Adsorption and Effect of Terminal Water Molecules: From Quantum Mechanics to Molecular Simulation. Langmuir 2010, 26, 8743−8750. (40) Yazaydın, A. O. z. r.; Benin, A. I.; Faheem, S. A.; Jakubczak, P.; Low, J. J.; Willis, R. R.; Snurr, R. Q. Enhanced CO2 Adsorption in Metal-Organic Frameworks via Occupation of Open-Metal Sites by Coordinated Water Molecules. Chem. Mater. 2009, 21, 1425−1430. (41) Cui, P.; Ma, Y.-G.; Li, H.-H.; Zhao, B.; Li, J.-R.; Cheng, P.; Balbuena, P. B.; Zhou, H.-C. Multipoint Interactions Enhanced CO2 Uptake: A Zeolite-like Zinc−Tetrazole Framework with 24-Nuclear Zinc Cages. J. Am. Chem. Soc. 2012, 134, 18892−18895. (42) He, W.-W.; Li, S.-L.; Li, W.-L.; Li, J.-S.; Yang, G.-S.; Zhang, S.R.; Lan, Y.-Q.; Shen, P.; Su, Z.-M. Multipoint interactions enhanced H2 storage and organosulfur removal in a microporous metal-organic framework. J. Mater. Chem. A 2013, 1, 11111−11116. (43) Liu, J.; Wang, Y.; Benin, A. I.; Jakubczak, P.; Willis, R. R.; LeVan, M. D. CO2/H2O Adsorption Equilibrium and Rates on Metal− Organic Frameworks: HKUST-1 and Ni/DOBDC. Langmuir 2010, 26, 14301−14307. (44) Deng, H.; Doonan, C. J.; Furukawa, H.; Ferreira, R. B.; Towne, J.; Knobler, C. B.; Wang, B.; Yaghi, O. M. Multiple Functional Groups of Varying Ratios in Metal-Organic Frameworks. Science 2010, 327, 846−850. (45) Eguchi, R.; Uchida, S.; Mizuno, N. Inverse and High CO2/C2H2 Sorption Selectivity in Flexible Organic−Inorganic Ionic Crystals. Angew. Chem., Int. Ed. 2012, 51, 1635−1639. (46) Khoerunnisa, F.; Fujimori, T.; Itoh, T.; Urita, K.; Hayashi, T.; Kanoh, H.; Ohba, T.; Hong, S. Y.; Choi, Y. C.; Santosa, S. J.; Endo, M.; Kaneko, K. Enhanced CO2 Adsorptivity of Partially Charged Single Walled Carbon Nanotubes by Methylene Blue Encapsulation. J. Phys. Chem. C 2012, 116, 11216−11222. (47) Luo, Y. L.; Li, B. Y.; Liang, L. Y.; Tan, B. E. Synthesis of costeffective porous polyimides and their gas storage properties. Chem. Commun. 2011, 47, 7704−7706. (48) Dawson, R.; Adams, D. J.; Cooper, A. I. Chemical tuning of CO2 sorption in robust nanoporous organic polymers. Chem. Sci. 2011, 2, 1173−1177. F

DOI: 10.1021/acs.inorgchem.6b02746 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry (49) Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648−5652. (50) Dolg, M.; Stoll, H.; Preuss, H.; Pitzer, R. M. Relativistic and correlation effects for element 105 (hahnium, Ha): a comparative study of M and MO (M = Nb, Ta, Ha) using energy-adjusted ab initio pseudopotentials. J. Phys. Chem. 1993, 97, 5852−5859. (51) Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (52) Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15−50. (53) Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M. R.; Singh, D. J.; Fiolhais, C. Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B: Condens. Matter Mater. Phys. 1992, 46, 6671−6687. (54) Perdew, J. P.; Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B: Condens. Matter Mater. Phys. 1992, 45, 13244−13249. (55) Vanderbilt, D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys. Rev. B: Condens. Matter Mater. Phys. 1990, 41, 7892−7895. (56) Monkhorst, H. J.; Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188−5192. (57) Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787−1799.

G

DOI: 10.1021/acs.inorgchem.6b02746 Inorg. Chem. XXXX, XXX, XXX−XXX