Article https://doi.org/10.1038/s41467-023-41227-4 Regulating adsorption performance of zeolites by pre-activation in electric fields Kaifei Chen 1, Zhi Yu1, Seyed Hesam Mousavi 1, Ranjeet Singh1, Qinfen Gu 2, Randall Q. Snurr 3, Paul A. Webley4 & Gang Kevin Li 1 While multiple external stimuli (e.g., temperature, light, pressure) have been reported to regulate gas adsorption, limited studies have been conducted on controlling molecular admission in nanopores through the application of electric fields (E-field). Here we show gas adsorption capacity and selectivity in zeolitemolecular sieves canbe regulatedby an external E-field. ThroughE-field pre-activation during degassing, several zeolites exhibited enhanced CO2 adsorption and decreased CH4 and N2 adsorptions, improving the CO2/CH4 and CO2/N2 separation selectivity by at least 25%. The enhanced separation performance of the zeolites pre-activated by E-fieldwasmaintained inmultiple adsorption/desorption cycles. Powder X-ray diffraction analysis and ab initio computational studies revealed that the cation relocation and framework expansion induced by the E-field accounted for the changes in gas adsorption capacities. These findings demonstrate a regulation approach to sharpen the molecular sieving capability by E-fields and open new avenues for carbon capture and molecular separations. Microporous materials generally refer to porous materials with pore diameters <20Å. Active regulation of the pore accessibility in micro- porous materials is the eventual goal of advanced adsorption-based molecule processing. It has been reported that some porous materials experience framework deformation or vibration of pore-keeping groups/ions in response to temperature, light, guest accommoda- tion, and other stimuli, thus altering the pore accessibility and improving gas separation, storage, and recognition1. Zeolites, composed of aluminosilicate crystalline frameworks and extra-framework cations, are widely used microporous materials in adsorption, molecular sieving, and catalytic processes. Although they are often regarded as rigid frameworks with limited flexibility, some small-pore zeolites, such as low silica potassium chabazite, K-KFI, and RHO, show a special molecular trapdoor mechanism that exclusively discriminates guest molecules2–4. Their pore accessibility can be regulated by the extra-framework cations which act as the “door- keeping” ions. The trapdoor cation allows the admission of gas moleculeswith strongpolarizability or polarity, such asCO2, which can sufficiently interact with the door-keeping cation to lower the energy barrier and cause the cation to deviate temporarily from the pore aperture center. It rejects “weak” gas molecules, such as CH4 and N2, which lack interaction ability with the cation and are unable to lower the energy barrier for pore admission5–7. Furthermore, a temperature increase amplifies the thermal oscillation of the trapdoor cation, increasing the admission chance of the guest molecules8. As a result, above a threshold temperature, the adsorption capacity of “weak” gases can anomalously increase9. Although the temperature is an effective stimulus to control pore accessibility, heating and coolingof thematerial are relatively slowdue to heat transfer limitations. Inspired by the recent reports on electric- field-regulated structural transition10–12 and molecule transport in microporous materials13,14, we hypothesize that the application of an electric field (E-field) gradient would facilitate cation relocation and influence the access of molecules to the internal pore space of Received: 27 December 2022 Accepted: 25 August 2023 Check for updates 1Department of Chemical Engineering, The University of Melbourne, Parkville, VIC 3010, Australia. 2Australian Synchrotron, ANSTO, 800 Blackburn Rd, Clayton, VIC 3168, Australia. 3Department of Chemical & Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL 60208, USA. 4Department of Chemical and Biological Engineering, Monash University, Clayton, VIC 3800, Australia. e-mail: Paul.Webley@monash.edu; li.g@unimelb.edu.au Nature Communications | (2023) 14:5479 1 12 34 56 78 9 0 () :,; 12 34 56 78 9 0 () :,; http://orcid.org/0000-0001-7523-1928 http://orcid.org/0000-0001-7523-1928 http://orcid.org/0000-0001-7523-1928 http://orcid.org/0000-0001-7523-1928 http://orcid.org/0000-0001-7523-1928 http://orcid.org/0000-0003-1546-3780 http://orcid.org/0000-0003-1546-3780 http://orcid.org/0000-0003-1546-3780 http://orcid.org/0000-0003-1546-3780 http://orcid.org/0000-0003-1546-3780 http://orcid.org/0000-0001-9209-4208 http://orcid.org/0000-0001-9209-4208 http://orcid.org/0000-0001-9209-4208 http://orcid.org/0000-0001-9209-4208 http://orcid.org/0000-0001-9209-4208 http://orcid.org/0000-0003-2925-9246 http://orcid.org/0000-0003-2925-9246 http://orcid.org/0000-0003-2925-9246 http://orcid.org/0000-0003-2925-9246 http://orcid.org/0000-0003-2925-9246 http://orcid.org/0000-0003-1700-297X http://orcid.org/0000-0003-1700-297X http://orcid.org/0000-0003-1700-297X http://orcid.org/0000-0003-1700-297X http://orcid.org/0000-0003-1700-297X http://crossmark.crossref.org/dialog/?doi=10.1038/s41467-023-41227-4&domain=pdf http://crossmark.crossref.org/dialog/?doi=10.1038/s41467-023-41227-4&domain=pdf http://crossmark.crossref.org/dialog/?doi=10.1038/s41467-023-41227-4&domain=pdf http://crossmark.crossref.org/dialog/?doi=10.1038/s41467-023-41227-4&domain=pdf mailto:Paul.Webley@monash.edu mailto:li.g@unimelb.edu.au microporous materials. Ideally, the E-field required for regulating the adsorption would be below the breakdown voltage of the respective gases. To date, understanding of E-field stimulation on gas adsorption is very limited, and experimental evidence is also rare. Although from the 1970s, it has been known that the Coulombic effect induced by external E-fields can influence the interaction between gas molecules and the surfaceof non-porousmetal oxides15,16, howgasadsorption in a porous medium can be affected by E-field remained unknown. Recently, some studies emerged in the area of applying E-fields to organic frameworks. A study on a membrane made of zeolitic imida- zolate framework ZIF-8 indicated that a 500V/mmE-field could reduce gas transport due to the transformation of polymorphs13. Some molecular simulation studies suggested certain metal organic frame- works (MOFs) with flexible frameworks and/or polar linkers, such as MIL-53 (Cr), IRMOF-1 and IRMOF-7, could undergo deformation by intense E-fields despite the concern of exceeding the breakdown voltage10–12,17,18. Our latest experiments demonstrated that a less intensive E-field can decrease the CO2 adsorption capacity in a narrow poreMIL-53 (Al) due to the reduced charge transfer under the E-field14. However, to the best of our knowledge, none of the studies could produce any positive effect by using E-fields on gas sorption i.e., enhanced capacity or selectivity, which largely compromises the potential of this technique in important applications such as gas separation and molecular sieving. In this study, we demonstrate how adsorption performance can be enhanced under stimulation by E-fields. Instead of using flexible MOFs, the subjects we chose are zeolites with rigid inorganic frameworks. We found that applying an E-field in situ during adsorption does not influence the adsorption capacity of zeolites. However, applying the E-field during pre-activation can remarkably change the subsequent adsorption capacities of CO2, CH4, and N2, in various zeolites. Importantly, such enhancement of CO2 adsorption capacity and selectivity can be maintained in multiple cycles of adsorption/desorption which is critical to pressure swing adsorption processes. We show that the reloca- tion of the cations through an E-field polarization during pre- activation can induce lattice expansion and improve the subsequent molecular separation ability of zeolites. The complex underlying mechanism was investigated by gas adsorption mea- surements, surface area analysis, synchrotron powder X-ray dif- fraction (PXRD) analysis, and ab initio density functional theory (DFT) calculation. To the best of our knowledge, this study is the first experimental confirmation that an external E-field during pre-activation can enhance the gas separation performance, which has significant implications for carbon capture and natural gas separations. Results E-field pre-activation of zeolites We used two parallel electrodes to apply an external E-field perpen- dicular to a 0.5mm zeolite layer (Fig. 1a), and the assembly was placed into a glass cell. The electrodes were separated by the insu- lated zeolite plate so that there was no current in the system and the external power source only provided a static E-field. The thin slice of the zeolite and the pair of steel plates (i.e., electrodes) allowed for fast dissipation of heat. We first tried to apply in situ E-fields during the adsorption processes of CO2 and CH4 but could not find any evidence that E-fields can directly control the gas uptake of zeolites (Supplementary Fig. 1) within the range of E-field gradients applied. We hypothesized that the application of the E-field during sample activation may be more successful in influencing the adsorption behavior. Accordingly, we developed a different pathway to pre- activate the sample by E-field before adsorption in which, after the zeolite sample was completely degassed to remove gas residue, an E-field was immediately applied to the highly evacuated sample at a high temperature (Fig. 1b). By this approach, a more intense E-field, up to 800V/mm, can be applied to the zeolite, as the high vacuum condition can avoid gas breakdown caused by the large voltage gradient at low partial pressures. The latter was commonly observed for the adsorption processwhen the voltagewas higher than 200V. A high temperature is necessary for E-field pre-activation, as the in situ PXRD analysis indicated that the shift in XRD reflex under an E-field only occurred at high temperatures (>353 K) (see below & Supple- mentary Fig. 11), beyond the temperature range that is favorable for adsorption. Gas P Thermostat Power source - + Clamp Glass Porous stainless steel, d = 1 mm Zeolite, d = 0.5 mm E-field pre-activation Adsorption T (K ) T U Degassing Time 573 0 400 623 U (V ) Tadsorption a b 30 min E-field Fig. 1 | The E-field activationof zeolitemolecular sieves. aThediagramof the experimental cell for applying the E-field on zeolites.bTheprocess of E-field pre-activation and gas adsorption of potassium chabazite. Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 2 Gas adsorption in r2KCHA regulated by E-field pre-activation Adsorption experiments with CO2, CH4, and N2 were performed on potassium chabazite with a Si/Al ratio of 2.2 (r2KCHA, Si2.43Al1.03O7K1.32). Before adsorption, the r2KCHAwaspre-activatedby an alternating E-field (4000Hz) at the intensity of 800V/mm at 573 K under vacuum for 30min. All isotherms were measured at least three times with independent samples to prevent random errors from being confused with E-field effects. They had good reproducibility with random errors of <2% (as shown by the error bars in Fig. 2d). Gas uptakes by r2KCHA without E-field pre-activation were in good accordance with literature values3. After the E-field pre-activa- tion, the adsorbed amount of CO2 at 273 K increased from4.27mmol/g to 4.42mmol/g (Fig. 2a). In contrast, the accessibility of r2KCHA for CH4 and N2 was significantly reduced, with an ~20% decrease in the adsorption capacity at 100 kPa (Fig. 2b, c). The E-field effect on a specific gas is consistent atdifferent temperatureswith small errors. As shown in Fig. 2d the error bars were much smaller than the uptake difference caused by the E-field pre-activation, excluding the possibi- lity that the effect of the E-fieldwas raised fromanalysis error. Through fitting the isotherms with theoretical adsorption models (Supple- mentary Fig. 5), we found the maximum adsorption capacity of CO2 was increased by 0.19mmol/g by the E-field pre-activation, and those of CH4 and N2 decreased by 0.27 and 0.61mmol/g, respectively. Nonetheless, the heat of adsorption for all three gases was decreased by the E-field (Supplementary Figs. 6–8). The opposite trends of gas uptake change induced by the E-field for CO2 and CH4 and N2 provide an opportunity to improve the gas separation selectivity. The pure component selectivities19 for component A against component B (A/B) was calculated by Eq. (1): Selectivity = xA=yA xB=yB ð1Þ where x is the adsorbed-phase concentration and y is the correspond- ing gas-phase concentration. At 100 kPa and 273 K, the pure component selectivities of unac- tivated r2KCHA for CO2/CH4 andCO2/N2 pairs are 3.69 and 6.83, which were increased to 4.64 and 9.73, respectively, after the pre-activation of an E-field, demonstrating impressive improvements for 26% and 42%. (Fig. 2e). Notably, the selectivity enhancement at low partial pressures was even more significant (Fig. 2e). After the thorough degassing, the E-field pre-activated r2KCHA can completely recover its adsorption capacity of CO2, and partially recover the adsorption capacities of CH4 and N2 (Fig. 2f). It is unclear why only CO2 can completely recover its adsorption capacity. We speculate that the CO2 molecules have a stronger interaction with the framework and cations than theN2 andCH4molecules. Therefore, CO2 molecules canhelp to lower the energy barrier for the relaxation of any structural transition introduced by the E-field pre-activation. In addition, a direct E-field (as opposed to an alternating field) with the sameE-field intensity of 800V/mmwas applied to the r2KCHA layer to investigate if the current direction influences the effect of E- fields, as the electric charge only flows in one direction in the direct current while in the alternating current, it changes direction periodi- cally. The uptake changes of CO2 and CH4 led by a direct E-field shared the same tendency as those resulting from an alternating E-field (Supplementary Fig. 9). This is to be expected as the zeolite layer was composedof isotropic powders randomly locatedwithin the layer and, on average, equally exposed to E-fields in different directions. Framework expansion of r2KCHA induced by E-fields The results above raise thequestionof how the applicationof anE-field during sample pre-activation causes the change in the adsorption capacities of multiple gases. A similar E-field induced reduction in gas permeance has been reported before in the context of gas transport through a ZIF-8 membrane13. It was attributed to the phase transfor- mation of the ZIF-8 lattice and the inhibition of linker rotations. The rotatable linkers inMOFs can respond to E-fields by changing the bond angles and linker orientations17. However, these explanations do not apply to zeolites, whose frameworks lack moveable linkages or func- tional groups20. Zeolite frameworks showed limited relaxation under pressure-induced hydration or variable temperature, and the changes in the unit cell were always accompanied by the cation relocation21,22. Therefore,we speculated that anE-fieldmay impose some influenceon the structure of zeolites as that observed in the above studies. To verify this conjecture, the r2KCHA pre-activated by the E-field was analyzed by PXRD before performing adsorption measurements to avoid possible relaxation. As a comparison, the normal r2KCHA without E-field pre-activationwas analyzed under the same conditions. The PXRD patterns of E-field pre-activated r2KCHA clearly showed a uniform shift of Bragg peaks to the smaller angle (Fig. 3a), suggesting an expansion of the framework23. This was verified by the cell volumes calculated using the cell refinement function in JADE (Fig. 3a). Syn- chrotron PXRD measurements with the in situ E-field application confirmed this phenomenon and further indicated that the E-field imposed a sustainable effect on r2KCHA, as removing the E-field did not immediately bring the Bragg peaks back to the original posi- tions (Fig. 3b). Cation relocation in r2KCHA induced by E-fields Weseek to explain themechanismof the E-field leading to a framework expansion and its association with the gas uptake changes. It has been widely reported that the migration of extra-framework cations can induce a relaxation of the aluminosilicate framework21. To date, researchers have only recognized the cation relocation in zeolites triggered by thermal oscillation9 or guest-host interaction3,7. No attempt has beenmade to apply E-fields on the zeolite to induce cation migration. We hypothesize that the framework expansion observed here is associated with a cation relocation caused by the E-field. To validate this hypothesis, we optimized the lattice structure of r2KCHA with the K+ ions located at different sites using ab initio DFT. Chabazite is a kind of small-pore zeolite consisting of double-six ring (D6R) prisms linked by the tilted four-membered ring (4MR)24. The three-dimensional structure constitutes a large supercavity accessed by six eight-membered rings (8MR). Each unit cell included three D6Rs and one and a half supercavities7. There are four general cation posi- tions in dehydrated chabazite25: site SI at the center of the D6R prism, site SII at the triad axis of the D6R prism but displaced towards the supercavity, site SIII in the supercavity above the 4MR, and site SIII’ in the 8MR (Fig. 4a). Therefore, there are three SI sites, six SII sites, nine SIII sites, and nine SIII’ sites in each unit cell. The cation number is dependent on the Si/Al ratio of zeolite, and in r2KCHA, there is a total of nine extra-framework K+ in each unit cell. Sincemonovalent cations, such as K+, energetically prefer site SIII’25,26, the stable initial structure of r2KCHA unit cell is composed of the framework and nine K+ that are all located at SIII’. By changing the unit cell volume, we plotted the relative energy profile with one K+ located at different sites through DFT calculation. When the system energy reaches the minimum, the structure of r2KCHA is the most stable. As shown in Fig. 4b, when all cations are located at SIII’, the unit cell volume of r2KCHA is 2152 Å3. A similar approach has been employed to calculate the unit cell volume of r2KCHA with one of the K+ relocated from SIII’ to SI, SII, and SIII, respectively (Supplementary Fig. 12). The unit cell expands if K+ relo- cates from the original SIII’ site to another site, which is well correlated with our observations in PXRD analysis. For instance, the unit cell volume increased to 2217 Å3 when one of the K+ ions relocated to the SI site (Fig. 4c, d). The simulated PXRDpatterns of r2KCHA also indicated significant left shifts that coincide with the experimental patterns (Fig. 4e). Meanwhile, the experimental PXRDpatterns indicated a long- Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 3 lasting peak shift (Fig. 3b), implying the cation can stay at the new site after the relocation. By employing the Berry phase method, the spontaneous polar- ization of r2KCHA inducedby the cation relocationwas calculated. The polarization difference resulting from the relocation of one K+ from SIII’ to SI was 1.2, 1.1, and 2.9 C/m2 along the x, y, and z axes, respectively. It was 10 times higher than the spontaneous polarization of sodalite at room temperature27, suggesting that an E-field pre-acti- vation can lead to a remnant polarization in r2KCHA. Experimental evidence for the cation relocation in r2KCHA is provided by the increase of BET surface area after the E-field pre- activation. As a typical trapdoor zeolite, the N2 surface area at 77 K of 0 20 40 60 80 100 120 0 1 2 3 4 5 Normal E-field pre-activated )g/lo m m( ekatp U P(kPa) CO2 0 20 40 60 80 100 120 0.0 0.2 0.4 0.6 0.8 1.0 1.2 Normal E-field pre-activated CH4 )g/lo m m( ekatp U P (kPa) 0 1 2 3 4 5 6 7 0.00 0.05 0.10 0.15 0.20 0.25 0.40 0.45 0.50 0.55 0.60 Normal E-field pre-activated Relaxed N2 20 40 60 80 100 0 5 10 15 20 25 30 35 CO2/CH4 Normal CO 2 /N 2 Normal CO2/CH4 E-field pre-activated ytivitceleS P (kPa) CO 2 /N 2 E-field pre-activated 0.8 0.9 1.0 1.1 1.2 1.3 CH4 ①③ ② ② ③ ① ① ③ ② 3.9 4.0 4.1 4.2 4.3 )g/lo m m( e katp U CO2 0.0 0.5 1.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 20 40 60 80 100 120 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 N2 )g/lo m m( ekatp U P (kPa) Normal E-field pre-activated 250 260 270 280 290 300 0.0 0.5 1.0 1.5 3.5 4.0 4.5 5.0 )g/lo m m( ekatp U Temperature (K) CO2 CH4 N2 E-field pre-activated Normal Normal E-field pre-activated Normal E-field pre-activated a b c d e f E-field Uptake Increase Uptake Decrease Uptake Decrease Error bar Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 4 r2KCHA is very low, as the K+ occupying the SIII’ site completely blocks the passage through the 8MR at low temperatures and hinders the diffusion of N2 3. Therefore, the N2 surface area of r2KCHA at 77 K is sensitive to the relocation of trapdoor cations. After being pre- activated by the E-field, the N2 adsorption of r2KCHA at 77 K sig- nificantly increased (Supplementary Fig. 13) and the BET surface area doubled from 11.73 to 21.67m2/g, confirming that some of the door- keeping cations have moved away to admit the originally rejected N2. According to the maximum adsorption amount of N2 in r2KCHA regardless of the trapdoor effect (Supplementary Table 3), the relo- catedK+ induced by the E-field pre-activationwas estimated to be ~12%. In conclusion, it is plausible that the E-field induced framework expansion of r2KCHA was consistent with K+ relocation from the ori- ginal SIII’ site. Gas admission regulated by the E-field induced cation relocation To explain how the cation relocation and framework expansion influ- ence gas adsorption, we placed the gas molecules into the supercavity of r2KCHA and calculated the adsorption energies with one K+ located at SIII’ and SI, respectively (other K+ were always at SIII’), as the K+ relocation to SI could induce a moderate expansion of the framework (Fig. 4c). The adsorption energy Eads is defined as the difference between the total energy of the r2KCHA-gas complex and the total energy of r2KCHA augmented by the total energy of an isolated gas molecule (Eq. (2)), Eads = Etotðr2KCHA+gasÞ � Eðr2KCHAÞ � EðgasÞ ð2Þ which represents the energy released during adsorption28. If we assume that one of the K+ ions in the r2KCHA unit cell relocates from the original SIII’ site to the SI site after the E-field pre-activation, then the change of adsorption energy Δ |Eads| caused by the E-field can be calculated by Eq. (3). ΔjEadsj= jEads SIðE�field pre�activatedÞj � jEads SIII’ðNormalÞj ð3Þ The Δ |Eads| of CH4 was −0.03 eV, suggesting the energy released during CH4 adsorption will be lower if the K+ relocates from SIII’ to SI, which is consistent with the decreased heat of adsorption obtained from experimental isotherms (Supplementary Figs. 6, 8 and Supple- mentary Table 2). The decrease in adsorption energy will reduce the gas adsorption capacity, which explains the situation observed in CH4 and N2 adsorption. As shown in Fig. 5a, b, the charge rearrangement (the enrichment or depletion of the electron density) after CH4 adsorption appears to be weakened when the K+ relocates from SIII’ (Fig. 5a) to SI (Fig. 5b). It suggested that compared with the K+ at SIII’, the CH4molecule had amore-limited interaction with the K+ located at SI, which contributed to the decreased adsorption energy after the K+ relocation. Meanwhile, the framework expansion induced by the K+ relocation (Fig. 3a) will also reduce the affinity between the host fra- mework and the guest molecules, as an expanded framework will decrease the electron density around the adsorption sites29. As “weak” gases with little interactionwith cations, CH4 andN2 are excluded by pore-blocking K+ when the temperatures are lower than the critical admission temperature (Tc) 8. Theoretically, the relocation of trapdoor K+ can lead to the successful passage of gas molecules through the 8MR, facilitating the adsorption of CH4 and N2. However, in our cases, such a gate-opening effect caused by the E-field induced cation relocation was not observed, as the adsorption capacities of CH4 and N2 after E-field pre-activation were decreased. It can be Fig. 2 | Effect of E-field pre-activation on gas adsorption in r2KCHA. Adsorption isotherms of (a) CO2, (b) CH4, and (c) N2 at 273 K without and with the E-field pre- activation. d Isobars of CO2, CH4, and N2 at 100 kPa. The error bars are marked on the graph. Dash lines are for the guidance of the eye only. Complete isotherms are shown in Supplementary Figs 2–4. e Selectivities of CO2/CH4 and CO2/N2 were increased by the E-field pre-activation. f Relaxation by degassing at 623 K of the E-field pre-activated r2KCHA for CO2 adsorption at 294K and 101 kPa, CH4 adsorption at 273K and 100 kPa, and N2 adsorption at 294K and 100 kPa. ① is the first adsorption measurement without E-field pre-activation. After ①, the sample was degassed, pre-activated with the E-field, and measured for the uptake, which was labeled as ②. After ②, the sample was degassed again, and the adsorption measurementwas conductedwithout E-fieldpre-activation, whichwas labeled as③. The E-field intensity is 800V/mm. Lines are guides for the eye. Source data are provided as a Source Data file. 5 10 15 20 25 30 35 40 45 50 ).u.a( ytisnetnI 2θ (degree) Normal E-field pre-activated Vol = 2470.85 Å3 Vol = 2499.87 Å3 a b 4.70 4.72 4.74 4.76 4.78 4.80 5th 0 V/mm 4th 267 V/mm 3rd 200 V/mm 2nd 100 V/mm 1st 0 V/mm ).u.a( ytisnetn I 2θ (degree) 353 K E-fieldLeft shift Fig. 3 | PXRDanalysis of r2KCHA. a PXRDpatterns of r2KCHAbefore and after the E-field pre-activation. The cell volumes were calculated by the cell refinement method of MDI Jade 6.0. The fitted profiles are shown in Supplementary Fig. 10. b PXRD patterns of r2KCHA with in situ E-field application at 353 K. The r2KCHA sample was successively exposed to E-fields of 0, 100, 200, 267, and 0V/mm. The E-field induced the left shift of the pattern. PXRD patterns at other temperatures are in Supplementary Fig. 11. Source data are provided as a Source Data file. Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 5 explained for two reasons. Firstly, the relocation of partial K+ cannot fully allow gas molecules to diffuse into the internal cavities, because other K+ that stay at SIII’will block the diffusion pathway. Although the BET surface area of r2KCHA increased after the E-field pre-activation, it was still far less than normalmicroporousmaterials, implying not all of the trapdoor K+ were relocated. Secondly, the adsorption tempera- tures (252, 273, and 294K) were beyond the Tc of CH4 (250–279 K) and N2 (243–266 K)8, atwhich the pore-blocking K+ has already temporarily deviated by thermal oscillation9. Therefore, the influence of the expanded framework and the decreased adsorption energy will dom- inate to decrease the adsorption capacity. Once the adsorption tem- perature is low enough, the gate-opening effect led by the cation relocation will be dominant to increase the adsorption capacity, as in the case of N2 adsorption in r2KCHA at 77 K (the BET surface area analysis). As for CO2, the |Eads| fromDFTwas decreased by0.013 eV after the E-field pre-activation, which was in agreement with the decrease in the heat of adsorption calculated from experimental isotherms (Supple- mentary Figs. 5, 8 and Supplementary Table 1). However, it seemed contradictory to the increased adsorption capacity of CO2 after the E-field pre-activation. Given this difference, an increase in adsorption sites due to the relocation of cationswas considered to account for the improvement of CO2 adsorption capacity. It has been reported that CO2moleculeswill prefer to occupy the center of 8MR in zeoliteswhile N2 cannot occupy that site 30. When an E-field was applied, some of the K+ located at SIII’ moved away, leaving a space for the CO2 molecule. We studied the system energies when the CO2 molecule passes through the unblocked 8MR and diffuses into the supercavity. When the CO2 occupies the 8MR, the system energy is the lowest (Fig. 5c), demonstrating that the SIII’ originally occupied by the K+ is a new adsorption site of CO2. In comparison, theCH4molecule cannot stay at the center of the 8MR, and its adsorption site remains inside the supercavity even if the SIII’ site is available (Fig. 5c). The reason is that the diameter of the 8MR in r2KCHA was ~3.8 Å, which is too small for CH4 (kinetic diameter = 3.8 Å) to locate, while CO2 (kinetic diameter = 3.3 Å) can more freely stay at the window site. Meanwhile, CO2 mole- cules have a stronger quadruple interaction with the framework than CH4molecules. Both factorsmade the center of the 8MR only ideal for CO2 30,31. Demonstration for other zeolites In addition to r2KCHA, some other zeolites were investigated, including trapdoor zeolites ZSM-25-Na (Na330[H30.7Al360.7Si1079.3O2880] 32, ZSM-25- K (K327.9Na2.1[H30.7Al360.7Si1079.3O2880]) 33, and the ionic liquid zeolite TMA-Y34,35 without trapdoor effect. The PXRD patterns of ZSM-25-K and TMA-Y under an E-field showed analogous shifts to the smaller angle as that of r2KCHA (Supplementary Figs. 14 and 15), suggesting that the E-field pre-activation can impose similar structural variations on differ- ent zeolites. A similar change in adsorption capacity induced by an E-field was also observed. Overall, the CO2 adsorption capacity of zeo- lites can be enhanced by E-fields while the adsorption of CH4 is com- monly suppressed as shown in Fig. 6a, b. This general pattern is suggestive of an attainable approach to improvingmolecular separation selectivity by E-fields. The exception was the CH4 uptake in ZSM-25-K at 252 K and 273 K, which significantly increased after E-field pre-activation (Fig. 6c and Supplementary Fig. 18). At 252 K and 273 K, which are below the Tc(CH4) of ZSM-25-K (275–350K), the K+ completely SI SII SIII SIII' 2120 2170 2220 2270 U ni t C el l V ol um e (Å 3 ) Cation Site SIII' SI SIII SII 1965 2006 2066 2127 2152 2190 2253 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1900 2000 2100 2200 2300 R el at iv e En er gy (e V) Unit Cell Volume (Å3) Relaxed unit cell SIII' a b c d e SIII' SIII SII SI 5 10 15 20 25 30 In te ns ity (a .u .) 2θ (degree) SIII' SIII SI Fig. 4 | Framework expansion of r2KCHA induced by cation relocation. a The schematic illustration of chabazite with four cation sites. Yellow: Potassium. Blue: Silicon or Aluminum. Red: Oxygen. b Changes of the relative energy over unit cell volume with all K+ located at SIII’. c The unit cell volumes of r2KCHA when one of the K+ locates at different cation sites. d The pathways for the K+ relocating from SIII’ to other cation sites. Yellow: Potassium. Blue: Silicon or Aluminum. Red: Oxy- gen. eThe simulatedPXRDpatterns of r2KCHAunit cellwith eight K+ ions locatedat SIII’ site and one K+ ion located at SI, SIII and SIII’ site, respectively. Source data are provided as a Source Data file. Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 6 blocked the 8MRof ZSM-25, leading to the lowadsorption capacity of CH4 33. Therefore, the relocation of trapdoor K+ induced by the E-field pre-activation can substantially reduce the energy barrier of CH4 molecules to pass through the 8MR, which dominated over the decreased adsorption energy caused by the framework expansion and finally increased theCH4 uptake. This phenomenonwas opposite to that observed in r2KCHA, in which the CH4 adsorption was reduced by the E-field pre-activation at the same temperatures. It is due to the particularity of trapdoor zeolites as the temperature is a crucial factor for oscillation and relocation of trapdoor cations9. The original state of trapdoor cations before E-field pre-activation will decide the priority of two opposite consequences brought by the cation relocation: the gate-opening effect that facilitates adsorption, and the decreased adsorption energy that impairs adsorption. For r2KCHA, the Tc(CH4) was 250–299 K8, implying that the trapdoor K+ was partially open above 250K, which diminished the gate-opening effect and emphasized the influence of the expanded framework after the E-field pre-activation. Based on the experimental phenomena and computational simulation,webelieve that an E-field can facilitate the cation relocation in zeolites. As demonstrated in previous studies, in the presence of an external stimulus, such as a high temperature, the cations will be able to overcome the energy barrier to relocate to a local minimum9. The E-field serves as the stimulus to energize the cation relocation in a way similar to a temperature-induced stimulus. When the cation relocates from its pristine location under anE-field, it canhold at its new location after the E-field was removed, because there will be an energy barrier for the cation to move back, as shown by the calculated energy profile in Supplementary Fig. 21. The E-field effect was reversible after the high-temperature degassing (Fig. 2f). However, degassing the zeolite at a low temperaturemayprolong the effectof E-field activation as themoderate temperature cannot facilitate the relocated cations to overcome the energy barrier to relax to its pristine location. This assumption was validated on TMA-Y. By degassing the E-field pre- activated TMA-Y at 393K, the E-field effect on TMA-Y can be retained for at least 5 cycles of adsorption-desorption, creating a high and long- lasting selectivity of CO2/CH4 (Fig. 6d). The decreased CO2 uptake in the third batch of adsorption resulted from the lower degassing tem- perature, which was the case in the TMA-Y without E-field pre-activa- tion as well. The lasting effect makes the E-field pre-activation applicable in the industrial adsorption process, such as pressure swing adsorption. The effect of E-field pre-activation on gas separation Column breakthrough experiments of TMA-Y for the mixture of CO2/CH4 (50:50, v/v) at 298 K at atmospheric pressure were con- ducted to evaluate the effect of E-field pre-activation on gas separation. During the breakthrough process, CH4 eluted at the start of the experiment while CO2 was detected later at the outlet gas. The E-field pre-activation of TMA-Y impaired the CH4 0 eV -0.047 eV ΔE R el at iv e en er gy le ve l 0 eV -0.042 eV -0.002 eV ΔE ba c Fig. 5 | Gas adsorption in r2KCHA simulated by the ab initio DFT calculation. a, b Electron density difference of the CH4-CHA complex with (a) all K+ located in SIII’ or (b) one of the K+ relocated from SIII’ to SI. The yellow iso-surfaces represent the electron accumulation region, and the blue iso-surfaces represent the electron depletion region. The iso-surfaces are plotted at 0.0001 e−/bohr3. c The relative energy level of r2KCHAwith oneCO2 (left) or CH4 (right)molecule passing through the unblocked 8MR. In r2KCHA, one of the K+ relocates from SIII’ to SI, leaving a space at the center of the 8MR. Purple: Potassium. Blue: Silicon. Silver: Aluminum. Red: Oxygen. Brown: Carbon. White: Hydrogen. Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 7 adsorption and accelerated its elution process (Fig. 7), which was consistent with the single-component adsorption isotherm. The breakthrough of CO2 did not show remarkable change after the E-field pre-activation of the adsorbent, which resulted from the competitive adsorption between CH4 and CO2, where the adsorp- tion capacity of theweak component CH4 is further depressed in the presence of the strong component CO2. The pre-activation by E-field amplified the suppression of CH4, leading to a sixfold increase of CO2/CH4 selectivity (calculated by themass ratio of CO2/ CH4 in the adsorbed gas mixture when the breakthrough CO2 con- centration was 5%). The CO2 adsorption capacity of TMA-Y calcu- lated based on the area integral under the breakthrough curve increased from 90.77 to 99.56mL/g after the E-field pre-activation. These results correspondedwith our prediction based on the single- component adsorption performance of TMA-Y and confirmed that the E-field pre-activation approach is effective in improving gas separation selectivity. Discussion We report an E-field pre-activation approach to regulating the gas adsorption capacity of zeolite molecular sieves. Such pre-activation involves applying an external E-field to zeolites during degassing, leading to the relocation of extra-framework cations and the expansion of the zeolite framework. After the E-field pre-activation, the vacant cation sites can be new adsorption sites for CO2 mole- cules, increasing the CO2 adsorption capacity. The framework expansion and the decrease of cation-guest interaction after the cation relocation resulted in the reduced adsorption energy of gas molecules, decreasing the adsorption capacities of CH4 and N2. Consequently, the adsorption selectivity of CO2 towards CH4 and N2 can be significantly improved by E-field pre-activation in various zeolites, including chabazite, zeolite ZSM-25, and zeolite TMA-Y. Meanwhile, the effect of E-field pre-activation can bemaintained for several cycles of adsorption-desorption by degassing the adsorbent at amoderate temperature. After a one-off E-field pre-activation, the a b 0 20 40 60 80 100 0.00 0.01 0.02 0.03 0.04 )g/lo m m( ekatp U P (kPa) Normal E-field pre-activated CH4 in ZSM-25-K 1 2 3 4 5 6 0.50 0.55 0.60 0.65 0.70 0.75 3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 393 K degas CH4 CO2 E-field pre-activated 473 K degas )g/lo m m( ekatp U Number of batch E-field pre-activated 2 3 4 5 6 6.5 7.0 7.5 8.0 E-field pre-activated C O 2/C H 4 ytivit celes Normal 2.0 2.2 2.4 3.8 4.0 4.2 TMA-YZSM-25-K )g/lo m m( ekatp U Normal E-field pre-activated ZSM-25-Na CO2 0.00 0.02 0.04 0.06 0.40 0.45 0.50 0.55 CH4 Normal E-field pre-activated )g/lo m m( ekatp U TMA-YZSM-25-KZSM-25-Na c d Fig. 6 | Effect of E-field pre-activation on gas adsorption in zeolites ZSM-25 and TMA-Y. Adsorption capacity of (a) CO2 and (b) CH4 in different zeolites (ZSM-25- Na, ZSM-25-K, and TMA-Y) at 101 kPa and 294K without and with the E-field pre- activation. The complete adsorption isotherms are in Supplementary Figs. 16–20. cAdsorption isothermof CH4 in ZSM-25-K at 252 Kwithout andwith the E-field pre- activation. d Repeated adsorption batches of CO2 and CH4 adsorptions in TMA-Y. The E-field pre-activation was only conducted before the second batch of adsorp- tion. The temperature represents the degassing temperature between two batches of adsorption trials. Source data are provided as a Source Data file. Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 8 CO2/CH4 selectivity of TMA-Y remained higher than that without E-field pre-activation in five subsequent adsorption tests. The gas separation selectivity of TMA-Y zeolite for CO2/CH4 in the binary breakthrough experiments showed a sixfold improvement by the E-field pre-activation. Notably, the E-field induced cation relocation can also provide a “gate-opening” effect to trapdoor zeolites, facilitating the admission of guest molecules that were originally blocked by the trapdoor ions. These findings reveal an E-field induced structural transition of zeolite that has never been dis- covered before. They also demonstrate how to use E-field to adjust the pore accessibility of microporous materials, which provides an approach to improving molecular separation and carbon capture. Methods Zeolite synthesis Three zeolites were prepared for this study: Chabazite, ZSM-25, and TMA-Y. Chabazite with the Si/Al of 2.2 was synthesized from HY zeolite (CBV400) following the reported procedure3. It was ion-exchanged by KCl solution twice to obtain fully exchanged potassium chaba- zite (r2KCHA). Specifically, 5 g of chabazite was refluxed in 200mL KCl (1 M) at 343 K for 24 h and then filtrated with deionized water 3 times. ZSM-25-Na zeolite was synthesized following the reported method36. ZSM-25-K zeolite was obtained by the ion exchange treat- ment of ZSM-25-Na. To be specific, 1 g of ZSM-25-Na was dispersed in 2mLof 1MKCl solution at roomtemperature for 2 h, and then filtrated with deionized water 4 times and dried at 353K overnight. The ion exchange procedure was repeated 6 times to get a fully exchanged ZSM-25-K. TMA-Y zeolite was prepared from NaY zeolite following the reported method, with ~30% of Na+ exchanged by the organic cation [C4H12N] +34. Sample degassing Before the E-field pre-activation and gas adsorption measurements, zeolite powderswere tableted into thin plates (0.5mm in thickness) by applying a pressure of 7000psi for 3min, cut into the size of 15mm (L) × 5mm (W), and heated stepwise to be thoroughly degassed under vacuum (<1.2 Pa) overnight. Specifically, the heating rate was 10K/min, and the sample was kept for 60min at 323 K, 373 K and 423 K, respectively before reaching the degassing temperature. The degas- sing temperatures for r2KCHA, ZSM-25, and TMA-Y were 623 K, 573 K, and 473 K, respectively, which were determined by their structural stabilities at high temperatures. E-field pre-activation of zeolites The E-field pre-activation of the zeolite samples was conducted during degassing. The zeolite layer was sandwiched between two porous stainless-steel plates which acted as electrodes and clipped by a fold- back clip which was separated from the electrodes by two pieces of glass. They were placed into a two-fold glass cell that can directly connect to the gas-dosing manifold of the adsorption measurement equipment. The glass cell was inserted with two metal hooks to hang the nichromewireswelded to the electrodes. The joint part of the glass cell was sealed with vacuum grease. The specific appearance of the sample and electrodes is shown in Fig. 1a. An external power source (Asterion Programmable AC/DC Power Source, California Instruments, USA) was connected to the metal hooks on the glass cell by wires to apply an E-field of 800V/mm (400V between 0.5mm) to the zeolite layer under high vacuum (0.01 Pa). For the E-field pre-activation, two batches of r2KCHA were exposed to two kinds of E-fields (a 4000Hz alternating current E-field and a direct current E-field) at 573K to investigate the influence of the current type. The TMA-Y was pre- activated by an alternating current E-field of 4000Hz at 473K. The ZSM-25 was pre-activated by a direct current E-field at 573 K. The procedure of the E-field pre-activation is shown in Fig. 1b. After 30min E-field pre-activation, the sample was rapidly cooled down under a vacuum by removing the heating jacket. The E-field was removed after the sample was cooled to room temperature. Gas adsorption measurements The adsorption isotherms of CO2, CH4 and N2 were recorded by a Micromeritics ASAP 2010 system. The equilibrium interval for CO2 was 60 s, and that for CH4 and N2 was 30 sec. During gas adsorption ana- lysis, the sample temperature was controlled by a Dewar filled with water for 294K, ice water for 273 K, and sodium chloride-ice mixture for 252 K. By supplementing the melted ice on time, the fluctuation of the temperature was strictly controlled within 1 K. Each set of adsorption measurements was conducted three times to quantify the error. Firstly, the normal adsorption isothermof the zeolite sample without an E-field pre-activation wasmeasured. Then the same sample was degassed again after adsorption and then pre-activated by the E-field before the second batch of adsorption measurements, which provided the adsorption isotherm with the E-field pre-activation. As a comparison, an in situ E-field of 400V/mmwas applied during the gas adsorptions of r2KCHA. Breakthrough gas separation The dynamic gas breakthrough experiment of TMA-Y for the dry mixture of CO2/CH4 (50:50, v/v) was carried out in a stainless-steel column (3.3mm inner diameter × 15.8mm) at 298K under 1 atm. The gas flow rates were all controlled by mass flow controllers (D07-7B, Beijing Sevenstar Flow Co., Ltd.). TMA-Y samples with or without E-field pre-activationwerepreviouslydegassed in the specialized tubes and then transferred into the breakthrough stainless steel column in the glove box. For a typical breakthrough experiment, 0.32 g of TMA-Y sample with a height of ~5.3mm were packed into the column and firstly purged with He flow for 1 h at 298K. Then the gas mixture of CO2/CH4 with the flow of 5mL/min was introduced into the column. The CO2 and CH4 concentrations in the outlet gas were detected using a gas analyzer and mass spectrometer (BELMass, BEL JAPAN. INC.) within 60min. 0 10 20 30 40 50 60 0.0 0.5 1.0 1.5 2.0 C t/C 0 Time (min/g) Without E-field pre-activation CH4 Without E-field pre-activation CO2 With E-field pre-activation CH4 With E-field pre-activation CO2 Fig. 7 | Experimental binary breakthrough curves for a dry gasmixture of CO2/ CH4 (50:50, v/v) on TMA-Y zeolites at 298K and 101.3 kPa. Ct and C0 denote the outlet and inlet concentrations, respectively. The red line denotes the TMA-Y sample has been pre-activated by an E-field of 800V/mm and 4000Hz, while the black line means the TMA-Y sample did not receive any E-field pre-activation. Source data are provided as a Source Data file. Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 9 Surface area analysis The specific surface area of r2KCHA without/with the E-field pre-acti- vation was determined by the Micromeritics 3Flex Adsorption Analy- zer with N2 at 77 K. After being pre-activated by the 800V/mm E-field, the r2KCHA was transferred into the sample tube of Micromeritics 3Flex in the N2 atmosphere for surface area analysis. The normal r2KCHA without E-field pre-activation was analyzed together as a comparison. The surface area was calculated by a standard Brunauer–Emmett–Teller (BET)method37 within the pressure range of 0.05 < P/P0 < 0.3. Powder X-ray diffraction analysis Powder X-ray diffraction (PXRD, Bruker D2 phaser) using Cu Kα radiation was used for the identification of possible framework trans- formations of r2KCHA after the E-field pre-activation. It was conducted in the air at room temperature. The lattice constants of the normal and the E-field pre-activated r2KCHA were obtained by cell refinement method using MDI Jade 6.0 software. The in situ synchrotron PXRD analysis was conducted at PD beamline, Australian Synchrotron, ANSTO with the E-field application during analysis to identify the instant response of zeolites to the E-field. The X-ray wavelength was calibrated to be 0.7732Å. The degassed zeolite powders were filled into a 0.7mm quartz capillary and then wax-sealed. Before analysis, the sealing wax was removed, and the capillarywasmounted on the PXRDdiffractometer to re-degas for 30min. During analysis, external E-fields with intensities of 100, 200, and 267 V/mm (with voltages of 120, 240, and 320V, respectively) were applied to the sample by the two-wire electrodes (1.2mm apart) parallel to the capillary. The PXRD data were collected over the tem- perature range of 203 ~ 493 K under vacuum, with a heating/cooling rate of 5 K/min. For each temperature, the sample was successively exposed to the E-fields with the voltages 0, 120, 240, 320, and 0V. For each setpoint, the samplewas held for 60 s and then scanned for 240 s. Ab initio density functional theory calculation The periodic density functional theory (DFT) calculation was con- ducted using the Vienna ab initio Simulation Package (VASP) V5.0.438. The initial crystal structure of r2KCHA was a rhombohedral lattice obtained from the Rietveld results of synchrotron PXRD data9. One r2KCHA unit cell is composed of three double six-ring prisms or one and a half supercavities. The generalized gradient approximation (GGA) and the projector augmented waves (PAW) approach39 were employed for the DFT calculation. The geometry optimization was performed with the Perdew-Burke-Ernzerhof (PBE) exchange- correlation functional. A gamma point only k-point mesh was used for a single unit cell. The cut-off energy of the plane wave basis-set was 405 eV, and a Gaussian smearing of 0.01 eV was applied. The atomic positions of r2KCHA were optimized until the forces acting on atoms were below 0.015 eV/Å, as suggested by the previous studies40. The DFT-D3 functional (IVDW= 11)41 was applied to account for the van der Waals interaction. The structure of r2KCHAwhen the K+ was located at different sites was not directly calculated by relaxing the unit cell, because the incomplete plane-wave basis set with respect to cell volume change will arouse the error of “Pulay stress”42. To reduce the Pulay stress, we performed the calculations at different volumes by changing the unit cell volume manually. The energy profile of the relative total energy versus the unit cell volume was plotted to obtain the optimized unit cell volume with the minimum total energy. The energy profile for the cation moving from SIII’ to SI was calculated by the nudged elastic band method43. The spontaneous polarization was calculated by using the Berry phase method44. The core electron density of the atoms was obtained by reconstructing all-electron charge density (the LAECHG tag) inVASP45, and the difference-electron density was calculated as the difference between the electron density of the r2KCHA-gas complex and the electron density of r2KCHA aug- mented by that of an isolated gas molecule28. Data availability The data generated in this study are provided in the Source Data file. Source data are provided with this paper. References 1. Chen, K. et al. Gating effect for gas adsorption in microporous materials-mechanisms and applications. Chem. Soc. Rev. 51, 1139–1166 (2022). 2. Remy, T. et al. Adsorption and separation of CO2 on KFI zeolites: effect of cation type and Si/Al ratio on equilibrium and kinetic properties. Langmuir 29, 4998–5012 (2013). 3. Shang, J. et al. Discriminative separation of gases by a “molecular trapdoor”mechanism in chabazite zeolites. J. Am. Chem. Soc. 134, 19246–19253 (2012). 4. Coudert, F. X. & Kohen, D. Molecular insight into CO2 “trapdoor” adsorption in zeolite Na-RHO. Chem. Mater. 29, 2724–2730 (2017). 5. Palomino, M., Corma, A., Jorda, J. L., Rey, F. & Valencia, S. Zeolite rho: A highly selective adsorbent for CO2/CH4 separation induced by a structural phase modification. Chem. Commun. 48, 215–217 (2012). 6. Lozinska, M. M. et al. Cation gating and relocation during the highly selective “trapdoor” adsorption of CO2 on univalent cation forms of zeolite RHO. Chem. Mater. 26, 2052–2061 (2014). 7. Shang, J., Li, G., Webley, P. A. & Liu, J. Z. A density functional theory study for the adsorption of various gases on a caesium-exchanged trapdoor chabazite. Comput. Mater. Sci. 122, 307–313 (2016). 8. Shang, J. et al. Temperature controlled invertible selectivity for adsorption of N(2) and CH(4) by molecular trapdoor chabazites. Chem. Commun. 50, 4544–4546 (2014). 9. Li, G. K. et al. Temperature-regulated guest admission and release in microporous materials. Nat. Commun. 8, 15777–15786 (2017). 10. Ghoufi, A., Benhamed, K., Boukli-Hacene, L. &Maurin,G. Electrically inducedbreathingof theMIL-53 (Cr)metal-organic framework.ACS Cent. Sci. 3, 394–398 (2017). 11. Schmid, R. An electric field induced breath for metal-organic fra- meworks. ACS Cent. Sci. 3, 369–371 (2017). 12. Kolesnikov, A. L., Budkov, Y. A., Möllmer, J., Kiselev, M. G. & Gläser, R. Metal–organic framework breathing in the electric field: a theo- retical study. J. Phys. Chem. C 123, 10333–10338 (2019). 13. Knebel, A. et al. Defibrillation of soft porous metal-organic frame- works with electric fields. Science 358, 347 (2017). 14. Chen, K. et al. Electrical regulation of CO2 adsorption in the metal- organic framework MIL-53. ACS Appl. Mater. Inter. 14, 13904–13913 (2022). 15. Ananchenko, G. S. et al. A molecular turnstile in para-octanoyl calix[4]arene nanocapsules. Chem. Commun. 7, 707–709 (2007). 16. Perreau, J. et al. Electricfield influenceon thebondingofCH2CHCN onto a nickel surface. Surf. Sci. 162, 776–781 (1985). 17. Durholt, J. P., Jahromi, B. F. & Schmid, R. Tuning the electric field response of MOFs by rotatable dipolar linkers. ACS Cent. Sci. 5, 1440–1448 (2019). 18. Namsani, S. & Yazaydin, A. O. Electric field induced rotation of halogenated organic linkers in isoreticular metal–organic frame- works for nanofluidic applications. Mol. Syst. Des. Eng. 3, 951–958 (2018). 19. Jensen, N. K. et al. Screening zeolites for gas separation applica- tions involving methane, nitrogen, and carbon dioxide. J. Chem. Eng. Data 57, 106–113 (2011). 20. Ilić, B. & Wettstein, S. G. A review of adsorbate and temperature- induced zeolite framework flexibility. Microporous Mesoporous Mater. 239, 221–234 (2017). Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 10 21. Lee, Y. et al. Non-framework cation migration and irreversible pressure-induced hydration in a zeolite. Nature 420, 485 (2002). 22. Kuznicki, S. M. et al. A titanosilicatemolecular sieve with adjustable pores for size-selective adsorption of molecules. Nature 412, 720–724 (2001). 23. Goetze, J., Yarulina, I., Gascon, J., Kapteijn, F. & Weckhuysen, B. M. Revealing lattice expansion of small-pore zeolite catalysts during the methanol-to-olefins process using combined operando x-ray diffraction and UV-Vis spectroscopy. ACS Catal. 8, 2060–2070 (2018). 24. Gavrilenko, P., Yakubovich, O., Massa, W. & Pekov, I. The crystal structure of chabazite-K. Crystallogr. Rep. 50, 544–553 (2005). 25. Smith, L. J., Eckert, H. &Cheetham,A. K. Potassiumcation effects on site preferences in themixed cation zeolite Li,Na−chabazite.Chem. Mater. 13, 385–391 (2001). 26. Saxton, C.G., Kruth, A., Castro,M.,Wright, P. A. &Howe, R. F. Xenon adsorption in synthetic chabazite zeolites. Microporous Mesopor- ous Mater. 129, 68–73 (2010). 27. Maeda, Y. et al. Improper ferroelectricity in stuffed aluminate sodalites for pyroelectric energy harvesting. Phys. Rev. Appl. 7, 034012 (2017). 28. Shang, J. et al. Adsorption of CO2, N2, and CH4 in Cs-exchanged chabazite: a combination of van der waals density functional theory calculations and experiment study. J. Chem. Phys. 140, 084705 (2014). 29. Rowsell, J. L., Spencer, E. C., Eckert, J., Howard, J. A. & Yaghi, O. M. Gas adsorption sites in a large-pore metal-organic framework. Sci- ence 309, 1350–1354 (2005). 30. Hudson, M. R. et al. Unconventional, highly selective CO2 adsorp- tion in zeolite SSZ-13. J. Am. Chem. Soc. 134, 1970–1973 (2012). 31. Krishna, R. & van Baten, J. M. Segregation effects in adsorption of CO2-containing mixtures and their consequences for separation selectivities in cage-type zeolites. Sep. Purif. Technol. 61, 414–423 (2008). 32. Guo, P. et al. A zeolite family with expanding structural complexity and embedded isoreticular structures. Nature 524, 74–78 (2015). 33. Zhao, J. et al. Nitrogen rejection from methane via a “trapdoor” K- ZSM-25 zeolite. J. Am. Chem. Soc. 143, 15195–15204 (2021). 34. Li G., May E. F.,Webley P. A., Huang S. H.-W. &Chan K. I. Method for gas separation. PCT patent WO2016094929A1 (2017). 35. Hu, G., May, E. F. & Li, K. G. We commercialized a methane capture technology in ten years—here’s how. Nature 604, 242–245 (2022). 36. Zhao, J. et al. Li+/ZSM-25 zeolite as a CO2 capture adsorbent with high selectivity and improved adsorption kinetics, showing CO2- induced framework expansion. J. Phys. Chem. C 122, 18933–18941 (2018). 37. Fagerlund, G. Determination of specific surface by the bet method. Mat.ériaux et. Constr. 6, 239–245 (1973). 38. Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996). 39. Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994). 40. Göltl, F. & Hafner, J. Alkane adsorption in Na-exchanged chabazite: the influence of dispersion forces. J. Chem. Phys. 134, 064102 (2011). 41. Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-d) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104 (2010). 42. Vanpoucke, D. E. P., Lejaeghere, K., Van Speybroeck, V., Waroquier, M. & Ghysels, A. Mechanical properties from periodic plane wave quantum mechanical codes: the challenge of the flexible nanoporous MIL-47 (V) framework. J. Phys. Chem. C 119, 23752–23766 (2015). 43. Jónsson H., Mills G. & Jacobsen K. W. Nudged elastic bandmethod for finding minimum energy paths of transitions. Classical and QuantumDynamics inCondensedPhaseSimulations. (eds. Berne, B. Ciccotti, G. & Coker D.) 385–404 (World Scientific, 1998). 44. King-Smith, R. D. & Vanderbilt, D. Theory of polarization of crystal- line solids. Phys. Rev. B Condens. Matter 47, 1651–1654 (1993). 45. Li, L. Effects of the interlayer interaction and electric field on the bandgap of polar bilayers: a case study of Sc2CO2. J. Phys. Chem.C 120, 24857–24865 (2016). Acknowledgements This work is sponsored by the Australia Research Council (no. DP190101336, G.K.L., P.A.W., and R.Q.S.). It was undertaken on the Powder X-ray Diffraction beamline at the Australian Synchrotron, part of ANSTO. The simulation was carried out at the Shanxi Supercomputing Center of China, and the calculations were performed on TianHe-2. Author contributions G.K.L. and K.C. conceived the idea. K.C. and Z.Y. conducted the experiments and analyzed the data. K.C. and S.H.M. conducted the computational simulation and analyzed the data. R.S. and Q.G. helped with characterizations. K.C. drafted the manuscript. This project was supervised by G.K.L., P.A.W., and R.Q.S. Competing interests The authors declare no competing interests. Additional information Supplementary information The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-41227-4. Correspondence and requests for materials should be addressed to Paul A. Webley or Gang Kevin Li. Peer review information Nature Communications thanks Alexander Knebel, Rochus Schmid and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available. Reprints and permissions information is available at http://www.nature.com/reprints Publisher’s note Springer Nature remains neutral with regard to jur- isdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/. © The Author(s) 2023 Article https://doi.org/10.1038/s41467-023-41227-4 Nature Communications | (2023) 14:5479 11 https://doi.org/10.1038/s41467-023-41227-4 http://www.nature.com/reprints http://creativecommons.org/licenses/by/4.0/ http://creativecommons.org/licenses/by/4.0/ Regulating adsorption performance of zeolites by pre-activation in electric fields Results E-field pre-activation of zeolites Gas adsorption in r2KCHA regulated by E-field pre-activation Framework expansion of r2KCHA induced by E-fields Cation relocation in r2KCHA induced by E-fields Gas admission regulated by the E-field induced cation relocation Demonstration for other zeolites The effect of E-field pre-activation on gas separation Discussion Methods Zeolite synthesis Sample degassing E-field pre-activation of zeolites Gas adsorption measurements Breakthrough gas separation Surface area analysis Powder X-ray diffraction analysis Ab initio density functional theory calculation Data availability References Acknowledgements Author contributions Competing interests Additional information