Biswajit Mohantya and
Gopal Avashthi*b
aDepartment of Chemistry, Central University of Haryana (CUH), Mahendergarh 123031, Haryana, India
bSchool of Sciences, P P Savani University, NH-8, GETCO, Near Biltech, Kosamba, 394125, Surat, Gujarat, India. E-mail: gopalavashthi@gmail.com
First published on 29th November 2022
The purification of light hydrocarbons is one of the most important chemical processes globally which consumes substantial energy. Porous materials are likely to improve the efficiency of the separation process by acting as regenerable solid adsorbents. To investigate such translational systems, the underlying mechanism of adsorption in the porous materials must be taken into account. Herein we report the adsorption and selective separation of C1–C4 hydrocarbons in the coinage metal-based macrocyclic metallocavitand Pillarplex, which exhibits excellent performance in the adsorption of CH4 at the ambient conditions with a binding energy of −17.9 kcal mol−1. In addition, the endohedral adsorption of C2–C4 hydrocarbon is impressive. The CH4, C2H4, C3H4, and 1,3-butadiene have potential uptake of 2.57, 4.26, 3.60, and 2.95 mmol g−1, respectively at ambient conditions are highest from their respective isomers. Selective separation of C1–C4 hydrocarbons is studied using ideal adsorption solution theory demonstrating its potential for one-step purification of C1–C3 hydrocarbons.
Molecular recognition supramolecular chemistry sets a new paradigm to study the adsorption and selective separation of guest molecules via the host–guest mechanism. The cucurbit[n]uril (n = 5,6,7,8,10) and pillar[n]arene (n = 5,6) are versatile macrocyclic molecules and ideal hosts in the supramolecular domain.7,8 Nau et al. studied the C1–C6 HCs within cucurbit[6]urils (CB[6]) for gas sensing and reveals that CB[6] is ideal for isobutene and cyclopentane adsorption.9 Coskun et al. studied the propane/methane separation using pillar[5]arene based conjugated microporous polymer via host–guest chemistry and observed the facile separation of propane from the natural gas mixture.10 Soft porous crystals of [2+3] imide based organic cage (NKPOC-1) and their HC adsorption phenomenon were investigated by Zhang et al. reveal the α-phase of NKOPC-1 selectively adsorbs propyne over propylene and propane at 298 K/1 bar condition.11
However, such organic hosts are restricted in large-scale employment due to their troublesome synthetic separation, purification, and post-synthetic modifications.12 Other than organic hosts, metal–organic frameworks (MOFs) are profound against guest binding using the supramolecular binding approach. The higher storage density of acetylene over ethylene has been studied by Schröder et al. using functionalized MOF by supramolecular binding and separation approach.13 Selective separation of C2–C4 HCs in the Zn2 (sdc)2(bpe) (sdc = 4,4′-stilbenedicarboxylate, bpe = 1,2-bis(4-pyridyl)ethane) MOF has been studied by Li et al., claimed that C2 HC has higher selectivity over C3 and C4.14 Although MOFs are an attractive candidate for HCs capture and separation, their stability remains questionable. The limited stability against reactive species such as water is a major concern as they can debase the MOF during the in-service events.15 The akin chemical properties, polarizability, and kinetic diameters of HCs (CH4, C2H6, C3H8, C2H2 have kinetic diameters of 3.8, 4.4, 4.9, 3.3, 4.2 Å respectively) coerce the complications for HC separations.15,16 As all the light HCs are adsorbent dependent, it is imperative to develop suitable adsorbents for effective HCs separation. Early reports on metal coordinated cavitands, macrocycles, and cages are the flawless supramolecular host with exceptional guest binding effect.17,18 The coinage metal-based organometallic host has potential in the post-synthetic modification, simple synthetic pathway, tunable cavity, tailor-made pore surface, and high stability against air and water.19 Pöthig and Altmann synthesized a well-defined tubular metallocavitand with Au(I) metal center; Pillarplex (PPX) analogous to pillar[5]arenes readily accommodate linear guests with high selectivity.20
To address this challenge, we have investigated the adsorption of small HCs potentially with low carbon content (C1–C4) within the cavity of the PPX. Density functional theory (DFT) is employed to determine the electronic structure and stability of the system. The natural bond orbital (NBO) charge transfer method is used to explore the electronic transfer within the host–guest system. The nature of the interaction between HCs and PPX is analyzed using a non-covalent index (NCI). The individual components of binding energy are computed using the energy decomposition analysis (EDA). The maximum uptake of HCs is studied using the Grand canonical Monte Carlo (GCMC) simulation at various temperatures from 77–340 K and the pressure range from 0.1–10 bar. The isosteric heat of adsorption is computed to determine the thermodynamic stability of the complex. Further, the selectivity of HCs from their counterpart mixtures is studied using ideal adsorption solution theory (IAST) for selective separation using PPX as a separator channel.
The depth of the cavity is 3.785 Å longer than the analogous pillar[5]arene (7.800 Å).20 The structural resemblance with pillar[5]arene rationale PPX as a potential host and required further investigations.
The low carbon content HCs starting from CH4 to C4H8 (C1–C4) and their isomers are studied extensively using PPX to understand the adsorption behavior. The stable HC@PPX complex is discerned by swapping the position and orientation of HCs within the PPX using the B97D2/LanL2DZ level of theory. The potential stable structure is shown in Fig. 1. Since the reports on the gold phosphine complexes were established as the potential catalyst, much attention has been gained in recent times to its structural conformations.21,22 It was found that the excellency of properties arises owing to Au⋯H–C distances which are remarkably short and even shorter than their van der Waals radii (2.86 Å) as well as hydrogen bonds.23
An investigation has been carried out comprising the Au⋯H–C interactions. In CH4@PPX, The CH4 resides closer to the Au atom with a distance of 3.033 Å. C2H2 was placed horizontally symmetrical fashion within the cavity with an Au⋯H–C distance of 2.937 Å. C2H4 and C2H6 molecules are in the horizontally inclined position with multiple Au⋯H–C interactions and a noticeable bond length of less than 3 Å is observed. In C3 HCs, Au⋯H–C interactions are impressive while in allene@PPX and C3H8@PPX some Au⋯H–C distances are found lesser than their van der Waals radii. C4 HCs are well buried inside the cavity of PPX and the majority of Au⋯H–C interactions are in the region of 2.721 to 3.054 Å which pronounce a weak hydrogen bond type interaction. In isobutylene and isobutane, the Au⋯H–C interactive distances are much lower compared to its isomers. The presence of such interactions influences the polarization effect by altering the bond length of HCs within the cavity. Trivial changes in bond lengths are observed for all the encapsulated HCs and details are given in Table S1.† Upon entrapment of HCs, the structural deformation is observed in the PPX, where the Au–Au distance of HC@PPX gets reduced to some extent. Conversely, the Au–Au bond length in allene@PPX (0.019 Å) and isobutane@PPX (0.055 Å) is higher than that of PPX as shown in Table S2.† The unusual Au–Au bond length in allene and isobutane is due to its stronger Au⋯H–C interaction and the distance is 2.727 and 2.721 Å, respectively which is lower than their van der Waal radii. Nonetheless, the N⋯H–C interactions can't be ignored which plays a crucial role in the HCs uptake via hydrogen bonding. Desiraju et al. established the N⋯H–C interaction as a hydrogen bond a weak one, largely electrostatic with a distance between 2.0–3.0 Å and it depends upon the C–H acidity and the –N basicity.24 In HC@PPX the more acidic π-character containing HCs such as allene (2.427 Å), C3H4 (2.616 Å), 1,2-butadiene (2.591 Å), C4H8 (2.547 Å), C4H10 (2.621 Å) and isobutane (2.348 Å) are shows hydrogen bonding while others are in electrostatic regime with N⋯H–C distances ca. 3.0 Å. The Au⋯H–C and N⋯H–C distances are shown in Fig. 2.
Fig. 2 The Au⋯H and N⋯H bond distances in the HC@PPX complex, The net binding energy values are in parenthesis inset in the figure. |
Notable interactions exist in the HC@PPX complex, it is obvious to prevail a structural distortion in both PPX and HCs. The deformation energy (DE) is calculated to determine the deformation of PPX upon complexation. The deformation energy (DE) is the energy difference between the total energy of the optimized guest and host molecule and energy after the conformation deformation occurred during the encapsulation process. DE of host and guest is calculated as follows,
DE (host) = E(host)optsp − E(host)opt |
DE (guest) = E(guest)optsp − E(guest)opt |
The DE of PPX is higher than the HCs given in Table 1. The DE of PPX in allene@PPX (3.7 kcal mol−1) is the highest and C4H8@PPX (0.4 kcal mol−1) is the least among the studied system. The DE of PPX, in C2H2@PPX, is 1.1 kcal mol−1 and 0.2 kcal mol−1 higher than that of C2H4@PPX and C2H6@PPX and the ease of π electron transfer from HC to the metal center is C2H2 > C2H4 > C2H6. The PPX in C3H8@PPX exhibits DE of 1.6 kcal mol−1 is 1.2, 1.7, and 2.1 kcal mol−1 lower than that of C3H6, C3H4, and allene, respectively due to the unavailability of π electrons. Owing to the presence of orthogonal cumulative π bonds in allene, its complex procures higher deformation. The unsaturated C3H6 and C3H4 amplify the significant π-electron transfer to PPX resulting in higher deformation of PPX in the complex. In C4 isomers the DE of 1,2-butadiene is 0.6 kcal mol−1 higher than 1,3-butadiene and 2-butyne is 2.0 kcal mol−1 higher than 1-butyne. The Z-isomer of C4H8@PPX experiences 0.1 kcal mol−1 higher DE than the E-isomer. The DE of isobutylene@PPX and isobutane@PPX is 3.0 and 1.2 kcal mol−1 higher than C4H8@PPX and C4H10@PPX isomers, respectively. DE unveils the structural distortion of PPX is simulated by the unsaturated HC species. The DE of HCs in the HC@PPX complex is lesser than 1.0 kcal mol−1. Isobutane undergoes maximum distortion with DE of 0.8 kcal mol−1 followed by C4H8 (0.594 kcal mol−1) and C2H2 is least given in Table S3.† As the maximum deformation occurs in guest C4H8, in its complex PPX deforms least. The higher distortion of isobutane is reflected by its bulkier size and larger van der Waals sphere whereas the C2H2 has the least due to its horizontal orientation within the cavity of PPX.
Gas | ΔE | BSSE | ΔE + ZPE | DE | ΔG298 | ΔH298 |
---|---|---|---|---|---|---|
a All the energy units are in kcal mol−1. | ||||||
CH4 | −17.9 | −18.0 | −17.2 | 1.0 | −8.4 | −16.8 |
C2H2 | −22.2 | −19.1 | −19.6 | 1.1 | −13.3 | −19.8 |
C2H4 | −28.5 | −23.5 | −26.9 | 0.9 | −15.6 | −26.8 |
C2H6 | −32.9 | −27.4 | −31.4 | 0.9 | −19.6 | −31.3 |
Allene | −30.4 | −26.3 | −28.1 | 3.7 | −14.9 | −28.3 |
C3H4 | −31.2 | −27.6 | −29.5 | 3.3 | −17.1 | −29.4 |
C3H6 | −37.9 | −34.1 | −36.1 | 2.8 | −23.9 | −36.0 |
C3H8 | −43.7 | −38.9 | −41.8 | 1.6 | −28.6 | −41.8 |
1,2-Butadiene | −42.6 | −38.2 | −40.9 | 1.6 | −28.1 | −40.7 |
1,3-Butadiene | −43.8 | −40.8 | −41.6 | 1.0 | −28.9 | −41.5 |
1-Butyne | −42.2 | −36.7 | −40.6 | 0.9 | −27.8 | −40.4 |
2-Butyne | −38.4 | −34.0 | −36.8 | 2.9 | −23.8 | −36.7 |
C4H8 | −48.4 | −42.9 | −46.1 | 0.4 | −32.5 | −46.1 |
(E)-C4H8 | −48.2 | −42.9 | −45.9 | 3.1 | −32.6 | −46.0 |
(Z)-C4H8 | −48.7 | −40.3 | −46.6 | 3.2 | −33.4 | −46.5 |
Isobutylene | −47.2 | −41.8 | −44.3 | 3.4 | −29.7 | −44.6 |
C4H10 | −54.7 | −48.6 | −52.6 | 1.5 | −38.9 | −52.6 |
Isobutane | −49.8 | −46.9 | −46.6 | 2.7 | −30.2 | −45.7 |
The binding energy (ΔE), BSSE, and ZPE corrected ΔE (ΔE + ZPE) are computed to determine the interaction strength between HC and PPX as given in Table 1. The ΔE is the energy difference between the HC@PPX and the sum of individual energy components of HC and PPX as expressed in eqn (1).
ΔE = EHC@PPX − [EHC + EPPX] | (1) |
The ΔE value portrays the interaction of HCs within PPX as fairly high owing to its soft interaction exerted between the Au(I) and HCs unravelled by the Pearson HSAB principle. The CH4 molecule strongly binds to PPX with a ΔE value of −17.9 kcal mol−1. The CH4 uptake capacity of PPX is exceptionally higher than many of the reported systems such as MOF-5 (−4.51 kcal mol−1),25 NI-MOF-74 (−7.78 kcal mol−1),26 Cu24(m-BDC)24(−8.27 kcal mol−1),27 TpPa1 COF (−5.39 kcal mol−1),28 Graphdiyne (−5.18 kcal mol−1)29 etc. The ultrahigh binding energy of CH4 in PPX is due to the strong interaction between the 6s levels of Au(I) and 2a1 orbital of CH4 and the charge transfer through level hybridization from CH4 to Au(I) and not via Au(I) unoccupied orbitals.30
An upraised ΔE value of C2 HCs observed with a lowering in the percentage of s character unveils the stability of C2H6 over C2H2 and C2H4 with ΔE of −32.9 kcal mol−1 is 10.7 and 4.4 kcal mol−1 higher than C2H2 and C2H4 respectively. The higher binding energy of the C2H6 is due to the mixing of the occupied π−(CH3) orbitals with the Au(I) orbitals resulting in an electron charge donation from C2H6 to Au(I) and a partial de-occupation of the π−(CH3).31 The partial de-occupation is reflected in the lowest unoccupied molecular orbital (LUMO) that is formed by the interaction of the occupied π−(CH3) of ethane and unoccupied s-orbitals of Au(I). The C2H4 stabilizes the π mode adsorption and strengthens the Au–C bond by readily balancing between donation and back-donation of electrons which strongly influence the molecular adsorption.32 The C2H2 is oriented in such a fashion that it experiences interaction with the Au(I) px orbital parallel to it via dative interaction by the Dewar–Chatt–Duncanson scheme resulting donation of filled πy orbitals of the ligands into the vacant s orbitals of the metal atoms and the back donation.33 In a confined PPX region, the overlap of orbitals experiences repulsions, hence a steady decline in energy is observed. The binding energy of C2 HCs in PPX is higher than reported benchmark study such as PCN-245 (C2H4:C2H6: −5.5 kcal mol−1:−6.5 kcal mol−1),34 HKUST-1 (C2H4: −5.5 kcal mol−1),35 IRMOF-8 (C2H4:C2H6: −8.4 kcal mol−1:−8.4 kcal mol−1),36 CuBTC(C2H4: −7.1 kcal mol−1),37 Ni(BODC)(TED) (C2H4:C2H6: −6.7 kcal mol−1:−7.9 kcal mol−1),38 RPM3-Zn (C2H6: −12.4 kcal mol−1),39 Fe-MOF-74(C2H2:C2H4:C2H6: −12.3 kcal mol−1:−12.3 kcal mol−1:−8.8 kcal mol−1),40 MUF-15 (C2H4:C2H6: −8.4 kcal mol−1:−8.8 kcal mol−1),41 Fe-MOF-74(C2H2:C2H4:C2H6: −6.3 kcal mol−1:−5.6 kcal mol−1:−6.1 kcal mol−1),42 etc.
The ΔE value of C3 HCs is also remarkable where C3H8 possesses maximum binding energy of −43.7 kcal mol−1 which is 13.3, 12.5 and 5.8 kcal mol−1 higher than of allene, C3H4, and C3H6, respectively. The absence of π-electron density in C3H8 portrays higher ΔE value while the existence of orthogonal π-electron density in allene reduce the binding strength to −30.4 kcal mol−1. The ΔE difference in allene/C3H4 (0.8 kcal mol−1), allene/C3H6 (7.5 kcal mol−1), C3H4/C3H6 (6.7 kcal mol−1) is substantial enough for facile separation of individual C3 HCs from their respective mixtures. The obtain ΔE is substantially higher than that of CuBTC (C3H8:C3H6: −10.8 kcal mol−1:−11.1 kcal mol−1),43 Fe-MOF-74 (C3H8:C3H6: −9.9 kcal mol−1:−14.3 kcal mol−1),40 pyr@Cu-BTC (C3H8:C3H6: −8.7 kcal mol−1:−11.6 kcal mol−1),44 UiO-66 (C3H8: 9.9 kcal mol−1),45 Cu-DTO (C3H8:C3H6: −14.0 kcal mol−1:−33.5 kcal mol−1),46 tpt-Mg-MOF-74 (C3H8: 8.15 kcal mol−1),47 WOFOUR-1-Ni (C3H8:C3H6: −11.9 kcal mol−1:−11.3 kcal mol−1)48 etc.
Furthermore, C4 HCs and their ten isomers are investigated as their contributions to the natural gas mixtures cannot be ignored. C4H10 interacts strongly with an effective binding energy of −54.7 kcal mol−1 whereas the binding strength of 2-butyne is least (−38.4 kcal mol−1). The notable energy difference of 6.3 kcal mol−1 between C4H10 and C4H8 is ideal for the mainstream separation of their mixture. 1,3-butadiene (−43.8 kcal mol−1) exhibits higher binding strength than 1,2 butadiene (−42.6 kcal mol−1) and the binding energy of 1-butyne is −42.2 kcal mol−1 is 3.8 kcal mol−1 higher than 2-butyne. The energy difference between E-C4H8/Z-C4H8 and isobutane/Isobutylene is noteworthy. In general, the binding energy values of HC@PPX are higher than the ever-studied system. Hence, PPX can be a good material for selective adsorption and separation of light HCs. To reduce the errors associated with total binding energy, BSSE and ZPE corrected binding (ΔE + ZPE) energy values are taken into account and no such perceived changes are noticed.
To explore the stability of PPX + HC → HC@PPX inclusion complex, Gibbs free energy (ΔG) has been computed for HC@PPX. The ΔG value for all the studied complexes is found to be negative, proclaiming HC@PPX is exergonic. The higher ΔH reveals the stability of the system is enthalpy driven. The calculated ΔG and ΔH values are given in Table 1. The ΔG values of C1–C4 HCs spaned from −8.4 to −38.9 kcal mol−1 while the ΔH values are in the range of −16.8 to −52.6 kcal mol−1. The entropy change (ΔS) has a negligible effect on the complexation given in Table S3.† Both the ΔG and ΔH values are the least for CH4@PPX while the maximum for C4H10@PPX. Noteworthy ΔG values in C2 and C3 suggest PPX as a material of interest for C2/C3 HCs separation. A trivial ΔG value is observed in C4 HCs, Thus the separation of C4 HC needs an external potential.
To understand the effect of the inter-ligand and cooperativity within the cavity, for example, we have considered some of the π-bonded molecules. It is noticed that the PPX host accommodates only two HC within the cavity while the third one residing the open face of the PPX. Within the cavity C2H2, C2H4, allene, C3H4, 1-butyne, and 2-butyne experience π⋯π interaction at a distance of 3.189, 3.340, 3.364, 3.280, 3.506 and 3.405 Å respectively. The large size of C3H6, and other C4 HC dimers exerted by CH⋯π interactions. The CH⋯π distance in C3H6 is 2.781 Å whereas in 1,2-butadiene, 1,3-butadiene, 2-butyne, C4H8, (E)-C4H8, (Z)-C4H8 and isobutylene, the distance is 3.022, 2.799, 2.526, 2.743, 2.517, 2.743 and 2.513 Å respectively. The π interaction impetus the enhanced adsorption of HCs within the confined cavity.
The cooperative effect is an interesting factor in host–guest supramolecular interaction; therefore, we have computed the cooperative adsorption binding energy (CABE) and adsorption binding energy (ABE). The CABE is calculated by taking out the HC from the PPX host cavity, and single point calculations were carried out with the same orientations as obtained from the PPX cavity and the ABE is evaluated by relaxing the structure obtained from the HC@PPX complex. The ABE of HC@PPX is higher than the CABE value which suggests a cooperative interaction exists in the host–guest system. The CABE, ABE, and C–C bond distance of relax and constrained dimer is shown in Fig. S2.† The relax dimer bond distance is higher than the constrained dimer which claims the cooperative interaction engendered within PPX. C3H4, 1-butyne, 2-butyne, and (E)-C4H8 exhibit higher cooperativity interaction with an energy value of 5.1, 6.5, 6.0, and 8.4 kcal mol−1. Thus, it is well observed that the HCs adsorption within the PPX is asserted by the cooperative effect.
Natural gas sweetening is another stumbling block to the industrial separation of natural gas and acidic gas such as CO2, and H2S. Hence we have studied CO2 and H2S adsorption by PPX using B97D2/LanL2DZ level of theory. The ΔE of CO2 and H2S is −17.6 and −21.9 kcal mol−1 respectively while CH4, C2H6, C3H8 and C4H10 is −17.9, −32.9, −43.7 and −54.7 kcal mo, respectively. The ΔE of an alkane is higher than acidic gas, thus alkane is the superior adsorbate over CO2 and CS2. The ΔG298 value for CO2, H2S is −9.0 and −10.0 kcal mol−1 whereas CH4, C2H6, C3H8 and C4H10 are −8.4, −19.6, −28.6 and −38.9 kcal mol−1, respectively. Hence it will be difficult to separate CH4 in presence of an excessive amount of CO2 and CS2 as the complexation is exergonic and the complexation will be CO2@PPX and CS2@PPX instead of CH4@PPX.
CH4@PPX + CO2 → CO2@PPX + CH4 ΔG = −0.6 kcal mol−1 |
CH4@PPX + CS2 → CS2@PPX + CH4 ΔG = −1.6 kcal mol−1 |
The C2–C4 hydrocarbon separation can be easily separated in the presence of CO2 and CS2 as the substitution reaction undergoes endergonic.
C2H6@PPX + CO2/CS2 → C2H6@PPX + CO2/CS2, ΔG = 10.6/9.6 kcal mol−1 |
C3H8@PPX + CO2/CS2 → C3H8@PPX + CO2/CS2, ΔG = 19.6/18.6 kcal mol−1 |
C4H10@PPX + CO2/CS2 → C4H10@PPX + CO2/CS2, ΔG = 29.9/28.8 kcal mol−1 |
The enthalpy contribution in CO2, CS2, CH4, C2H6, C3H8 and C4H10 is −16.3, −21.3, −19.8, −28.3, −40.7, and −52.6 kcal mol−1, respectively imparts the reaction is enthalpy driven and is ideal to use PPX for the separation of HCs except for CH4 in presence of acidic gases.
The adsorption isotherm plot of 298 K is shown in Fig. 5(a) while the 77, 273, and 340 K are shown in Fig. S4(a), S5(a), and S6(a),† respectively. The adsorption is higher at 77 K due to the lowering of intermolecular repulsion between the CH4 molecule. 340 K results decrease in CH4 uptake, this can be attributed to high temperature favors the CH4 in the free gaseous state rather than adsorbed state, and also in high temperature the pore volume of the substrate decreases rendering lower adsorption.
Fig. 5 (a) Represents the adsorption isotherm of C1/C2 HCs (b) represents the adsorption isotherm of C3 HCs and (C) depicts C4 HCs at 298 K temperature and pressure 0.1–10 bar. |
The uptake of C2H6 (3.0–4.1 mmol g−1) is high compared to its isomers in the low-pressure region (298 K/0.1–1 bar). At high pressure, the C2H4 adsorption is prominent followed by C2H6 and C2H2. The gravimetric uptake of C2H2, C2H4, and C2H6 are 2.80, 4.26, and 4.22 mmol g−1 respectively at ambient conditions while at 10 bar the uptake increases to 4.37, 4.74, and 4.50 mmol g−1 respectively. At 77 K/1 bar the uptake varies C2H2 > C2H4 > C2H6 whereas at 273 K the adsorption strength is C2H4 > C2H6 > C2H2 up to 4 bar and thereafter C2H4 > C2H2 > C2H6. The 340 K is contrary to 273 K i.e. C2H6 > C2H4 > C2H2. This reveals at low-temperature the uptake of unsaturated HC is more effective whereas at high temperatures the C2H6 adsorption is significant. The adsorption isotherms plot of C2 HC at ambient condition is shown in Fig. 5(a) while Fig. S4(a), S5(a), and S6(a),† represents the isotherm plot of 77, 273, and 340 K respectively.
A type-IV adsorption isotherm has been observed for allene shows significant uptake of 3.76 mmol g−1 at 298 K/2 bar whereas other C3 isomers exhibit type-I at 298 K. At ambient conditions, the uptake of allene, C3H4, C3H6, and C3H8 is 3.51, 3.60, 3.57 and 3.12 mmol g−1, respectively are shown in Fig. 5(b), increasing pressure up to 1.5 bar results higher allene uptake followed by C3H4, C3H6, and C3H8. The gravimetric capacity of C3H8 is lower than other C3 isomers due to a larger kinetic diameter of 4.9 Å. At 77 K [Fig. S4(b)†], the adsorption isotherm of C3H4 is moderately higher than allene with increasing pressure and C3H6 adsorption is notably higher than C3H8. At STP [Fig. S5(b)†] the allene uptake is very high and high pressure portrays the C3H4 and allene have an equal propensity for adsorption. Thus, 273 K and 0.1–4 bar could be the ideal setup for selective adsorption of C3 HCs. The adsorption of allene (3.29 mmol g−1), C3H4 (3.20 mmol g−1), and C3H6 (3.33 mmol g−1) swing trivially at 340 K [Fig. S6(b)†], and the uptake of C3H8 (2.90 mmol g−1) is lowest. Increasing pressure reduces the C3H6 adsorption while inducing allene uptake. Unusual behavior is observed in C3H4 uptake capacity at a high-pressure and overall range of temperature which coincides with the allene adsorption isotherm. This can be explained by the activation of π-bond of allene and alkyne by Au(I) attributed to relativistic effects. Though Au(I) is coordinated with the electron-rich NHC groups which propel the π-activation more significantly with C3H4 rather than allene in high-pressure regions.
C4 HC uptake by PPX is in the ideal range, the highest adsorption is observed for 1,3-butadiene and isobutane is least at ambient condition and 1,2-butadiene follows type-IV adsorption while others are of type-I shown in Fig. 5(c). A preference analysis of C4H10, isobutane, 1,3-butadiene and isobutylene over other C4 HCs is investigated owing to its presence in the natural gas mixture. At ambient setup, the uptake capacity of C4H10, isobutane, 1,3-butadiene, and isobutylene are 2.32, 1.93, 2.95, and 2.36 mmol g−1, respectively. The remaining HC adsorption is in the ideal range. Except for 1,3-butadiene (3.08 mmol g−1), other HCs do not have notable uptake in cryogenic condition/1 bar as shown in Fig. S4(c).† An unrelieved uptake of HC has been observed at ambient and STP conditions [Fig. S5(c)†]. At 340 K/1 bar, the adsorption of 1,2-butadiene is high and with increasing pressure, the isobutylene shows steep growth in the adsorption as shown in Fig. S6(c).† The other C4 HC uptake strength is monotonous due to their equal kinetic diameter. Thus, the simulation study reveals that PPX might be a good candidate for the capture of C1–C4 HC under ambient conditions. Further, we have performed the density distribution of the center of mass of HC molecules in GCMC simulation as shown in Fig. S7.† The results from GCMC simulation at ambient condition indicate that the guest HC molecules lies at the center and also adhere to the outer open site of PPX due to the presence of the –N atom in the imidazolylidene motif. C1 and C2 HC show higher density coverage than C3 and C4 due to smaller size and lower kinetic diameter. The thermodynamic view of HC adsorption using PPX at 298 K has been determined using the isosteric heat of adsorption (Qst) and is shown in Fig. 6.
Fig. 6 (a) Represents the isosteric heat of adsorption (Qst) of C1/C2 HCs (b) depicts Qst of C3 HCs and (c) portrays Qst of C4 HCs at 298 K temperature. |
In C1/C2 HC, the Qst of CH4 and C2H2 are constant over the entire pressure range (0.1–10 bar) whereas C2H4 and C2H6 increase slowly from 0.1–1 bar and afterward remain constant. The increasing Qst value depicts the adsorbate interaction with the PPX and lateral interaction between the adsorbate molecules which is corroborated by the DFT and cooperative binding study. The constant Qst portrays the adsorption that occurs in the energetically homogenous adsorbent which is due to the effect of capillary condensation in the mesoporous material. The onset of capillary condensation results in constant Qst induced by confinement within the pore network. Akin to C2H4 and C2H6, C3 and C4 HCs exhibit the increase in Qst value from 0.1–1 bar enunciates the adsorbate–adsorbent interaction and constant from 1–10 bar attributed by capillary condensation.
Fig. 7 Ideal adsorption solution theory adopted selectivity of HCs by PPX (a) selectivity of C1/C2 HC (b) C3 HCs (C) C4 HCs at 298 K and 0.1–10 bar. |
The adsorption of C3 HCs is steady as compared to the other reported results. The C3H4 adsorption is less studied and only MFM-300(In)62 adsorption capacity is 6.30 mmol g−1 which is higher than PPX. The C3H6 adsorption in PPX is 3.57 mmol g−1 is higher than PCN-207-FA (2.17 mmol g−1),66 UMCM-151(3.07 mmol g−1),51 MOF-801(3.50 mmol g−1).67 The C3H8 adsorption in PPX is 3.12 mmol g−1, moreover, the reported MOFs show promising C3H8 adsorption compared to the PPX. Although the uptake of C3 HCs in PPX is quite low but the selectivity is well noticeable. The C4 HC uptake is more propitious and it is observed that the uptake capacity of (E)-C4H8, (Z)-C4H8, 1,3-butadiene, C4H8, isobutylene, C4H10, and isobutene using PPX is higher than the M-Gallate (M = Ni/Mg/Co),68 ZU-36-CO,69 TMOF-1/ZU-619,70 M2(dobdc)/M2(m-dobdc) (M = Mn/Fe/Co/Ni)71 MOFs shown in Table S6.† Only FJI-H19 (ref. 72) and NKM-101a73 show better uptake than PPX.
The energy decomposition analysis (EDA) has been carried out using dispersion corrected hybrid functional B3LYP-D3 with triple-zeta two polarization function TZ2P implemented in the ADF (2017) program.79–82 Zeroth order relativistic approximation (ZORA) method has been incorporated for computation as gold exhibits a large relativistic effect.83 In EDA, the fragment-based computations are invoked to treat PPX as one fragment and the guest molecules are other fragments. The total energy is decomposed into four major energy components viz. electrostatic (Eelct), dispersion (Edisp), orbital (Eorb), and the Pauli term (Epau). The former three energy terms are attractive while the latter one is repulsive. The Eelct term represents the classical electrostatic interaction that exists between the two charged species namely host and guest after complexation. Edisp term arises due to the electronic fluctuation between two adjacent molecules attributed to temporary dipoles in the complex. Eorb is the orbital interactions initiated by the overlapping of occupied and unoccupied atomic orbitals of the host–guest system. The Epau is the Paulis repulsive term exerted by the overlapping of an identical quantum number between the two interacting fragments.
The real-space visualization of interaction has been studied using the NCI index which identifies noncovalent interactions solely based on electron density and its derivatives.84 NCI index is based on the 2D plot of reduced density gradient (RDG), s as a function of electron density ρ(r):
(2) |
The plot where RDG value approaches zero is identified as the noncovalent region.84 Analysis of electron density determines the binding strength of interacting molecules. However, both attractive and repulsive interactions have a positive correlation with electron density in the corresponding region.85 A real space function sign (λ2) multiplied by electron density ρ(r) is indeed to define the attractive and repulsive term where the curvature λ2 is the second eigenvalue of the Hessian density matrix.86 A negative and positive value of λ2 is identified as strong attractive and repulsive interactions respectively while for van der Waals interaction the value of λ2 approaches to zero.86,87 Multiwfn program has been used to compute the NCI index whereas Chemcraft for the gradient isosurface visualization, where the color codes decide the value of sign (λ2)ρ(r).88,89 Blue and green color represents the attractive hydrogen bonding and van der Waals interaction whereas red color depicts the repulsive strain ensure by steric crowding.
Grand canonical Monte Carlo (GCMC) simulation has been carried out to determine the maximum HCs uptake capacity of PPX at an adsorption temperature of 77, 273, 298, and 340 K and pressure from 0.1 to 10 bar. For each simulation 1 × 105 GCMC cycles were used for equilibration and 1 × 105 cycles were used to assess the ensemble average properties. The simulations were performed using the sorption module of Material Studio.90,91 The universal force field (UFF) parameter has been used to describe the interatomic interaction with a cut-off distance of 12.5 Å has been used for all Lennard–Jones interactions.92 Long-range electrostatic interactions were computed using the Ewald summation method.93 The simulation has been carried out in a 1 × 1 × 2 supercell with imposed periodic boundary conditions in all directions.
The adsorbate uptake capacity using PPX crystal obtained from the GCMC simulation is further utilized for the computation of selective separations. Langmuir adsorption isotherm model has been used to examine the selectivity of individual gas components over multi-component using ideal adsorption solution theory (IAST).94
(3) |
Footnote |
† Electronic supplementary information (ESI) available: Includes physical parameters of pristine and encapsulated HCs, DE of HCs, charge transfer properties, Au–Au bond distances, change in entropy, energy decomposition values and their percentage contribution, cooperative effect bond distances, NCI-RDG plot, GCMC adsorption isotherm plot, selectivity graph, and Cartesian coordinates, etc. See DOI: https://doi.org/10.1039/d2ra07183e |
This journal is © The Royal Society of Chemistry 2022 |