Intermolecular interactions governing adsorption of bulky molecules on ZIF-8: Insight from adsorption kinetics of benzene and 6-membered ring alicyclic hydrocarbons

Intermolecular interactions governing adsorption of bulky molecules on ZIF-8: Insight from adsorption kinetics of benzene and 6-membered ring alicyclic hydrocarbons | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Intermolecular interactions governing adsorption of bulky molecules on ZIF-8: Insight from adsorption kinetics of benzene and 6-membered ring alicyclic hydrocarbons Takahiro Ueda, Yuta Yamada, Kota Fujii, Ryota Mihara, Taku Iiyama, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-7276006/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 19 Mar, 2026 Read the published version in Adsorption → Version 1 posted 11 You are reading this latest preprint version Abstract This study examined the adsorption kinetics of benzene and six-membered ring alicyclic hydrocarbons on Zeolitic imidazolate framework-8 (ZIF-8) to elucidate the interactions between the adsorbed molecules and linkers on the adsorption behavior of bulky molecules. The temperature dependence of the Fickian diffusion coefficient was analyzed to determine the activation entropy and energy of diffusion. The significant negative values observed for the activation entropy, from − 172 to − 217 J K − 1 mol − 1 , indicate a restriction on the degree of freedom of the adsorbate molecules within the diffusion transition state, thereby reducing the diffusion coefficients. The ZIF-8 differentiates between 1,3- and 1,4-cyclohexadiene based on kinetic mechanisms, indicating its potential as an isomeric molecular sieve. The molecular orientations of the adsorbate molecules traversing the ZIF-8 aperture were investigated using Monte Carlo simulations. Order parameter analysis revealed that the molecular orientation was constrained to minimize the energy barrier with passage through the aperture. Fourier-transform infrared (FT-IR) spectroscopy and 13 C cross-polarization magic angle spinning (CP/MAS) nuclear magnetic resonance (NMR) measurements supported the CH-π interactions between the adsorbate molecules and 2-methylimidazolate linkers. This interaction maintained the orientation of the adsorbate molecules, facilitating their passage through the ZIF-8 apertures. CH-π interactions influence molecular orientation control in the adsorption and diffusion of large molecules containing π electrons in porous materials with organic ligands. Adsorption kinetics Molecular sieving Diffusion Activation entropy CH-π interaction Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 1 Introduction The growing interest in hydrogen as a prospective energy source has led to increased demand for hydrogen as a renewable resource in recent years. To establish a dependable hydrogen supply, there is an urgent need to construct large-scale storage facilities, and economically viable transportation systems [ 1 – 3 ]. The organic hydride process is a promising strategy for hydrogen supply [ 4 , 5 ] involving the storage of hydrogen in saturated hydrocarbons such as methylcyclohexane, which can be dehydrogenated to produce unsaturated hydrocarbons like toluene for hydrogen extraction [ 6 , 7 ]. The efficiency of the hydrogen extraction process can be improved by continuously removing toluene from the reaction system. Therefore, the development of molecular sieve technologies capable of effectively separating toluene from methylcyclohexane–toluene mixtures while minimizing energy consumption is essential. Beyond the methylcyclohexane/toluene system, the advancement of molecular sieve materials for the selective separation and adsorption of saturated and unsaturated compounds in bulky hydrocarbons with high hydrogen densities is anticipated to further enhance the organic hydride method [ 8 – 12 ]. Zeolitic Imidazolate Framework 8 (ZIF-8), a type of metal-organic framework, is one of the leading candidates for such molecular sieve materials. The ZIF-8 is represented by the formula [Zn(C 4 H 5 N 2 ) 2 ] n , consisting of Zn 2+ ions and the 2-methylimidazolate anion (C 4 H 5 N 2 − , MeIm − ) [ 13 ]. It possesses a sodalite structure characterized by micropores with a diameter of 11.4 Å, interconnected by eight 6-membered ring (6MR) openings measuring 3.4 Å in diameter [ 14 ]. This material exhibits significant potential for gas adsorption and separation applications [ 15 , 16 ], which can be attributed to its exceptional molecular sieving properties stemming from its flexible framework. The rotational motion of the bridging ligand modifies the apparent size of the 6MR opening, allowing the capture of molecules larger than 3.4 Å [ 17 , 18 ]. This phenomenon, referred to as "gate adsorption," results from the framework's flexibility [ 19 – 21 ]. The unique molecular sieving behavior of ZIF-8 is characterized by a nearly linear relationship between the activation energy for gas molecular diffusion and the size of the adsorbed molecules [ 18 ], although the effective diffusion coefficient derived from the adsorption rate decreases with the kinetic diameter of the adsorbed molecules. This diffusion behavior is linked to the swing effect, which is induced by the torsional motion of the imidazole ligand around the Zn–Im–Zn axis [ 22 – 24 ]. Significant fluctuations in the bridging ligands have been confirmed using various spectroscopic techniques, including far-infrared absorption spectroscopy [ 25 ], inelastic neutron scattering [ 26 ], terahertz spectroscopy [ 27 ], and solid-state NMR [ 28 , 29 ] methods. However, bulky molecules that exhibit significant molecular shape anisotropy may adopt an energetically favorable orientation for navigating through the 6MR aperture in an orientation that minimizes the molecular cross-sectional area. In previous studies, it was observed that for benzene and its derivatives (toluene and xylene), the minimum molecular cross-section rather than the average molecular diameter predominantly influenced the molecular orientation of the adsorbed molecules as they traversed the 6MR aperture [ 30 , 31 ]. Recently, Poryvaev et al. reported that benzene was adsorbed onto ZIF-8 with an adsorption selectivity of approximately 270% from a benzene/cyclohexane solution containing > 90% cyclohexane [ 32 ]. Their findings indicated that the interaction between benzene and cyclohexane with the 6MR aperture of ZIF-8 was governed by the molecular thickness of the adsorbed molecules, suggesting that the molecular orientation was critical for the passage of benzene and cyclohexane through the 6MR aperture. Thus, for adsorbed molecules with bulky and highly anisotropic geometries, both the intermolecular interactions with the 6MR aperture (enthalpic contributions) and the molecular orientation during passage through the aperture (entropic contributions) were expected to significantly influence their diffusion behavior. To sustain a favorable molecular orientation for passage through the 6MR aperture, highly directional interactions likely govern the orientation of the adsorbed molecules, as opposed to isotropic intermolecular interactions such as dispersion forces. Identifying the intermolecular interactions that contribute to the maintenance of molecular orientation is anticipated to enhance the theoretical understanding of the adsorption behavior of bulky molecules on ZIF-8, particularly regarding diffusion phenomena. Potential candidates for such interactions include hydrogen bonds, tetrel bonds, CH-π interactions, and halogen-π interactions [ 33 , 34 ]. We hypothesized that the double bond of 2-methylimidazolate, the bridging ligand of ZIF-8, would facilitate CH-π interactions. This study aimed to elucidate the molecular-level adsorption mechanism of bulky hydrocarbons in ZIF-8. Initially, we focused on the adsorption kinetics of benzene and several 6-membered ring alicyclic hydrocarbons on ZIF-8, with an emphasis on evaluating the influence of intermolecular interactions on diffusion behavior. The selected alicyclic compounds (1,3-cyclohexadiene, 1,4-cyclohexadiene, cyclohexene, and cyclohexane, summarized in Scheme S1 in supporting information) had a similar projected molecular diameter to benzene but differed in molecular thickness and π-electron content (see Table S1 , in SI). These variations facilitate the analysis of the diffusion coefficients and activation parameters in relation to the molecular dimensions and the presence of π-electrons. This study employed Monte Carlo simulations to investigate the critical process of molecules traversing the 6MR aperture of ZIF-8, with particular emphasis on the orientation of bulky molecules during this process. Additionally, Fourier-transform infrared absorption spectroscopy and 13 C cross-polarization magic-angle sample spinning (CPMAS) NMR techniques were utilized to assess adsorbate-adsorbent interactions, aiming to identify specific intermolecular interactions that influence molecular orientation as the adsorbate molecules navigate through the 6MR aperture. These research strategies collectively offer molecular-level understanding of the adsorption mechanism of bulky hydrocarbon molecules on ZIF-8. 2 Experimental 2.1 Chemicals The ZIF-8 (Basolite Z1200, BASF) and 1,3-cyclohexadiene (97% purity) were obtained from Sigma-Aldrich. Benzene (purity > 99.5%), cyclohexane (purity > 99.5%), 1,4-cyclohexadiene (purity > 95.0%), and cyclohexene (purity > 97.0%) were purchased from Fujifilm Wako Pure Chemicals Co., Ltd. The reagents utilized were Wako Special Grade for benzene and cyclohexane and Wako 1st Grade for 1,4-cyclohexadiene and cyclohexene. All the reagents were used without further purification. For this study, benzene, 1,3-cyclohexadiene, 1,4-cyclohexadiene, cyclohexene, and cyclohexane are referred to as BZ, 1,3-CHD, 1,4-CHD, CHEN, and CHEX, respectively. 2.2 Characterization 2.2.1 Nitrogen adsorption isotherm Prior to use, ZIF-8 was activated by heating under vacuum (0.1 Torr) at 150°C for 12 h. Nitrogen adsorption isotherms were recorded using a Gemini 2375 apparatus (Shimadzu Corporation) to determine the pore volume and specific surface area of ZIF-8. The isotherms were recorded over a relative pressure ( p / p 0 ) range of 0.0001–0.98, where p 0 represents the saturated vapor pressure of nitrogen at its boiling point of liquid nitrogen (77 K). The porosity of ZIF-8 was evaluated from the nitrogen adsorption isotherm (Figure S1 , Supplementary Information), which exhibited a characteristic Type I (a) pattern according to the IUPAC classification. The BET surface area was calculated to be 1074 m 2 g –1 , with a corresponding pore volume of 0.61 cm 3 g –1 . 2.2.2 Scanning electron microscopy (SEM) Scanning electron microscopy (SEM) was performed using a field-emission scanning electron microscope SU6600 (Hitachi High-Technologies Co.) to analyze the size distribution and morphology of the crystallites. The SEM images are shown in Figure S2 (Supplementary Information). The average particle diameter, based on the crystalline size distribution (Figure S3 ), was determined to be 450 nm, corresponding to a radius of 225 nm. 2.2.3 Fourier transferred infrared (FT-IR) absorption spectroscopy Fourier-transform infrared (FT-IR) measurements were performed using an FT/IR-4600AC spectrometer (JASCO Co.) over the wavenumber range of 400–4000 cm –1 at room temperature with a resolution of 1 cm –1 . To prevent the amorphization of the ZIF-8 sample due to uneven pressure during mortar grinding, KBr pellet samples were prepared by compressing the powder sample between two KBr plates. Spectra were recorded in transmittance units and subsequently converted to absorbance values. 2.2.4 Solid-state 13 C NMR spectroscopy High-resolution solid-state 13 C NMR spectra were obtained using a Bruker Avance III 400WB spectrometer equipped with a superconducting magnet (9.4 T), which provided resonance frequencies of 400.25 and 100.64 MHz for the 1 H and 13 C nuclei, respectively. Free induction decay (FID) signals for 13 C nuclei were acquired using a single-pulse sequence with 1 H decoupling under magic-angle sample spinning (MAS) conditions and a 1 H decoupling MAS cross-polarization (CP) pulse sequence. The powdered samples were packed in a 4 mm-diameter cylindrical zirconia rotor and spun at 10 kHz for MAS. Adamantane was used as an external standard for the chemical shift, with a CH 2 peak referenced at 38.5 ppm (relative to tetramethylsilane). 2.3 Measurement of time dependence of adsorption uptake The ZIF-8 powder sample underwent pretreatment (heating at 150°C in a vacuum of 0.1 Torr for 12 h). The samples with an average weight of 7–15 mg were placed in an aluminum pan (5 mm × 3 mm depth). The aluminum pan was placed in a hermetically sealed vessel containing 3 mL of the appropriate solvent. The presence of air in the container had a minimal effect on solvent vapor adsorption, as nitrogen (N 2 ) and oxygen (O 2 ) behave as supercritical gases at ambient temperatures. The adsorbates used in this study were BZ, 1,3-CHD, 1,4-CHD, CHEN, and CHEX. The vessel was affixed to a base constructed from thermally insulating styrofoam and placed within a temperature-controlled enclosure with an accuracy of ± 0.5°C. Vapor adsorption was conducted at 30, 40, and 50°C under saturated vapor pressures of the solvent. The temperature dependence of the saturated vapor pressure is shown in Figure S4 (Supplementary Information). At a constant temperature, BZ, 1,3-CHD, and CHEX exhibited similar vapor pressures within ± 1 kPa, whereas CHEN and 1,4-CHD displayed slightly lower vapor pressures, with differences of 1.5 kPa for CHEN and 5 kPa for 1,4-CHD. Consequently, the influence of the vapor pressure of each solvent on the adsorption rate was expected to be negligible. The samples were weighed at intervals until a stable value was reached, with weight measurements exhibiting an accuracy of ± 5%. 2.4 Monte Carlo simulations A Monte Carlo simulation was performed using proprietary software at a constant temperature of 298 K under NVT conditions, with a minimum of 100,000 trials. The force field parameters —including Lennard-Jones parameters and partial charges—assigned to the atoms forming the 6MR aperture were consistent with those used in the molecular dynamics (MD) simulations of ZIF-8 conducted by Parkes et al. [ 35 , 36 ]. The configuration of the adsorbed molecule was determined by optimizing its structure using density functional theory (DFT) calculations. The Gaussian 16 software package [ 37 ] was employed, utilizing the B3LYP functional and 6-31G(d) basis set for the calculations. In the Monte Carlo simulations, the OPLS-AA force field parameters [ 38 , 39 ] were used to define the Lennard-Jones parameters and partial charges of the adsorbed molecules. The 6MR structure employed in the simulation was modeled by extracting one-eighth of the unit cells dimension ( a /2, a /2, a/ 2) from the ZIF-8 crystal (lattice constant a ). The vector from the lattice point (0, 0, 0) to (1/2, 1/2, 1/2) was designated as the Z -axis in the reference coordinate system. The origin of this reference coordinate system was set at the center of the 6MR plane (lattice point (1/4, 1/4, 1/4)). The direction toward lattice point (0, 0, 0) was defined as negative (–), and that toward (1/2, 1/2, 1/2) as positive (+). The 6MR structure is shown in Figure S5 . At a specified Z value, the mass center of the adsorbed molecule was displaced over a 4 Å × 4 Å plane centered on the reference coordinate axis. The most energetically favorable configuration was then determined using Monte Carlo simulation. The opening of the 6MR aperture was adjusted by systematically rotating the tilt angle of the 2-methylimidazole bridging ligand along its N-N axis. The tilt angle was defined as the angle between the molecular plane of the bridging ligand and the 6MR aperture plane. Variations in tilt were expressed as positive (+Δ θ ) or negative (–Δ θ ) deviations from the ligand’s original orientation of the crystal lattice. The adsorbed molecule was moved along the Z-axis in the reference coordinate system from + 4 Å to − 4 Å in 0.2 Å increments, and the potential energy between the adsorbed molecule and the 6MR for each Z value was evaluated by varying Δ θ . For the order parameter analysis, the simulations were carried out in the Z value from + 5 Å to − 5 Å. 3 Results and Discussion 3.1 Time dependence of adsorption uptake Figure 1 depicts the time-dependence of the adsorption uptake of individual hydrocarbons on ZIF-8 at 303 K, illustrating the relationship between the amount of adsorption and the square root of the adsorption time. Initially, all the vapor species exhibited a linear increase in the adsorption amount, which was proportional to the square root of time, suggesting that micropore diffusion was the primary mechanism governing the adsorption of individual solvent vapors on ZIF-8. These observations can be characterized using the Fickian diffusion model, which is applicable to spherical particles [ 40 ]. Assuming that the concentration of adsorbed molecules at the center of the particle is negligible (i.e., in the short diffusion time region), the model can be approximated as $$\:m\left(t\right)={m}_{\text{e}\text{x}\text{t}}+{m}_{\infty\:}\frac{6}{\sqrt{\pi\:}}\sqrt{\frac{{D}_{\text{f},\text{i}\text{n}\text{i}}t}{{r}_{\text{a}\text{v}\text{e}}^{2}}}$$ 1 where D f,ini denotes the Fickian diffusion coefficient, which is derived from the initial slope of the time-dependent adsorption curve; r ave represents the average radius of the spherical particles; m ∞ indicates the saturation adsorption amount; and m ext refers to the adsorption amount present on the surface of the particle [ 41 ]. The diffusion coefficient D f,ini can be calculated from the slope of the linear portion depicted in Fig. 1 , assuming that the saturation adsorption amount is known. The data indicated that BZ and 1,4-CHD achieved saturation during adsorption, whereas 1,3-CHD approached saturation. In contrast, the CHEN and CHEX did not reach saturation during this period. For convenience, the adsorption quantity recorded at the longest adsorption time was designated as the saturation adsorption quantity, m ∞ , at equilibrium. The D f,ini values for each adsorbed molecule were extracted from the linear segments, as shown in Fig. 1 . Figure S6 (SI) presents the time-dependent adsorption amounts at 313 and 323 K. Table 1 summarizes the D f,ini values obtained at 303, 313, and 323 K. Notably, despite being structural isomers with identical molecular weights, 1,3-CHD and 1,4-CHD exhibited distinct diffusion coefficients. This finding suggests that ZIF-8 can differentiate between structural isomers based on their respective adsorption behavior through kinetic mechanisms. This observation is particularly significant, as it implies that ZIF-8 can potentially serve as a molecular sieve for structural isomers—an ability with important implications for separation applications. Table 1 Diffusion coefficients of adsorbate molecule determined from the initial slope of adsorption uptake curve Adsorbate D f,ini / m 2 s − 1 303 K 313 K 323 K BZ 1.3 (1) × 10 ‒19 2.2 (1) × 10 ‒19 2.8 (1) × 10 ‒19 1,4-CHD 5.6 (1) × 10 ‒20 7.9 (1) × 10 ‒20 9.5 (1) × 10 ‒20 1,3-CHD 3.2 (1) × 10 ‒20 3.5 (1) × 10 ‒20 5.4 (1) × 10 ‒20 CHEN 1.2 (1) × 10 ‒20 1.8 (1) × 10 ‒20 2.2 (1) × 10 ‒20 CHEX (chair) 8.2 (1) × 10 ‒21 1.1 (1) × 10 ‒20 1.7 (1) × 10 ‒20 The member in parentheses presents the experimental error. When larger molecules are adsorbed onto ZIF-8, they must navigate through the 6MR aperture, which consists of six zinc (Zn) atoms and 2-methylimidazole bridging ligands. A previous study established that, for benzene and its derivatives adsorbed on ZIF-8 [ 30 ], traversal through the 6MR aperture represents the rate-limiting step in the diffusion process. Furthermore, the movement of bulky molecules through the 6MR aperture is primarily governed by configurational diffusion, which is characterized by a restricted molecular orientation and rotation. Consequently, the temperature dependence of D f,ini should be assessed using the following equation derived from Eyring's transition state theory [ 42 ] rather than the traditional Arrhenius equation: $$\:{D}_{\text{f},\text{i}\text{n}\text{i}}=\text{e}{\lambda\:}^{2}\frac{kT}{h}{\text{e}}^{\varDelta\:{S}^{\ddagger}/R}{\text{e}}^{{-E}_{a}/RT}$$ 2 where E a signifies the diffusion activation energy, Δ S ‡ represents the diffusion activation entropy, h denotes Planck's constant, k is Boltzmann's constant, and λ indicates the distance between the molecular jumps. During the diffusion process, molecules jump between adjacent lattices via 6MR apertures with λ = 3 1/2 a /2. The lattice constant a is 1.701 nm at 298 K [ 22 ], resulting in a λ of 1.473 nm. Figure 2 illustrates a logarithmic representation of D f,ini ꞏ T ‒1 versus the reciprocal temperature (Eyring plot). The slope and intercept of the graph correspond to the activation energy and entropy, respectively. The activation parameters for each adsorbed molecule are listed in Table 2 . Table 2 Diffusion coefficients at 303 K and activation parameters determined from the temperature dependence of the diffusion coefficients. Adsorbate σ average /nm E a / kJ mol − 1 Δ S ‡ / J K − 1 mol − 1 BZ 0.525 28 (2) ‒172 (5) 1,4-CHD 0.537 19 (2) ‒208 (10) 1,3-CHD 0.538 18 (2) ‒217 (10) CHEN 0.548 24 (2) ‒205 (6) CHEX (chair) 0.558 27 (1) ‒200 (3) ‒ Figure 3 presents the molecular size dependence of D f,ini at 303 K and the activation parameters E a and ‒Δ S ‡ for individual hydrocarbons. The semi-logarithmic plot of D f,ini shows a linear decrease with increasing molecular size, whereas the standard plots of E a and ‒Δ S ‡ exhibited a V-shaped dependence on the molecular size. Specifically, E a reached its minimum for 1,4-CHD and 1,3-CHD, which possessed intermediate molecular sizes, and attained its maximum values for BZ and CHEX, representing the smallest and largest molecular sizes, respectively. This observation suggests that the passage of adsorbed molecules through 6MR aperture is governed by the combination of molecular size and specific host-guest interactions at the 6MR aperture site. The activation energy E a can be lowered by stabilizing the transition state or raised by stabilizing the adsorbed molecules within the pores. A lower E a enhances the likelihood of adsorbed molecules overcoming the diffusion energy barrier at a specific temperature. Conversely, Δ S ‡ exerts an opposing influence on D f,ini . The substantial negative activation entropy indicates a pronounced restriction on the degrees of freedom of the adsorbed molecule within the diffusion transition state, which leads to a reduced diffusion frequency. In BZ, the effect of Δ S ‡ —despite having the smallest absolute value—surpasses that of the largest E a , resulting in the highest D f,ini value among the compounds. For 1,4-CHD and 1,3-CHD, the E a values are less than 20 kJ mol ‒1 ; however, their Δ S ‡ values reach up to ca. 210 J K ‒1 mol ‒1 . The competing effects of these opposing factors may result in the contribution of activation entropy surpassing that of activation energy, thereby causing a decline in D f,ini as the molecular size increases. In BZ, 1,4-CHD, and 1,3-CHD, E a and Δ S ‡ exhibited a trade-off relationship. In contrast, for CHEN and CHEX, both the high E a values and the relatively large absolute values of Δ S ‡ decrease the D f,ini values. In the following sections, we examine the pronounced effects of molecular orientation on the diffusion of adsorbed molecules through the 6MR and identify the specific intermolecular interactions that significantly influence the E a . 3.2 Monte Carlo simulation for 6MR passage of adsorbate molecules Monte Carlo simulations were employed to elucidate the potential energy profile and molecular orientation of the adsorbed molecules as they traverse the 6MR aperture. In this context, the 6MR aperture, characterized by a tilted angle (Δ θ = 0) as derived from the crystal structure, presents an activation energy required for molecular passage with activation energies ranging from several hundred to several thousand kJ mol ‒1 . Consequently, it is nearly impossible for molecules to pass through the aperture using only the thermal energy available at room temperature. This observation implies that the 6MR aperture must undergo enlargement to facilitate the passage of target molecules. To quantify this, we evaluated the activation energy required for molecular passage as a function of the tilted angle of the bridging ligand by systematically varying this angle. By employing polynomial interpolation on the acquired data, we identified a tilted angle that aligned with the experimentally measured values. Figure 4 illustrates the energy profile associated with passage through the 6MR aperture, corresponding to the experimental activation energy for each adsorbate, as determined by the time dependence of the adsorption uptake. Two minima were observed on either side of the 6MR aperture. The energy minimum at the CH 3 neck ( Z /Å ~ − 2) was deeper than that at the CH neck ( Z /Å ~ +3). This finding suggests that the passage of molecules trapped in the CH 3 neck through the 6MR aperture is energetically unfavorable. Consequently, the passage of molecules through the 6MR aperture is only considered when they approach from the CH neck side [ 30 ]. In this case, the energy barrier encountered by an adsorbed molecule navigating through the 6MR aperture is defined as the difference (Δ E = E CH,max − E CH,min ) between the maximum energy ( E CH,max ) located near the CH neck (+ 0.5 ≤ Z /Å ≤ +1) and the minimum energy ( E CH,min ) at the CH neck ( Z /Å ~ +3). Notably, in the case of CHEX, an additional energy maximum was detected on the outside of the CH neck ( Z /Å ~ +2). The presence of the dual maxima in the energy profile aligns with the findings from previous studies. This behavior supports the use of a model that incorporates pre-equilibrium adsorption at the CH neck, represented as: A(g) + M(s) ⇌ AM*(s) ⇌ AM(s) [ 30 ]. In this model, the apparent activation energy can be expressed as Δ E = E a,a – E a,d + E a,in , where E a,a denotes the energy barrier for pre-adsorption at the 6MR aperture, E a,d is the energy barrier for the desorption of the pre-adsorbed molecule, and E a,in is the energy barrier for passage through the CH neck. The activation energy evaluated by MC simulation, the resulting tilted angle of the bridging ligand, and the aperture opening diameter are listed in Table 3 . The correlation between the tilted angle and the opening diameter of the 6MR aperture, which reproduces the experimental values, appears to align more closely with MIN-1 than with MIN-2, as proposed by Webster et al. [ 43 ] (see Scheme S1 and Table S1 in the SI). This suggests that the energy barrier for the passage of the adsorbed molecule through the 6MR aperture is significantly influenced by the cross-sectional area of the adsorbed molecule, indicating that the molecular orientation is highly constrained during this process. This conclusion was consistent with the findings of our previous study. Table 3 Experimental and simulated activation energies of adsorbed molecules passing through the 6MR aperture, tilted angles of the bridging ligand, and average diameters of the 6MR aperture. Adsorbate Δ E / kJ mol − 1 Δ θ c / deg. d p,ave / nm MIN-1 b) /nm MIN-2 b) /nm BZ 27.6 17.0 0.47 0.3277 0.6628 1,4-CHD 19.3 21.1 0.49 0.3845 0.6612 1,3-CHD 17.5 19.9 0.49 0.3845 0.6612 CHEN 24.2 20.55 0.49 0.4414 0.6595 CHEX (chair) 27.0 a) 24.05 0.52 0.4982 0.6580 a. This value corresponds to the apparent activation energy when adsorption pre-equilibration occurs: Δ E = E a,a – E a,d + E a,in . b. Critical dimensions of adsorbate molecule for entry into zeolite pores. In the slit-shaped pores, the size of the adsorption in the minimum dimension, MIN-1, was determined. In the cylindrical pores, the size of the molecule in two directions is considered: the minimum dimension, MIN-1, and the next smallest dimension, MIN-2. To further analyze the orientation of the adsorbed molecules, we defined two order parameters: \(\:{S}_{\perp\:}=\frac{1}{2}\left(3⟨{\text{cos}}^{2}\alpha\:⟩-1\right)\) (3a) \(\:{S}_{\parallel\:}=\frac{1}{2}\left(3⟨{\text{cos}}^{2}\beta\:⟩-1\right)\) (3b) where α and β defined the angles formed by the vectors perpendicular and parallel to the molecular plane, respectively, with respect to the perpendiculars passing through the center of the 6MR aperture. The characterization of these vectors and their orientation within the reference coordinate system are shown in Figure S7 (see SI). Figure 5 presents the order parameter for each molecule as it traverses the 6MR aperture. An animation depicting the molecular passage through the 6MR is provided in the SI. When the molecular orientation is isotropically distributed, the order parameter is zero. However, as each molecule passed through the aperture, both S ⊥ and S || deviated from zero, indicating that the orientation of the molecules was constrained by the aperture. Notably, the negative value of S ⊥ suggests that the molecules navigate the aperture in an orientation in which the molecular plane is nearly perpendicular to the aperture plane. This behavior was exhibited by all the adsorbate molecules. This indicated that the passage of adsorbed molecules through the 6MR aperture was accompanied by significant orientational constraints. In other words, the translational and rotational degrees of freedom of the adsorbed molecules were considerably restricted during their traversal through the 6MR aperture. Consequently, the observed significant negative activation entropy can be attributed to the transition of molecules in the vapor to a diffusion transition state, in which the molecular orientation is markedly constrained at the 6MR aperture. Conversely, S || represents the variation in the rotational orientation of the molecules around an axis perpendicular to the molecular plane and exhibited different behaviors for each adsorbate molecule in the range − 3 ≤ Z /Å ≤ +3. For BZ, 1,4-CHD, CHEN, and CHEX, the values ranged from − 0.5 < S || < +0.5, whereas for 1,3-CHD, 0 < S || < +1. In BZ and 1,4-CHD, S || exhibited a bipolar change centered at Z = 0, indicating that the molecule turned its orientation around an axis perpendicular to the molecular plane and then turned it back. In the CHEN and CHEX, the variation in S || showed a similar trend, although the range of the Z -axis was narrow. This feature implies that these molecules undergo twisting and wiggling around the axis perpendicular to the molecular plane when traversing the 6MR aperture (see the animation provided in the SI). In contrast, for 1,3-CHD, S || varied continuously from 0 to + 1 during the passage through 6MR, suggesting that the molecular orientation around the perpendicular axis changes progressively throughout the diffusion process. That is, the 1,3-CHD molecule rotates around the axis as it moves through the aperture (see the animation in the SI). This rotation is believed to minimize steric hindrance within the 6MR aperture by reducing friction between the molecule and the framework. This behavior differs from that of other molecules and may partly explain why 1,3-CHD exhibits a lower activation energy for diffusion than the other molecules. Distinguishing between 1,3-CHD and 1,4-CHD is crucial for understanding the factors influencing the diffusion coefficient of ZIF-8. The differences in the E a and S || behaviors observed in the MC simulations between 1,3-CHD and 1,4-CHD may reflect the slight differences in the molecular cross-sections and arrangement of π bonds within the molecules. 3.3 Analysis of adsorbate-adsorbent interactions by spectroscopies 3.3.1 FT-IR spectroscopy Figure 6 displays the symmetric C–H stretching bands associated with the 2-methylimidazole group. The full wavenumber regions of the infrared spectra are shown in Figure S8 (SI). Notably, no shift was observed in the band upon cyclohexane adsorption, while a pronounced redshift was observed over the wavenumber range of 1–2 cm ‒1 upon benzene and cyclic alkene adsorption. This shift was more pronounced for 1,3-CHD, 1,4-CHD, and BZ than for CHEN. This observation suggests the elongation of C–H bond, which is attributed to the attractive interaction between the π-electrons of the adsorbed molecules and the C–H bond of 2-methylimidazole [ 31 ]. One plausible attractive interaction is the CH-π interaction, which has been reported for various organic compounds, metal complexes, and molecular assemblies [ 33 ]. It has been established that CH-π bonding induces stretching of the C–H bonds, resulting in a redshift of the CH stretching band [ 44 , 45 ]. Additionally, the protons of the methyl group may function as proton donors in CH-π interactions, thereby leading to an anticipated redshift of the CH symmetric stretching band in the methyl group of the 2-methylimidazolate moiety [ 31 ]. However, the symmetric CH stretching band of the methyl group overlaps with the CH stretching bands of the adsorbate molecules. Deconvolution of the infrared bands potentially confirmed the redshift of the CH symmetric stretching band of the methyl group. However, the accuracy of the resultant redshift from this deconvolution was insufficient for a meaningful discussion of the minor wavenumber shifts observed in this study. The elongation of the C–H bonds and the corresponding redshift of the CH symmetric stretching band were also qualitatively corroborated by the molecular assembly of 2-methylimidazole with benzene, 1,3-CHD, and 1,4-CHD, optimized using DFT calculations (Figure S9, Figure S10, and Table S3 ). Generally, intermolecular CH-π interactions occur within 3 Å [ 33 ] and exhibit a high degree of directionality, with the orientation of the C–H bond aligned with that of the π-orbital. Consequently, the C–H and π bonds tend to be oriented nearly perpendicular to each other. Following the formation of the CH-π bond, the adsorbed molecules were arranged such that their molecular planes were positioned at an angle that was nearly perpendicular to the 6MR aperture plane. This orientation facilitates the ingress of adsorbed molecules into the 6MR apertures. Thus, the CH-π interaction between the 2-methylimidazole ring and the adsorbed molecules is crucial for governing the molecular orientation necessary for the effective passage through the 6MR aperture for benzene and cyclic alkenes. 3.3.2 13 C CP/MAS NMR spectroscopy The 13 C CP/MAS NMR spectra of ZIF-8, both in its pure form and with various adsorbates are shown in Fig. 7 . Cross-polarization (CP) measurements are particularly effective for highlighting immobile molecular segments within solid matrices while attenuating signals from more mobile components. Consequently, the highly mobile adsorbed molecules appeared as minor signals in the spectra. The predominant peaks are attributed to the 2-methylimidazole moiety, suggesting that the ZIF-8 framework exhibited greater rigidity and reduced mobility than the adsorbate molecules. Specifically, the resonance lines at approximately 14, 124, and 151 ppm correspond to the methyl group, CH carbon of the imidazole ring, and quaternary carbon of the 2-methylimidazolate moiety, respectively. Figures 8 (a), 8(b), and 8(c) provide enlarged views of these spectral components. Notably, the adsorption of benzene and cyclic alkenes resulted in downfield shifts of 0.5–1.3 ppm in the resonance lines associated with the methyl group and CH carbon of the imidazole ring, whereas the adsorption of cyclohexane did not elicit such shifts. The resonance lines corresponding to quaternary carbons remained largely unchanged after adsorption. These observations suggest a significant interaction between the π-electrons of the adsorbed molecules and both the methyl group and the CH carbon of the 2-methylimidazolate moiety. Our previous studies indicated that both the methyl and CH protons of the imidazole ring can act as proton donors in CH-π bond. According to Pople's point magnetic dipole model [ 46 , 47 ], the shielding effect of carbon along a perpendicular line through the center of the molecular plane of the benzene ring is described by the equation Δ σ C = 51.96/ R 3 , where R represents the average distance between the carbon nucleus of interest and the molecular plane in angstroms (Å). For the optimized structure of the 2-methylimidazole-benzene complex derived from DFT calculations (see Figure S9), R was determined to be 3.27 Å, which yielded a shielding effect of approximately 1.5 ppm. In the case of CH-π bond, the C-H bond is anticipated to elongate owing to the attraction of the proton by the π electrons, resulting in an anti-shielding effect on the 13 C nucleus. Scheiner reported an anti-shielding effect of 200 ppm/Å for the tetrel bond of a methyl group [ 48 ]. To account for the observed downfield shift, an average elongation of approximately 0.01–0.015 Å in the C-H bond is necessary, which is consistent with the expected bond-length increase associated with a strong CH-π bond. In summary, the downfield shifts in the resonance lines resulting from the adsorption of benzene and cyclic alkenes can be attributed to the CH-π bond between the CH protons of the methyl group and imidazole ring and π electrons of the adsorbed molecules. 3.4 Local structure of molecular passage through 6MR This study elucidates the intermolecular interactions that dictate molecular orientation during the adsorption of bulky molecules onto ZIF-8, particularly as these molecules navigate through 6MR apertures. The results suggest a specific configuration for benzene and six-membered-ring alicyclic hydrocarbons as they traverse the 6MR aperture (Fig. 9 ). The CH-π bond are pivotal in influencing the molecular orientation of benzene and cyclic alkenes during their passage through the aperture. As these molecules approach the CH bottleneck of the 6MR, they form CH-π bonds with the CH groups of the 2-methylimidazole ring at this bottleneck. This highly directional CH-π bonding results in an nearly perpendicular (T-shaped) alignment of the imidazole ring and adsorbed molecule. This arrangement ensured that the adsorbed molecules maintained a configuration that minimized their molecular cross-sectional areas (molecular planes) upon entering the 6MR aperture. After successfully passing through the aperture, the adsorbed molecule is stabilized by forming an additional CH-π bond with the methyl group on the methyl bottleneck side (Figure S10). In contrast, cyclic alkanes such as cyclohexane primarily engage in dispersion-force interactions. The optimal molecular orientation for enhancing these interactions occurs when the plane of the adsorbed molecule is parallel to the 2-methylimidazole ring, as this alignment maximizes the contact area between them. Consequently, a molecule approaching the CH bottleneck side of 6MR is captured by dispersion force interactions with the methyl group and quaternary carbon in the 2-methylimidazolate moieties. This results in the adsorbed molecule orienting its minimum molecular cross-section (molecular face) towards the 6MR aperture, thereby facilitating its entry. After passing through the 6MR aperture, the molecule stabilizes through dispersive force interactions with the three methyl groups on the methyl bottleneck side. Cyclohexene, being a cyclic alkene, also exhibits significant contributions from dispersion force interactions due to its substantial proportion of saturated hydrocarbon components, in addition to the CH-π interactions (see Figure S10). Finally, ZIF-8 exhibited remarkably different diffusion coefficients for 1.3-CHD and 1,4-CHD. This result is of particular interest because it highlights the potential of ZIF-8 as a molecular sieve for isomer separation. Further investigation is necessary to elucidate these detailed mechanisms. However, the isomer distinction may involve the symmetry of the phonon modes associated with the swing effect of the bridging ligands that constitute the 6MR, as well as the symmetry of the adsorbate molecule. Additionally, considering the aspect of configuration diffusion, this isomer distinction might be linked to the internal arrangement of π electrons through the CH-π bond, which controls the molecular orientation with respect to the 6MR aperture. 4 Conclusion In this study, we investigated the adsorption kinetics of benzene and six-membered-ring alicyclic hydrocarbons on ZIF-8, to clarify the interactions between the adsorbed molecules and linkers in relation to the adsorption behavior of larger molecules. Monte Carlo simulations were employed to elucidate the local structural dynamics of the molecules as they navigated through the six-membered ring aperture of ZIF-8. Additionally, the FT-IR and 13 C cross-polarization magic angle spinning nuclear magnetic resonance (CP/MAS NMR) were used to identify the adsorbate-adsorbate interactions that influence molecular orientation during transit through the aperture. The findings indicate that the diffusion of benzene and other six-membered ring hydrocarbons through the six-membered ring aperture is predominantly governed by a significant negative activation entropy ranging from − 172 to − 217 J K − 1 mol − 1 . This observation suggests a pronounced effect of the molecular orientation on the transit of adsorbate molecules through the apertures. Monte Carlo simulations further elucidated this phenomenon, revealing that the aperture diameter significantly affected the activation energy required for the molecular passage. Order parameter analysis revealed that the molecular orientation was tightly constrained and aligned in a manner that minimizes the energy barrier for traversing the aperture. Spectroscopic analyses employing FT-IR and 13 C CP/MAS NMR techniques identified a specific CH-π interaction between the π electrons of the adsorbed molecules and the methyl group and imidazole ring CH of the 2-methylimidazolate linker. For unsaturated cyclic hydrocarbons possessing π electrons, the CH-π interaction predominantly governs the molecular orientation during aperture transit. Conversely, saturated cyclic hydrocarbons are primarily influenced by dispersion force interactions with the methyl groups of the bridging ligand. These orientation-controlling interactions are anticipated to play a critical role in governing molecular adsorption and diffusion behaviors of other metal-organic frameworks (MOFs), porous coordination polymers (PCPs), and porous organic materials (POMs) composed of organic materials. Finally, the distinct diffusion coefficients of 1.3-CHD and 1,4-CHD are of particular interest. This indicated the potential of ZIF-8 as a molecular sieve for structural isomers. Declarations Acknowledgements The authors thank Dr. Naoya Inazumi and Dr. Yasuto Todokoro of the Analytical Instrument Facility, Graduate School of Science, Osaka University, for their helpful advice, guidance, and instructions regarding the solid-state NMR measurements. We would like to thank Editage (www.editage.jp) for English language editing. Author Contributions Takahiro Ueda: Conceptualization, methodology, validation, resources, software, formal analysis, visualization, writing–original draft, review, editing, supervision, funding acquisition, and project administration. Yuta Yamada: Methodology, investigation, data curation, and visualization. Kota Fujii: Investigation, DFT calculations, data curation, and visualization. Ryota Mihara: Investigation, Monte Carlo simulation, data curation, and visualization. Taku Iiyama: Methodology, investigation, software, data curation, validation, visualization, review, and editing. Yasutaka Hamada: Investigation, DFT calculations, data curation, and visualization. Mitsutaka Okumura: Methodology, software, and supervision. Funding This research was supported by JSPS KAKENHI Grant-in-Aid for Scientific Research (C) (Grant Numbers JP21K04979, JP24K8360). This work was the result of using research equipment shared in the MEXT Project to promote the public utilization of advanced research infrastructure (Program for supporting the construction of core facilities), Grant Numbers JPMXS0441200021, JPMXS0441200023, and JPMXS0441200024. Availability of data and materials The data presented in this article will be available from the corresponding author on reasonable request. Ethical Approval Not applicable. Conflicts of interest The authors declare no conflicts of interest. References Fan, L., Tu, Z., Chan, S.H.: Recent development of hydrogen and fuel cell technologies: A review. Energy Rep. 7 , 8421–8446 (2021). https://doi.org/10.1016/j.egyr.2021.08.003 Crandall, B.S., Brix, T., Weber, R.S., Jiao, F.: Techno-Economic Assessment of Green H 2 Carrier Supply Chains. Energy Fuels. 37 , 1441–1450 (2023). https://doi.org/10.1021/acs.energyfuels.2c03616 Salman, M.S., Rambhujun, N., Pratthana, C., Srivastava, K., Aguey-Zinsou, K.-F.: Catalysis in Liquid Organic Hydrogen Storage: Recent Advances, Challenges, and Perspectives. Ind. Eng. Chem. Res. 61 , 6067–6105 (2022). https://doi.org/10.1021/acs.iecr.1c03970 Swesi, Y., Kerleau, P., Pitault, I., Heurtaux, F., Ronze, D.: Purification of hydrogen from hydrocarbons by adsorption for vehicles application. Sep. Purif. Technol. 56 , 25–37 (2007). https://doi.org/10.1016/j.seppur.2006.09.023 Acharya, D., Ng, D., Xie, Z.: Recent Advances in Catalysts and Membranes for MCH Dehydrogenation: A Mini Review. Membranes. 11 , 955 (2021). https://doi.org/10.3390/membranes11120955 Meng, J., Zhou, F., Ma, H., Yuan, X., Wang, Y., Zhang, J.: A Review of Catalysts for Methylcyclohexane Dehydrogenation. Top. Catal. 64 , 509–520 (2021). https://doi.org/10.1007/s11244-021-01465-6 Chen, L., Verma, P., Hou, K., Qi, Z., Zhang, S., Liu, Y.-S., Guo, J., Stavila, V., Allendorf, M.D., Zheng, L., Salmeron, M., Prendergast, D., Somorjai, G.A., Su, J.: Reversible dehydrogenation and rehydrogenation of cyclohexane and methylcyclohexane by single-site platinum catalyst. Nat. Commun., 13, 1092 (2022). https://www.nature.com/articles/s41467-022-28607-y Zhokh, A., Serebrii, T., Strizhak, P.: Diffusion-controlled adsorption–desorption of C6 hydrocarbons by multiwall carbon nanotubes. Chem. Phys. Lett. 805 , 139935 (2022). https://doi.org/10.1016/j.cplett.2022.139935 Mukherjee, S., Sensharma, D., Qazvini, O.T., Dutta, S., Macreadie, L.K., Ghosh, S.K., Babarao, R.: Advances in adsorptive separation of benzene and cyclohexane by metal-organic framework adsorbents. Coord. Chem. Rev. 437 , 213852 (2021). https://doi.org/10.1016/j.ccr.2021.213852 Graczová, E., Molnár, D., Steltenpohl, P.: Separation of Methylcyclohexane-Toluene Mixture Using Imidazolium Ionic Liquids. Chem. Eng. Trans. 94 , 745–750 (2022). https://doi.org/10.3303/CET2294124 Valencia, A., Muñiz-Valencia, R., Ceballos-Magaña, S.G., Rojas-Mayorga, C.K., Bonilla-Petriciolet, A., González, J., Aguayo-Villarreal, I.A.: Cyclohexane and benzene separation by fixed-bed adsorption on activated carbons prepared from coconut shell. Environ. Technol. Innov. 25 , 102076 (2022). https://doi.org/10.1016/j.eti.2021.102076 Chen, C., Guan, H., Li, H., Zhou, Y., Huang, Y., Wei, W., Hong, M., Wu, M.: A Noncovalent π-Stacked Porous Organic Molecular Framework for Selective Separation of Aromatics and Cyclic Aliphatics. Angew Chem. Int. Ed. 61 , e202201646 (2022). https://doi.org/10.1002/anie.202201646 Phan, A., Doonan, C.J., Uribe-Romo, F.J., Knobler, C.B., O’Keeffe, M., Yaghi, O.Y.: Synthesis, Structure, and Carbon Dioxide Capture Properties of Zeolitic Imidazolate Frameworks. Acc. Chem. Res. 43 , 58–67 (2010). https://doi.org/10.1021/ar900116g Banerjee, R., Phan, A., Wang, B., Knobler, C., Furukawa, H., O’Keeffe, M., Yaghi, O.M.: High-Throughput Synthesis of Zeolitic Imidazolate Frameworks and Application to CO 2 Capture. Science. 319 , 939–943 (2008). https://doi.org/10.1126/science.1152516 Cheng, P., Hu, Y.H.: H 2 O-Functionalized Zeolitic Zn(2-methylimidazole) 2 Framework (ZIF-8) for H 2 Storage. J. Phys. Chem. C. 118 , 21866–21872 (2014). https://pubs.acs.org/doi/ 10.1021/jp507030g Böhme, U., Barth, B., Paula, C., Kuhnt, A., Schwieger, W., Mundstock, A., Caro, J., Hartmann, M.: Ethene/Ethane and Propene/Propane Separation via the Olefin and Paraffin Selective Metal–Organic Framework Adsorbents CPO-27 and ZIF-8. Langmuir. 29 , 8592–8600 (2013). https://doi.org/10.1021/la401471g Zhang, C., Lively, R.P., Zhang, K., Johnson, J.R., Karvan, O., Koros, W.J.: Unexpected Molecular Sieving Properties of Zeolitic Imidazolate Framework-8. J. Phys. Chem. Lett. 3 , 2130–2134 (2012). https://doi.org/10.1021/jz300855a Zhang, K., Lively, R.P., Zhang, C., Chance, R.R., Koros, W.J., Sholl, D.S., Nair, S.: Exploring the Framework Hydrophobicity and Flexibility of ZIF-8: From Biofuel Recovery to Hydrocarbon Separations. J. Phys. Chem. Lett. 4 , 3618–3622 (2013). https://doi.org/10.1021/jz402019d Fairen-Jimenez, D., Moggach, S.A., Wharmby, M.T., Wright, P.A., Parsons, S., Düren, T.: Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations. J. Am. Chem. Soc. 133 , 8900–8902 (2011). https://doi.org/10.1021/ja202154j Tian, T., Wharmby, M.T., Parra, J.B., Ania, C.O., Fairen-Jimenez, D.: Role of crystal size on swing-effect and adsorption induced structure transition of ZIF-8. Dalton Trans. 45 , 6893–6900 (2016). https://doi.org/10.1039/C6DT00565A Coudert, F.-X.: Molecular Mechanism of Swing Effect in Zeolitic Imidazolate Framework ZIF-8: Continuous Deformation upon Adsorption. ChemPhysChem. 18 , 2732–2738 (2017). https://doi.org/10.1002/cphc.201700463 Morris, W., Stevens, C.J., Taylor, R.E., Dybowski, C., Yaghi, O.M., Garcia-Garibay, M.A.: NMR and X-ray Study Revealing the Rigidity of Zeolitic Imidazolate Frameworks. J. Phys. Chem. C. 116 , 13307–13312 (2012). https://doi.org/10.1021/jp303907p Tanaka, H., Ohsaki, S., Hiraide, S., Yamamoto, D., Watanabe, S., Miyahara, M.T.: Adsorption-Induced Structural Transition of ZIF-8: A Combined Experimental and Simulation Study. J. Phys. Chem. C. 118 , 8445–8454 (2014). https://doi.org/10.1021/jp500931g Ania, C.O., García-Pérez, E., Haro, M., Gutiérrez-Sevillano, J.J., Valdés-Solís, T., Parra, J.B., Calero, S.: Understanding Gas-Induced Structural Deformation of ZIF-8. J. Phys. Chem. Lett. 3 , 1159–1164 (2012). https://doi.org/10.1021/jz300292y Ryder, M.R., Civalleri, B., Bennett, T.D., Henke, S., Rudić, S., Cinque, G., Fernandez-Alonso, F., Tan, J.-C.: Identifying the Role of Terahertz Vibrations in Metal-Organic Frameworks: From Gate-Opening Phenomenon to Shear-Driven Structural Destabilization. Phys. Rev. Lett. 113 , 215502 (2014). https://doi.org/10.1103/PhysRevLett.113.215502 Casco, M.E., Cheng, Y.Q., Daemen, L.L., Fairen-Jimenez, D., Ramos-Fernández, E.V., Ramirez-Cuesta, A.J., Silvestre-Albero, J.: Gate-opening effect in ZIF-8: the first experimental proof using inelastic neutron scattering. Chem. Commun. 52 , 3639–3642 (2016). https://doi.org/10.1039/C5CC10222G Tan, N.Y., Ruggiero, M.T., Orellana-Tavra, C., Tian, T., Bond, A.D., Korter, T.M., Fairen-Jimenez, D., Zeitler, J.A.: Investigation of the terahertz vibrational modes of ZIF-8 and ZIF-90 with terahertz time-domain spectroscopy. Chem. Commun. 51 , 16037–16040 (2015). https://doi.org/10.1039/C5CC06455D Kolokolov, D.I., Stepanov, A.G., Jobic, H.: Mobility of the 2-Methylimidazolate Linkers in ZIF-8 Probed by 2 H NMR: Saloon Doors for the Guests. J. Phys. Chem. C. 119 , 27512–27520 (2015). https://doi.org/10.1021/acs.jpcc.5b09312 Ueda, T., Nakai, M., Yamatani, T.: A solid-state 1 H-NMR study of the dynamic structure of ZIF-8 and its role in the adsorption of bulky molecules. Adsorption. 23 , 887–901 (2017). https://doi.org/10.1007/s10450-017-9902-5 Ueda, T., Yamatani, T., Okumura, M.: Dynamic Gate Opening of ZIF-8 for Bulky Molecule Adsorption as Studied by Vapor Adsorption Measurements and Computational Approach. J. Phys. Chem. C. 123 , 27542–27553 (2019). https://doi.org/10.1021/acs.jpcc.9b07239 Yagi, R., Ueda, T.: Substitution (CH 3 , Cl, or Br) effects of the imidazolate linker on benzene adsorption kinetics for the zeolitic imidazolate framework (ZIF)-8. Phys. Chem. Chem. Phys. 25 , 20585–20596 (2023). https://doi.org/10.1039/D3CP01662E Poryvaev, A.S., Yazikova, A.A., Polyukhov, D.M., Fedin, M.V.: Ultrahigh selectivity of benzene/cyclohexane separation by ZIF-8 framework: Insights from spin-probe EPR spectroscopy. Micropor Mesopor Mater. 330 , 111564 (2022). https://doi.org/10.1016/j.micromeso.2021.111564 Nishio, M.: CH/π hydrogen bonds in crystals. CrystEngComm. 6 , 130–158 (2004). https://doi.org/10.1039/B313104A Scheiner, S.: Origins and properties of the tetrel bond. Phys. Chem. Chem. Phys. 23 , 5702–5717 (2021). https://doi.org/10.1039/D1CP00242B Parkes, M.V., Demir, H., Teich-McGoldrick, S.L., Sholl, D.S., Greathouse, J.A., Allendorf, M.D.: Molecular dynamics simulation of framework flexibility effects on noble gas diffusion in HKUST-1 and ZIF-8. Micropor Mesopor Mater. 194 , 190–199 (2014). https://doi.org/10.1016/j.micromeso.2014.03.027 Zheng, B., Sant, M., Demontis, P., Suffritti, G.B.: Force Field for Molecular Dynamics Computations in Flexible ZIF-8 Framework. J. Phys. Chem. C. 116 , 933–938 (2012). https://doi.org/10.1021/jp209463a Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Scalmani, G., Barone, V., Petersson, G.A., Nakatsuji, H., Li, X., Caricato, M., Marenich, A.V., Bloino, J., Janesko, B.G., Gomperts, R., Mennucci, B., Hratchian, H.P., Ortiz, J.V., Izmaylov, A.F., Sonnenberg, J.L., Williams-Young, D., Ding, F., Lipparini, F., Egidi, F., Goings, J., Peng, B., Petrone, A., Henderson, T., Ranasinghe, D., Zakrzewski, V.G., Gao, J., Rega, N., Zheng, G., Liang, W., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Throssell, K., Montgomery, J.A. Jr., Peralta, J.E., Ogliaro, F., Bearpark, M.J., Heyd, J.J., Brothers, E.N., Kudin, K.N., Staroverov, V.N., Keith, T.A., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A.P., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Millam, J.M., Klene, M., Adamo, C., Cammi, R., Ochterski, J.W., Martin, R.L.: (2019). Morokuma, K., Farkas, O., Foresman, J. B., Fox, D. J.: Gaussian 16, Revision C.01, Gaussian, Inc., Wallingford CT Jorgensen, W.L., Tirado-Rives, J.: The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J. Am. Chem. Soc. 110 , 1657–1666 (1988). https://pubs.acs.org /doi/10.1021/ja00214a001 Jorgensen, W.L., Maxwell, D.S., Tirado-Rives, J.: Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 118 , 11225–11236 (1996). https://doi.org/10.1021/ja9621760 Kärger, J., Ruthven, D.M., Theodorou, D.N.: Diffusion in nanoporous materials. Wiley-VCH, Weinheim, Germany (2012) Breck, D.W.: Zeolite Molecular Sieves, pp. 671–688. John Willey & Sons, New York (1974) Glasstone, S., Laidler, K.J., Eyring, H.: The Theory of Rate Processes: The Kinetics of Chemical Reactions, Viscosity, Diffusion and Electrochemical Phenomena. McGraw-Hill, New York, USA (1941) Webster, C.E., Drago, R.S., Zerner, M.C.: Molecular dimensions for adsorptives. J. Am. Chem. Soc. 120 , 5509–5516 (1998). https://doi.org/10.1021/ja973906m Fujii, A., Morita, S., Miyazaki, M., Ebata, T., Mikami, N.: A Molecular Cluster Study on Activated CH/π Interactions: Infrared Spectroscopy of Aromatic Molecule – Acetylene Clusters. J. Phys. Chem. A. 108 (14), 2652–2658 (2004). https://doi.org/10.1021/jp049946b Tsuzuki, S., Fujii, A.: Nature and physical origin of CH/π interaction: significant difference from conventional hydrogen bonds. Phys. Chem. Chem. Phys. 10 , 2584–2594 (2008). https://doi.org/10.1039/B718656H Pople, J.A.: Proton Magnetic Resonance of Hydrocarbons. J. Chem. Phys. 24 , 1111 (1956). https://doi.org/10.1063/1.1742701 Sahakyan, A.B., Vendruscolo, M.: Analysis of the Contributions of Ring Current and Electric Field Effects to the Chemical Shifts of RNA Bases. J. Phys. Chem. B. 117 , 1989–1998 (2013). https://doi.org/10.1021/jp3057306 Scheiner, S.: Dependence of NMR chemical shifts upon CH bond lengths of a methyl group involved in a tetrel bond. Chem. Phys. Lett. 714 , 61–64 (2019). https://doi.org/10.1016/j.cplett.2018.10.069 Additional Declarations No competing interests reported. Supplementary Files SupportinginformationTUsubmitted.docx GraphicalAbstractTU.tif CHEX6MR.mp4 13CHD6MR.mp4 CHEN6MR.mp4 14CHD6MR.mp4 BZ6MR.mp4 Cite Share Download PDF Status: Published Journal Publication published 19 Mar, 2026 Read the published version in Adsorption → Version 1 posted Editorial decision: Revision requested 14 Sep, 2025 Reviews received at journal 04 Sep, 2025 Reviews received at journal 01 Sep, 2025 Reviews received at journal 27 Aug, 2025 Reviewers agreed at journal 20 Aug, 2025 Reviewers agreed at journal 19 Aug, 2025 Reviewers agreed at journal 19 Aug, 2025 Reviewers invited by journal 19 Aug, 2025 Editor assigned by journal 04 Aug, 2025 Submission checks completed at journal 04 Aug, 2025 First submitted to journal 02 Aug, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-7276006","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":503826320,"identity":"c7df745c-b6b4-468f-a1ff-26ed4999e9fd","order_by":0,"name":"Takahiro Ueda","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABG0lEQVRIiWNgGAWjYHACA4YPDAcgzAQQwQOXOYBTC+MMkrUw86BI8uBSCFN/vHnjZ5uaO3K6DewPPzxgsLHn7zn8AGiInTwD41ms1hicOVYsnXPsmbHZAR5jiQSGtMQZZ9tA9iYbNjCcS8Cq5UaOgXRuw+HEbQd42IB+OZzAcJ7B/DcPAzNQ+RkDrFruvzH+bdlwuH7bAfZnIC328ufZPwBtqcet5QaPmTRjw+EEswMMZiAtjBvO9oAcdhinFskzaWWWPccOG247DPKLQVrixjNnChjnGBw3bMPhF77jhzff+FFzWN7sePvDjz8qbOzlzqRvYHhTUS3PL4E9xBCAGexOuIMZGNgkzuDXgQXw95CsZRSMglEwCoYlAADzsF/d0bxu3AAAAABJRU5ErkJggg==","orcid":"","institution":"The University of Osaka","correspondingAuthor":true,"prefix":"","firstName":"Takahiro","middleName":"","lastName":"Ueda","suffix":""},{"id":503826321,"identity":"6204bace-9baa-4973-8580-d9cf59c604f3","order_by":1,"name":"Yuta Yamada","email":"","orcid":"","institution":"The University of Osaka","correspondingAuthor":false,"prefix":"","firstName":"Yuta","middleName":"","lastName":"Yamada","suffix":""},{"id":503826322,"identity":"e7af3413-c4b5-41fc-a905-c07b7893cfc7","order_by":2,"name":"Kota Fujii","email":"","orcid":"","institution":"The University of Osaka","correspondingAuthor":false,"prefix":"","firstName":"Kota","middleName":"","lastName":"Fujii","suffix":""},{"id":503826323,"identity":"b611cd60-81a0-42c6-89cd-0611d884df9d","order_by":3,"name":"Ryota Mihara","email":"","orcid":"","institution":"Shinshu University","correspondingAuthor":false,"prefix":"","firstName":"Ryota","middleName":"","lastName":"Mihara","suffix":""},{"id":503826324,"identity":"76ced65e-08c0-418f-a24f-7c0caf88c5ac","order_by":4,"name":"Taku Iiyama","email":"","orcid":"","institution":"Shinshu University","correspondingAuthor":false,"prefix":"","firstName":"Taku","middleName":"","lastName":"Iiyama","suffix":""},{"id":503826325,"identity":"4915f6e8-510b-4b62-8d16-37ce2444ccd3","order_by":5,"name":"Yasutaka Hamada","email":"","orcid":"","institution":"The University of Osaka","correspondingAuthor":false,"prefix":"","firstName":"Yasutaka","middleName":"","lastName":"Hamada","suffix":""},{"id":503826326,"identity":"c1ad66a4-3652-4c69-9fd4-ddd964410a0d","order_by":6,"name":"Mitsutaka Okumura","email":"","orcid":"","institution":"The University of Osaka","correspondingAuthor":false,"prefix":"","firstName":"Mitsutaka","middleName":"","lastName":"Okumura","suffix":""}],"badges":[],"createdAt":"2025-08-02 06:38:14","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-7276006/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-7276006/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s10450-026-00675-2","type":"published","date":"2026-03-19T15:58:35+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":90037511,"identity":"8e629cd2-b52c-4168-a6b7-43fe253990b1","added_by":"auto","created_at":"2025-08-27 16:08:02","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1296947,"visible":true,"origin":"","legend":"\u003cp\u003eTime dependence of the relative mass of benzene and six-membered ring alicyclic hydrocarbons adsorbed on ZIF-8 at 303 K. The solid lines represent the results of fitting Equation (1) to the linear section of each dataset.\u003c/p\u003e","description":"","filename":"image1.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/8be3c2da42e78c213fa9b8cf.png"},{"id":90037512,"identity":"9c80b86a-ed7c-461f-862a-a86e0a716bf2","added_by":"auto","created_at":"2025-08-27 16:08:03","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1004990,"visible":true,"origin":"","legend":"\u003cp\u003eEyring plot of the diffusion coefficients of all hydrocarbon molecules\u003c/p\u003e","description":"","filename":"image2.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/5b4a8b159c04a96f6059a27b.png"},{"id":90038177,"identity":"0ce9fcd4-c35d-4600-be8d-7eed8835d768","added_by":"auto","created_at":"2025-08-27 16:16:03","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":938540,"visible":true,"origin":"","legend":"\u003cp\u003eMolecular size dependence of Fickian diffusion coefficients, \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e, for adsorbate molecules at 303 K and the activation parameters, \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e (open symbols) and ‒Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e‡\u003c/sup\u003e (filled symbols)\u003c/p\u003e","description":"","filename":"image3.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/3bc3390100a824d5c15b6317.png"},{"id":90037515,"identity":"6149a16a-fbc4-4509-8c6d-b05c194a7e3e","added_by":"auto","created_at":"2025-08-27 16:08:03","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":844835,"visible":true,"origin":"","legend":"\u003cp\u003ePotential energy profiles obtained through Monte Carlo simulation for the adsorbed molecules as they pass through the 6MR aperture\u003c/p\u003e","description":"","filename":"image4.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/3c2654ca8afa563b1f5ee774.png"},{"id":90037524,"identity":"fea45dab-6e7e-47b0-838a-3cc77a963d3d","added_by":"auto","created_at":"2025-08-27 16:08:03","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":2146242,"visible":true,"origin":"","legend":"\u003cp\u003eOrder parameters for adsorbate molecules while passing through the 6MR aperture\u003c/p\u003e","description":"","filename":"image5.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/df2cdf7e4701ad0aa56d2a47.png"},{"id":90038184,"identity":"85441ab8-6c4b-424e-a1a8-ce66add52fea","added_by":"auto","created_at":"2025-08-27 16:16:03","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":510418,"visible":true,"origin":"","legend":"\u003cp\u003eFT-IR spectra of the symmetric C–H stretching band of the 2-methylimidazolate group in ZIF-8 without guests and after saturation with each hydrocarbon\u003c/p\u003e","description":"","filename":"image6.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/acfe8eae4783c1f9388d2718.png"},{"id":90037532,"identity":"115ecd62-dc83-4353-acc2-ab21acdd60e2","added_by":"auto","created_at":"2025-08-27 16:08:03","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":535842,"visible":true,"origin":"","legend":"\u003cp\u003e\u003csup\u003e13\u003c/sup\u003eC CP/MAS NMR spectra of ZIF-8 and ZIF-8 with various adsorbate molecules. The asterisk indicates the spinning sideband. The minor peaks correspond to the adsorbate molecules\u003c/p\u003e","description":"","filename":"image7.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/930267d08d36ab632d88b0a8.png"},{"id":90037519,"identity":"ac54d053-68fb-44c1-abdb-6947ae4c18f8","added_by":"auto","created_at":"2025-08-27 16:08:03","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":536366,"visible":true,"origin":"","legend":"\u003cp\u003eExtended view of \u003csup\u003e13\u003c/sup\u003eC CPMAS NMR spectra in the chemical shift region of the imidazolate quaternary carbon (a), imidazolate ring CH carbon (b), and methyl carbon (c). The adsorption of benzene and cyclic alkenes induced low-field shifts of 0.5‒1.3 ppm in the resonance lines of the methyl group and CH carbon of the imidazole ring\u003c/p\u003e","description":"","filename":"image8.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/ebd09279bc24343ae788e07c.png"},{"id":90038181,"identity":"a20bd43e-0458-40ae-bb4c-05b1329b4cbd","added_by":"auto","created_at":"2025-08-27 16:16:03","extension":"png","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":1265337,"visible":true,"origin":"","legend":"\u003cp\u003eModel of the configuration of benzene and six-membered-ring alicyclic hydrocarbons traversing the 6MR aperture in ZIF-8\u003c/p\u003e","description":"","filename":"image9.png","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/39de55c0b5770d6c31f778df.png"},{"id":105223411,"identity":"a9dc59cd-9ffa-4b27-ba50-c60bef97d069","added_by":"auto","created_at":"2026-03-23 16:05:51","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":10184314,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/0862a838-c6e6-4cd9-b8cd-37d5342b1b7b.pdf"},{"id":90037523,"identity":"adbf0a6d-a755-41b6-9d48-91ec7aa9beb5","added_by":"auto","created_at":"2025-08-27 16:08:03","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":2503358,"visible":true,"origin":"","legend":"","description":"","filename":"SupportinginformationTUsubmitted.docx","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/1f23bc1a5b735721680e2a99.docx"},{"id":90039069,"identity":"0c7d7b20-ea5d-4acd-a5d8-dea5a1fa7dba","added_by":"auto","created_at":"2025-08-27 16:24:04","extension":"tif","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":24224030,"visible":true,"origin":"","legend":"","description":"","filename":"GraphicalAbstractTU.tif","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/020052c44d61f970bd2cabf9.tif"},{"id":90038178,"identity":"d4db653f-970f-4aec-a9cd-71be3f449e44","added_by":"auto","created_at":"2025-08-27 16:16:03","extension":"mp4","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":4146719,"visible":true,"origin":"","legend":"","description":"","filename":"CHEX6MR.mp4","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/46afea8631e78ab2e05d7295.mp4"},{"id":90038179,"identity":"8cd76917-d5a0-4fab-8bf7-b39cf0655eff","added_by":"auto","created_at":"2025-08-27 16:16:03","extension":"mp4","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":4106108,"visible":true,"origin":"","legend":"","description":"","filename":"13CHD6MR.mp4","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/dfebf8c0432c55ed94176582.mp4"},{"id":90037537,"identity":"7e5f716d-ecc3-40c2-bc8c-6690c15d5d47","added_by":"auto","created_at":"2025-08-27 16:08:03","extension":"mp4","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":4017145,"visible":true,"origin":"","legend":"","description":"","filename":"CHEN6MR.mp4","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/bf0661ea99b53b93a1be1936.mp4"},{"id":90038180,"identity":"f03a818f-d025-4985-9fda-3be61f336a58","added_by":"auto","created_at":"2025-08-27 16:16:03","extension":"mp4","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":4080833,"visible":true,"origin":"","legend":"","description":"","filename":"14CHD6MR.mp4","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/3de817d3ff3a2c7170ab78af.mp4"},{"id":90038188,"identity":"41629202-0962-41db-a085-c8facf0fdabb","added_by":"auto","created_at":"2025-08-27 16:16:03","extension":"mp4","order_by":7,"title":"","display":"","copyAsset":false,"role":"supplement","size":3869392,"visible":true,"origin":"","legend":"","description":"","filename":"BZ6MR.mp4","url":"https://assets-eu.researchsquare.com/files/rs-7276006/v1/484d9d21c7bf8b3add02d09d.mp4"}],"financialInterests":"No competing interests reported.","formattedTitle":"Intermolecular interactions governing adsorption of bulky molecules on ZIF-8: Insight from adsorption kinetics of benzene and 6-membered ring alicyclic hydrocarbons","fulltext":[{"header":"1 Introduction","content":"\u003cp\u003eThe growing interest in hydrogen as a prospective energy source has led to increased demand for hydrogen as a renewable resource in recent years. To establish a dependable hydrogen supply, there is an urgent need to construct large-scale storage facilities, and economically viable transportation systems [\u003cspan additionalcitationids=\"CR2\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. The organic hydride process is a promising strategy for hydrogen supply [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e] involving the storage of hydrogen in saturated hydrocarbons such as methylcyclohexane, which can be dehydrogenated to produce unsaturated hydrocarbons like toluene for hydrogen extraction [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e, \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. The efficiency of the hydrogen extraction process can be improved by continuously removing toluene from the reaction system. Therefore, the development of molecular sieve technologies capable of effectively separating toluene from methylcyclohexane\u0026ndash;toluene mixtures while minimizing energy consumption is essential. Beyond the methylcyclohexane/toluene system, the advancement of molecular sieve materials for the selective separation and adsorption of saturated and unsaturated compounds in bulky hydrocarbons with high hydrogen densities is anticipated to further enhance the organic hydride method [\u003cspan additionalcitationids=\"CR9 CR10 CR11\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e].\u003c/p\u003e\u003cp\u003eZeolitic Imidazolate Framework 8 (ZIF-8), a type of metal-organic framework, is one of the leading candidates for such molecular sieve materials. The ZIF-8 is represented by the formula [Zn(C\u003csub\u003e4\u003c/sub\u003eH\u003csub\u003e5\u003c/sub\u003eN\u003csub\u003e2\u003c/sub\u003e)\u003csub\u003e2\u003c/sub\u003e]\u003csub\u003e\u003cem\u003en\u003c/em\u003e\u003c/sub\u003e, consisting of Zn\u003csup\u003e2+\u003c/sup\u003e ions and the 2-methylimidazolate anion (C\u003csub\u003e4\u003c/sub\u003eH\u003csub\u003e5\u003c/sub\u003eN\u003csub\u003e2\u003c/sub\u003e\u003csup\u003e\u0026minus;\u003c/sup\u003e, MeIm\u003csup\u003e\u0026minus;\u003c/sup\u003e) [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. It possesses a sodalite structure characterized by micropores with a diameter of 11.4 \u0026Aring;, interconnected by eight 6-membered ring (6MR) openings measuring 3.4 \u0026Aring; in diameter [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. This material exhibits significant potential for gas adsorption and separation applications [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e], which can be attributed to its exceptional molecular sieving properties stemming from its flexible framework. The rotational motion of the bridging ligand modifies the apparent size of the 6MR opening, allowing the capture of molecules larger than 3.4 \u0026Aring; [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. This phenomenon, referred to as \"gate adsorption,\" results from the framework's flexibility [\u003cspan additionalcitationids=\"CR20\" citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. The unique molecular sieving behavior of ZIF-8 is characterized by a nearly linear relationship between the activation energy for gas molecular diffusion and the size of the adsorbed molecules [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e], although the effective diffusion coefficient derived from the adsorption rate decreases with the kinetic diameter of the adsorbed molecules. This diffusion behavior is linked to the swing effect, which is induced by the torsional motion of the imidazole ligand around the Zn\u0026ndash;Im\u0026ndash;Zn axis [\u003cspan additionalcitationids=\"CR23\" citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. Significant fluctuations in the bridging ligands have been confirmed using various spectroscopic techniques, including far-infrared absorption spectroscopy [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e], inelastic neutron scattering [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e], terahertz spectroscopy [\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e], and solid-state NMR [\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e] methods.\u003c/p\u003e\u003cp\u003eHowever, bulky molecules that exhibit significant molecular shape anisotropy may adopt an energetically favorable orientation for navigating through the 6MR aperture in an orientation that minimizes the molecular cross-sectional area. In previous studies, it was observed that for benzene and its derivatives (toluene and xylene), the minimum molecular cross-section rather than the average molecular diameter predominantly influenced the molecular orientation of the adsorbed molecules as they traversed the 6MR aperture [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e, \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. Recently, Poryvaev et al. reported that benzene was adsorbed onto ZIF-8 with an adsorption selectivity of approximately 270% from a benzene/cyclohexane solution containing\u0026thinsp;\u0026gt;\u0026thinsp;90% cyclohexane [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]. Their findings indicated that the interaction between benzene and cyclohexane with the 6MR aperture of ZIF-8 was governed by the molecular thickness of the adsorbed molecules, suggesting that the molecular orientation was critical for the passage of benzene and cyclohexane through the 6MR aperture. Thus, for adsorbed molecules with bulky and highly anisotropic geometries, both the intermolecular interactions with the 6MR aperture (enthalpic contributions) and the molecular orientation during passage through the aperture (entropic contributions) were expected to significantly influence their diffusion behavior.\u003c/p\u003e\u003cp\u003eTo sustain a favorable molecular orientation for passage through the 6MR aperture, highly directional interactions likely govern the orientation of the adsorbed molecules, as opposed to isotropic intermolecular interactions such as dispersion forces. Identifying the intermolecular interactions that contribute to the maintenance of molecular orientation is anticipated to enhance the theoretical understanding of the adsorption behavior of bulky molecules on ZIF-8, particularly regarding diffusion phenomena. Potential candidates for such interactions include hydrogen bonds, tetrel bonds, CH-π interactions, and halogen-π interactions [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. We hypothesized that the double bond of 2-methylimidazolate, the bridging ligand of ZIF-8, would facilitate CH-π interactions.\u003c/p\u003e\u003cp\u003eThis study aimed to elucidate the molecular-level adsorption mechanism of bulky hydrocarbons in ZIF-8. Initially, we focused on the adsorption kinetics of benzene and several 6-membered ring alicyclic hydrocarbons on ZIF-8, with an emphasis on evaluating the influence of intermolecular interactions on diffusion behavior. The selected alicyclic compounds (1,3-cyclohexadiene, 1,4-cyclohexadiene, cyclohexene, and cyclohexane, summarized in Scheme S1 in supporting information) had a similar projected molecular diameter to benzene but differed in molecular thickness and π-electron content (see Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e, in SI). These variations facilitate the analysis of the diffusion coefficients and activation parameters in relation to the molecular dimensions and the presence of π-electrons. This study employed Monte Carlo simulations to investigate the critical process of molecules traversing the 6MR aperture of ZIF-8, with particular emphasis on the orientation of bulky molecules during this process. Additionally, Fourier-transform infrared absorption spectroscopy and \u003csup\u003e13\u003c/sup\u003eC cross-polarization magic-angle sample spinning (CPMAS) NMR techniques were utilized to assess adsorbate-adsorbent interactions, aiming to identify specific intermolecular interactions that influence molecular orientation as the adsorbate molecules navigate through the 6MR aperture. These research strategies collectively offer molecular-level understanding of the adsorption mechanism of bulky hydrocarbon molecules on ZIF-8.\u003c/p\u003e"},{"header":"2 Experimental","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\u003ch2\u003e2.1 Chemicals\u003c/h2\u003e\u003cp\u003eThe ZIF-8 (Basolite Z1200, BASF) and 1,3-cyclohexadiene (97% purity) were obtained from Sigma-Aldrich. Benzene (purity\u0026thinsp;\u0026gt;\u0026thinsp;99.5%), cyclohexane (purity\u0026thinsp;\u0026gt;\u0026thinsp;99.5%), 1,4-cyclohexadiene (purity\u0026thinsp;\u0026gt;\u0026thinsp;95.0%), and cyclohexene (purity\u0026thinsp;\u0026gt;\u0026thinsp;97.0%) were purchased from Fujifilm Wako Pure Chemicals Co., Ltd. The reagents utilized were Wako Special Grade for benzene and cyclohexane and Wako 1st Grade for 1,4-cyclohexadiene and cyclohexene. All the reagents were used without further purification. For this study, benzene, 1,3-cyclohexadiene, 1,4-cyclohexadiene, cyclohexene, and cyclohexane are referred to as BZ, 1,3-CHD, 1,4-CHD, CHEN, and CHEX, respectively.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec4\" class=\"Section2\"\u003e\u003ch2\u003e2.2 Characterization\u003c/h2\u003e\u003cdiv id=\"Sec5\" class=\"Section3\"\u003e\u003ch2\u003e2.2.1 Nitrogen adsorption isotherm\u003c/h2\u003e\u003cp\u003ePrior to use, ZIF-8 was activated by heating under vacuum (0.1 Torr) at 150\u0026deg;C for 12 h. Nitrogen adsorption isotherms were recorded using a Gemini 2375 apparatus (Shimadzu Corporation) to determine the pore volume and specific surface area of ZIF-8. The isotherms were recorded over a relative pressure (\u003cem\u003ep\u003c/em\u003e/\u003cem\u003ep\u003c/em\u003e\u003csub\u003e0\u003c/sub\u003e) range of 0.0001\u0026ndash;0.98, where \u003cem\u003ep\u003c/em\u003e\u003csub\u003e0\u003c/sub\u003e represents the saturated vapor pressure of nitrogen at its boiling point of liquid nitrogen (77 K). The porosity of ZIF-8 was evaluated from the nitrogen adsorption isotherm (Figure \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e, Supplementary Information), which exhibited a characteristic Type I (a) pattern according to the IUPAC classification. The BET surface area was calculated to be 1074 m\u003csup\u003e2\u003c/sup\u003e g\u003csup\u003e\u0026ndash;1\u003c/sup\u003e, with a corresponding pore volume of 0.61 cm\u003csup\u003e3\u003c/sup\u003e g\u003csup\u003e\u0026ndash;1\u003c/sup\u003e.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec6\" class=\"Section3\"\u003e\u003ch2\u003e2.2.2 Scanning electron microscopy (SEM)\u003c/h2\u003e\u003cp\u003eScanning electron microscopy (SEM) was performed using a field-emission scanning electron microscope SU6600 (Hitachi High-Technologies Co.) to analyze the size distribution and morphology of the crystallites. The SEM images are shown in Figure \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e (Supplementary Information). The average particle diameter, based on the crystalline size distribution (Figure \u003cspan refid=\"MOESM3\" class=\"InternalRef\"\u003eS3\u003c/span\u003e), was determined to be 450 nm, corresponding to a radius of 225 nm.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec7\" class=\"Section3\"\u003e\u003ch2\u003e2.2.3 Fourier transferred infrared (FT-IR) absorption spectroscopy\u003c/h2\u003e\u003cp\u003eFourier-transform infrared (FT-IR) measurements were performed using an FT/IR-4600AC spectrometer (JASCO Co.) over the wavenumber range of 400\u0026ndash;4000 cm\u003csup\u003e\u0026ndash;1\u003c/sup\u003e at room temperature with a resolution of 1 cm\u003csup\u003e\u0026ndash;1\u003c/sup\u003e. To prevent the amorphization of the ZIF-8 sample due to uneven pressure during mortar grinding, KBr pellet samples were prepared by compressing the powder sample between two KBr plates. Spectra were recorded in transmittance units and subsequently converted to absorbance values.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec8\" class=\"Section3\"\u003e\u003ch2\u003e2.2.4 Solid-state \u003csup\u003e13\u003c/sup\u003eC NMR spectroscopy\u003c/h2\u003e\u003cp\u003eHigh-resolution solid-state \u003csup\u003e13\u003c/sup\u003eC NMR spectra were obtained using a Bruker Avance III 400WB spectrometer equipped with a superconducting magnet (9.4 T), which provided resonance frequencies of 400.25 and 100.64 MHz for the \u003csup\u003e1\u003c/sup\u003eH and \u003csup\u003e13\u003c/sup\u003eC nuclei, respectively. Free induction decay (FID) signals for \u003csup\u003e13\u003c/sup\u003eC nuclei were acquired using a single-pulse sequence with \u003csup\u003e1\u003c/sup\u003eH decoupling under magic-angle sample spinning (MAS) conditions and a \u003csup\u003e1\u003c/sup\u003eH decoupling MAS cross-polarization (CP) pulse sequence. The powdered samples were packed in a 4 mm-diameter cylindrical zirconia rotor and spun at 10 kHz for MAS. Adamantane was used as an external standard for the chemical shift, with a CH\u003csub\u003e2\u003c/sub\u003e peak referenced at 38.5 ppm (relative to tetramethylsilane).\u003c/p\u003e\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Sec9\" class=\"Section2\"\u003e\u003ch2\u003e2.3 Measurement of time dependence of adsorption uptake\u003c/h2\u003e\u003cp\u003eThe ZIF-8 powder sample underwent pretreatment (heating at 150\u0026deg;C in a vacuum of 0.1 Torr for 12 h). The samples with an average weight of 7\u0026ndash;15 mg were placed in an aluminum pan (5 mm \u0026times; 3 mm depth). The aluminum pan was placed in a hermetically sealed vessel containing 3 mL of the appropriate solvent. The presence of air in the container had a minimal effect on solvent vapor adsorption, as nitrogen (N\u003csub\u003e2\u003c/sub\u003e) and oxygen (O\u003csub\u003e2\u003c/sub\u003e) behave as supercritical gases at ambient temperatures. The adsorbates used in this study were BZ, 1,3-CHD, 1,4-CHD, CHEN, and CHEX. The vessel was affixed to a base constructed from thermally insulating styrofoam and placed within a temperature-controlled enclosure with an accuracy of \u0026plusmn;\u0026thinsp;0.5\u0026deg;C. Vapor adsorption was conducted at 30, 40, and 50\u0026deg;C under saturated vapor pressures of the solvent. The temperature dependence of the saturated vapor pressure is shown in Figure \u003cspan refid=\"MOESM4\" class=\"InternalRef\"\u003eS4\u003c/span\u003e (Supplementary Information). At a constant temperature, BZ, 1,3-CHD, and CHEX exhibited similar vapor pressures within \u0026plusmn;\u0026thinsp;1 kPa, whereas CHEN and 1,4-CHD displayed slightly lower vapor pressures, with differences of 1.5 kPa for CHEN and 5 kPa for 1,4-CHD. Consequently, the influence of the vapor pressure of each solvent on the adsorption rate was expected to be negligible. The samples were weighed at intervals until a stable value was reached, with weight measurements exhibiting an accuracy of \u0026plusmn;\u0026thinsp;5%.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e\u003ch2\u003e2.4 Monte Carlo simulations\u003c/h2\u003e\u003cp\u003eA Monte Carlo simulation was performed using proprietary software at a constant temperature of 298 K under NVT conditions, with a minimum of 100,000 trials. The force field parameters \u0026mdash;including Lennard-Jones parameters and partial charges\u0026mdash;assigned to the atoms forming the 6MR aperture were consistent with those used in the molecular dynamics (MD) simulations of ZIF-8 conducted by Parkes et al. [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]. The configuration of the adsorbed molecule was determined by optimizing its structure using density functional theory (DFT) calculations. The Gaussian 16 software package [\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e] was employed, utilizing the B3LYP functional and 6-31G(d) basis set for the calculations. In the Monte Carlo simulations, the OPLS-AA force field parameters [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e, \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e] were used to define the Lennard-Jones parameters and partial charges of the adsorbed molecules. The 6MR structure employed in the simulation was modeled by extracting one-eighth of the unit cells dimension (\u003cem\u003ea\u003c/em\u003e/2, \u003cem\u003ea\u003c/em\u003e/2, \u003cem\u003ea/\u003c/em\u003e2) from the ZIF-8 crystal (lattice constant \u003cem\u003ea\u003c/em\u003e). The vector from the lattice point (0, 0, 0) to (1/2, 1/2, 1/2) was designated as the \u003cem\u003eZ\u003c/em\u003e-axis in the reference coordinate system. The origin of this reference coordinate system was set at the center of the 6MR plane (lattice point (1/4, 1/4, 1/4)). The direction toward lattice point (0, 0, 0) was defined as negative (\u0026ndash;), and that toward (1/2, 1/2, 1/2) as positive (+). The 6MR structure is shown in Figure \u003cspan refid=\"MOESM5\" class=\"InternalRef\"\u003eS5\u003c/span\u003e. At a specified \u003cem\u003eZ\u003c/em\u003e value, the mass center of the adsorbed molecule was displaced over a 4 \u0026Aring; \u0026times; 4 \u0026Aring; plane centered on the reference coordinate axis. The most energetically favorable configuration was then determined using Monte Carlo simulation. The opening of the 6MR aperture was adjusted by systematically rotating the tilt angle of the 2-methylimidazole bridging ligand along its N-N axis. The tilt angle was defined as the angle between the molecular plane of the bridging ligand and the 6MR aperture plane. Variations in tilt were expressed as positive (+Δ\u003cem\u003eθ\u003c/em\u003e) or negative (\u0026ndash;Δ\u003cem\u003eθ\u003c/em\u003e) deviations from the ligand\u0026rsquo;s original orientation of the crystal lattice. The adsorbed molecule was moved along the Z-axis in the reference coordinate system from +\u0026thinsp;4 \u0026Aring; to \u0026minus;\u0026thinsp;4 \u0026Aring; in 0.2 \u0026Aring; increments, and the potential energy between the adsorbed molecule and the 6MR for each \u003cem\u003eZ\u003c/em\u003e value was evaluated by varying Δ\u003cem\u003eθ\u003c/em\u003e. For the order parameter analysis, the simulations were carried out in the \u003cem\u003eZ\u003c/em\u003e value from +\u0026thinsp;5 \u0026Aring; to \u0026minus;\u0026thinsp;5 \u0026Aring;.\u003c/p\u003e\u003c/div\u003e"},{"header":"3 Results and Discussion","content":"\u003cdiv id=\"Sec12\" class=\"Section2\"\u003e\u003ch2\u003e3.1 Time dependence of adsorption uptake\u003c/h2\u003e\u003cp\u003eFigure \u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e depicts the time-dependence of the adsorption uptake of individual hydrocarbons on ZIF-8 at 303 K, illustrating the relationship between the amount of adsorption and the square root of the adsorption time. Initially, all the vapor species exhibited a linear increase in the adsorption amount, which was proportional to the square root of time, suggesting that micropore diffusion was the primary mechanism governing the adsorption of individual solvent vapors on ZIF-8. These observations can be characterized using the Fickian diffusion model, which is applicable to spherical particles [\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e]. Assuming that the concentration of adsorbed molecules at the center of the particle is negligible (i.e., in the short diffusion time region), the model can be approximated as\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\n$$\\:m\\left(t\\right)={m}_{\\text{e}\\text{x}\\text{t}}+{m}_{\\infty\\:}\\frac{6}{\\sqrt{\\pi\\:}}\\sqrt{\\frac{{D}_{\\text{f},\\text{i}\\text{n}\\text{i}}t}{{r}_{\\text{a}\\text{v}\\text{e}}^{2}}}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e1\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003ewhere \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e denotes the Fickian diffusion coefficient, which is derived from the initial slope of the time-dependent adsorption curve; \u003cem\u003er\u003c/em\u003e\u003csub\u003eave\u003c/sub\u003e represents the average radius of the spherical particles; \u003cem\u003em\u003c/em\u003e\u003csub\u003e\u0026infin;\u003c/sub\u003e indicates the saturation adsorption amount; and \u003cem\u003em\u003c/em\u003e\u003csub\u003eext\u003c/sub\u003e refers to the adsorption amount present on the surface of the particle [\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. The diffusion coefficient \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e can be calculated from the slope of the linear portion depicted in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, assuming that the saturation adsorption amount is known. The data indicated that BZ and 1,4-CHD achieved saturation during adsorption, whereas 1,3-CHD approached saturation. In contrast, the CHEN and CHEX did not reach saturation during this period. For convenience, the adsorption quantity recorded at the longest adsorption time was designated as the saturation adsorption quantity, \u003cem\u003em\u003c/em\u003e\u003csub\u003e\u0026infin;\u003c/sub\u003e, at equilibrium. The \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e values for each adsorbed molecule were extracted from the linear segments, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. Figure \u003cspan refid=\"MOESM6\" class=\"InternalRef\"\u003eS6\u003c/span\u003e (SI) presents the time-dependent adsorption amounts at 313 and 323 K. Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e summarizes the \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e values obtained at 303, 313, and 323 K. Notably, despite being structural isomers with identical molecular weights, 1,3-CHD and 1,4-CHD exhibited distinct diffusion coefficients. This finding suggests that ZIF-8 can differentiate between structural isomers based on their respective adsorption behavior through kinetic mechanisms. This observation is particularly significant, as it implies that ZIF-8 can potentially serve as a molecular sieve for structural isomers\u0026mdash;an ability with important implications for separation applications.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003eDiffusion coefficients of adsorbate molecule determined from the initial slope of adsorption uptake curve\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"4\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eAdsorbate\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colspan=\"3\" nameend=\"c4\" namest=\"c2\"\u003e\u003cp\u003e\u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e / m\u003csup\u003e2\u003c/sup\u003e s\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e303 K\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e313 K\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e323 K\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eBZ\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e1.3 (1) \u0026times; 10\u003csup\u003e‒19\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e2.2 (1) \u0026times; 10\u003csup\u003e‒19\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e2.8 (1) \u0026times; 10\u003csup\u003e‒19\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e1,4-CHD\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e5.6 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e7.9 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e9.5 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e1,3-CHD\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e3.2 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e3.5 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e5.4 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCHEN\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e1.2 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e1.8 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e2.2 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCHEX (chair)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e8.2 (1) \u0026times; 10\u003csup\u003e‒21\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e1.1 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e1.7 (1) \u0026times; 10\u003csup\u003e‒20\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003eThe member in parentheses presents the experimental error.\u003c/p\u003e\u003cp\u003eWhen larger molecules are adsorbed onto ZIF-8, they must navigate through the 6MR aperture, which consists of six zinc (Zn) atoms and 2-methylimidazole bridging ligands. A previous study established that, for benzene and its derivatives adsorbed on ZIF-8 [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e], traversal through the 6MR aperture represents the rate-limiting step in the diffusion process. Furthermore, the movement of bulky molecules through the 6MR aperture is primarily governed by configurational diffusion, which is characterized by a restricted molecular orientation and rotation. Consequently, the temperature dependence of \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e should be assessed using the following equation derived from Eyring's transition state theory [\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e] rather than the traditional Arrhenius equation:\u003cdiv id=\"Equ2\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ2\" name=\"EquationSource\"\u003e\n$$\\:{D}_{\\text{f},\\text{i}\\text{n}\\text{i}}=\\text{e}{\\lambda\\:}^{2}\\frac{kT}{h}{\\text{e}}^{\\varDelta\\:{S}^{\\ddagger}/R}{\\text{e}}^{{-E}_{a}/RT}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e2\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003ewhere \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e signifies the diffusion activation energy, Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e represents the diffusion activation entropy, \u003cem\u003eh\u003c/em\u003e denotes Planck's constant, \u003cem\u003ek\u003c/em\u003e is Boltzmann's constant, and λ indicates the distance between the molecular jumps. During the diffusion process, molecules jump between adjacent lattices via 6MR apertures with λ\u0026thinsp;=\u0026thinsp;3\u003csup\u003e1/2\u003c/sup\u003e\u003cem\u003ea\u003c/em\u003e/2. The lattice constant \u003cem\u003ea\u003c/em\u003e is 1.701 nm at 298 K [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e], resulting in a λ of 1.473 nm.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eFigure \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e illustrates a logarithmic representation of \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003eꞏ\u003cem\u003eT\u003c/em\u003e \u003csup\u003e‒1\u003c/sup\u003e versus the reciprocal temperature (Eyring plot). The slope and intercept of the graph correspond to the activation energy and entropy, respectively. The activation parameters for each adsorbed molecule are listed in Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e.\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003eDiffusion coefficients at 303 K and activation parameters determined from the temperature dependence of the diffusion coefficients.\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"4\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eAdsorbate\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003e\u003cem\u003eσ\u003c/em\u003e\u003csub\u003eaverage\u003c/sub\u003e /nm\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003e\u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e / kJ mol\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e\u003cp\u003eΔ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e/ J K\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e mol\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eBZ\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e0.525\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e28 (2)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e‒172 (5)\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e1,4-CHD\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e0.537\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e19 (2)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e‒208 (10)\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e1,3-CHD\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e0.538\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e18 (2)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e‒217 (10)\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCHEN\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e0.548\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e24 (2)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e‒205 (6)\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCHEX (chair)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e\u003cp\u003e0.558\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e27 (1)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e‒200 (3)\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003ctfoot\u003e\u003ctr\u003e\u003ctd colspan=\"4\"\u003e‒\u003c/td\u003e\u003c/tr\u003e\u003c/tfoot\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eFigure \u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e presents the molecular size dependence of \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e at 303 K and the activation parameters \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e and ‒Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e for individual hydrocarbons. The semi-logarithmic plot of \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e shows a linear decrease with increasing molecular size, whereas the standard plots of \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e and ‒Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e exhibited a V-shaped dependence on the molecular size. Specifically, \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e reached its minimum for 1,4-CHD and 1,3-CHD, which possessed intermediate molecular sizes, and attained its maximum values for BZ and CHEX, representing the smallest and largest molecular sizes, respectively. This observation suggests that the passage of adsorbed molecules through 6MR aperture is governed by the combination of molecular size and specific host-guest interactions at the 6MR aperture site. The activation energy \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e can be lowered by stabilizing the transition state or raised by stabilizing the adsorbed molecules within the pores. A lower \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e enhances the likelihood of adsorbed molecules overcoming the diffusion energy barrier at a specific temperature. Conversely, Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e exerts an opposing influence on \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e. The substantial negative activation entropy indicates a pronounced restriction on the degrees of freedom of the adsorbed molecule within the diffusion transition state, which leads to a reduced diffusion frequency. In BZ, the effect of Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e\u0026mdash;despite having the smallest absolute value\u0026mdash;surpasses that of the largest \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e, resulting in the highest \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e value among the compounds. For 1,4-CHD and 1,3-CHD, the \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e values are less than 20 kJ mol\u003csup\u003e‒1\u003c/sup\u003e; however, their Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e values reach up to ca. 210 J K\u003csup\u003e‒1\u003c/sup\u003e mol\u003csup\u003e‒1\u003c/sup\u003e. The competing effects of these opposing factors may result in the contribution of activation entropy surpassing that of activation energy, thereby causing a decline in \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e as the molecular size increases. In BZ, 1,4-CHD, and 1,3-CHD, \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e and Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e exhibited a trade-off relationship. In contrast, for CHEN and CHEX, both the high \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e values and the relatively large absolute values of Δ\u003cem\u003eS\u003c/em\u003e\u003csup\u003e\u0026Dagger;\u003c/sup\u003e decrease the \u003cem\u003eD\u003c/em\u003e\u003csub\u003ef,ini\u003c/sub\u003e values. In the following sections, we examine the pronounced effects of molecular orientation on the diffusion of adsorbed molecules through the 6MR and identify the specific intermolecular interactions that significantly influence the \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec13\" class=\"Section2\"\u003e\u003ch2\u003e3.2 Monte Carlo simulation for 6MR passage of adsorbate molecules\u003c/h2\u003e\u003cp\u003eMonte Carlo simulations were employed to elucidate the potential energy profile and molecular orientation of the adsorbed molecules as they traverse the 6MR aperture. In this context, the 6MR aperture, characterized by a tilted angle (Δ\u003cem\u003eθ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0) as derived from the crystal structure, presents an activation energy required for molecular passage with activation energies ranging from several hundred to several thousand kJ mol\u003csup\u003e‒1\u003c/sup\u003e. Consequently, it is nearly impossible for molecules to pass through the aperture using only the thermal energy available at room temperature. This observation implies that the 6MR aperture must undergo enlargement to facilitate the passage of target molecules. To quantify this, we evaluated the activation energy required for molecular passage as a function of the tilted angle of the bridging ligand by systematically varying this angle. By employing polynomial interpolation on the acquired data, we identified a tilted angle that aligned with the experimentally measured values. Figure\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e illustrates the energy profile associated with passage through the 6MR aperture, corresponding to the experimental activation energy for each adsorbate, as determined by the time dependence of the adsorption uptake. Two minima were observed on either side of the 6MR aperture. The energy minimum at the CH\u003csub\u003e3\u003c/sub\u003e neck (\u003cem\u003eZ\u003c/em\u003e/\u0026Aring; ~ \u0026minus;\u0026thinsp;2) was deeper than that at the CH neck (\u003cem\u003eZ\u003c/em\u003e/\u0026Aring; ~ +3). This finding suggests that the passage of molecules trapped in the CH\u003csub\u003e3\u003c/sub\u003e neck through the 6MR aperture is energetically unfavorable. Consequently, the passage of molecules through the 6MR aperture is only considered when they approach from the CH neck side [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. In this case, the energy barrier encountered by an adsorbed molecule navigating through the 6MR aperture is defined as the difference (Δ\u003cem\u003eE\u003c/em\u003e\u0026thinsp;=\u0026thinsp;\u003cem\u003eE\u003c/em\u003e\u003csub\u003eCH,max\u003c/sub\u003e \u0026minus; \u003cem\u003eE\u003c/em\u003e\u003csub\u003eCH,min\u003c/sub\u003e) between the maximum energy (\u003cem\u003eE\u003c/em\u003e\u003csub\u003eCH,max\u003c/sub\u003e) located near the CH neck (+\u0026thinsp;0.5\u0026thinsp;\u0026le;\u0026thinsp;\u003cem\u003eZ\u003c/em\u003e/\u0026Aring; \u0026le; +1) and the minimum energy (\u003cem\u003eE\u003c/em\u003e\u003csub\u003eCH,min\u003c/sub\u003e) at the CH neck (\u003cem\u003eZ\u003c/em\u003e/\u0026Aring; ~ +3). Notably, in the case of CHEX, an additional energy maximum was detected on the outside of the CH neck (\u003cem\u003eZ\u003c/em\u003e/\u0026Aring; ~ +2). The presence of the dual maxima in the energy profile aligns with the findings from previous studies. This behavior supports the use of a model that incorporates pre-equilibrium adsorption at the CH neck, represented as: A(g)\u0026thinsp;+\u0026thinsp;M(s) ⇌ AM*(s) ⇌ AM(s) [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. In this model, the apparent activation energy can be expressed as Δ\u003cem\u003eE\u003c/em\u003e\u0026thinsp;=\u0026thinsp;\u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,a\u003c/sub\u003e \u0026ndash; \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,d\u003c/sub\u003e + \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,in\u003c/sub\u003e, where \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,a\u003c/sub\u003e denotes the energy barrier for pre-adsorption at the 6MR aperture, \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,d\u003c/sub\u003e is the energy barrier for the desorption of the pre-adsorbed molecule, and \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,in\u003c/sub\u003e is the energy barrier for passage through the CH neck. The activation energy evaluated by MC simulation, the resulting tilted angle of the bridging ligand, and the aperture opening diameter are listed in Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e. The correlation between the tilted angle and the opening diameter of the 6MR aperture, which reproduces the experimental values, appears to align more closely with MIN-1 than with MIN-2, as proposed by Webster et al. [\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e] (see Scheme S1 and Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e in the SI). This suggests that the energy barrier for the passage of the adsorbed molecule through the 6MR aperture is significantly influenced by the cross-sectional area of the adsorbed molecule, indicating that the molecular orientation is highly constrained during this process. This conclusion was consistent with the findings of our previous study.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003eExperimental and simulated activation energies of adsorbed molecules passing through the 6MR aperture, tilted angles of the bridging ligand, and average diameters of the 6MR aperture.\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"6\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eAdsorbate\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003eΔ\u003cem\u003eE\u003c/em\u003e / kJ mol\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003eΔ\u003cem\u003eθ\u003c/em\u003e\u003csub\u003ec\u003c/sub\u003e / deg.\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e\u003cp\u003e\u003cem\u003ed\u003c/em\u003e\u003csub\u003ep,ave\u003c/sub\u003e / nm\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c5\"\u003e\u003cp\u003eMIN-1\u003csup\u003eb)\u003c/sup\u003e/nm\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c6\"\u003e\u003cp\u003eMIN-2 \u003csup\u003eb)\u003c/sup\u003e/nm\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eBZ\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e27.6\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e17.0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.47\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.3277\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e0.6628\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e1,4-CHD\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e19.3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e21.1\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.49\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.3845\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e0.6612\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e1,3-CHD\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e17.5\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e19.9\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.49\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.3845\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e0.6612\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCHEN\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e24.2\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e20.55\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.49\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.4414\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e0.6595\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCHEX (chair)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e27.0 \u003csup\u003ea)\u003c/sup\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e24.05\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.52\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.4982\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e\u003cp\u003e0.6580\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003ea. This value corresponds to the apparent activation energy when adsorption pre-equilibration occurs: Δ\u003cem\u003eE\u003c/em\u003e\u0026thinsp;=\u0026thinsp;\u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,a\u003c/sub\u003e \u0026ndash; \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,d\u003c/sub\u003e + \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea,in\u003c/sub\u003e.\u003c/p\u003e\u003cp\u003eb. Critical dimensions of adsorbate molecule for entry into zeolite pores. In the slit-shaped pores, the size of the adsorption in the minimum dimension, MIN-1, was determined. In the cylindrical pores, the size of the molecule in two directions is considered: the minimum dimension, MIN-1, and the next smallest dimension, MIN-2.\u003c/p\u003e\u003cp\u003eTo further analyze the orientation of the adsorbed molecules, we defined two order parameters:\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"No\" id=\"Taba\" border=\"1\"\u003e\u003ccolgroup cols=\"2\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{S}_{\\perp\\:}=\\frac{1}{2}\\left(3⟨{\\text{cos}}^{2}\\alpha\\:⟩-1\\right)\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003e(3a)\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{S}_{\\parallel\\:}=\\frac{1}{2}\\left(3⟨{\\text{cos}}^{2}\\beta\\:⟩-1\\right)\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003e(3b)\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003ewhere \u003cem\u003eα\u003c/em\u003e and \u003cem\u003eβ\u003c/em\u003e defined the angles formed by the vectors perpendicular and parallel to the molecular plane, respectively, with respect to the perpendiculars passing through the center of the 6MR aperture. The characterization of these vectors and their orientation within the reference coordinate system are shown in Figure S7 (see SI). Figure\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e presents the order parameter for each molecule as it traverses the 6MR aperture. An animation depicting the molecular passage through the 6MR is provided in the SI. When the molecular orientation is isotropically distributed, the order parameter is zero. However, as each molecule passed through the aperture, both \u003cem\u003eS\u003c/em\u003e\u003csub\u003e\u0026perp;\u003c/sub\u003e and \u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e deviated from zero, indicating that the orientation of the molecules was constrained by the aperture. Notably, the negative value of \u003cem\u003eS\u003c/em\u003e\u003csub\u003e\u0026perp;\u003c/sub\u003e suggests that the molecules navigate the aperture in an orientation in which the molecular plane is nearly perpendicular to the aperture plane. This behavior was exhibited by all the adsorbate molecules. This indicated that the passage of adsorbed molecules through the 6MR aperture was accompanied by significant orientational constraints. In other words, the translational and rotational degrees of freedom of the adsorbed molecules were considerably restricted during their traversal through the 6MR aperture. Consequently, the observed significant negative activation entropy can be attributed to the transition of molecules in the vapor to a diffusion transition state, in which the molecular orientation is markedly constrained at the 6MR aperture.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eConversely, \u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e represents the variation in the rotational orientation of the molecules around an axis perpendicular to the molecular plane and exhibited different behaviors for each adsorbate molecule in the range \u0026minus;\u0026thinsp;3\u0026thinsp;\u0026le;\u0026thinsp;\u003cem\u003eZ\u003c/em\u003e/\u0026Aring; \u0026le; +3. For BZ, 1,4-CHD, CHEN, and CHEX, the values ranged from \u0026minus;\u0026thinsp;0.5\u0026thinsp;\u0026lt;\u0026thinsp;\u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e \u0026lt; +0.5, whereas for 1,3-CHD, 0\u0026thinsp;\u0026lt;\u0026thinsp;\u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e \u0026lt; +1. In BZ and 1,4-CHD, \u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e exhibited a bipolar change centered at \u003cem\u003eZ\u003c/em\u003e\u0026thinsp;=\u0026thinsp;0, indicating that the molecule turned its orientation around an axis perpendicular to the molecular plane and then turned it back. In the CHEN and CHEX, the variation in \u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e showed a similar trend, although the range of the \u003cem\u003eZ\u003c/em\u003e-axis was narrow. This feature implies that these molecules undergo twisting and wiggling around the axis perpendicular to the molecular plane when traversing the 6MR aperture (see the animation provided in the SI). In contrast, for 1,3-CHD, \u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e varied continuously from 0 to +\u0026thinsp;1 during the passage through 6MR, suggesting that the molecular orientation around the perpendicular axis changes progressively throughout the diffusion process. That is, the 1,3-CHD molecule rotates around the axis as it moves through the aperture (see the animation in the SI). This rotation is believed to minimize steric hindrance within the 6MR aperture by reducing friction between the molecule and the framework. This behavior differs from that of other molecules and may partly explain why 1,3-CHD exhibits a lower activation energy for diffusion than the other molecules. Distinguishing between 1,3-CHD and 1,4-CHD is crucial for understanding the factors influencing the diffusion coefficient of ZIF-8. The differences in the \u003cem\u003eE\u003c/em\u003e\u003csub\u003ea\u003c/sub\u003e and \u003cem\u003eS\u003c/em\u003e\u003csub\u003e||\u003c/sub\u003e behaviors observed in the MC simulations between 1,3-CHD and 1,4-CHD may reflect the slight differences in the molecular cross-sections and arrangement of π bonds within the molecules.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec14\" class=\"Section2\"\u003e\u003ch2\u003e3.3 Analysis of adsorbate-adsorbent interactions by spectroscopies\u003c/h2\u003e\u003cdiv id=\"Sec15\" class=\"Section3\"\u003e\u003ch2\u003e3.3.1 FT-IR spectroscopy\u003c/h2\u003e\u003cp\u003eFigure \u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e displays the symmetric C\u0026ndash;H stretching bands associated with the 2-methylimidazole group. The full wavenumber regions of the infrared spectra are shown in Figure S8 (SI). Notably, no shift was observed in the band upon cyclohexane adsorption, while a pronounced redshift was observed over the wavenumber range of 1\u0026ndash;2 cm\u003csup\u003e‒1\u003c/sup\u003e upon benzene and cyclic alkene adsorption. This shift was more pronounced for 1,3-CHD, 1,4-CHD, and BZ than for CHEN. This observation suggests the elongation of C\u0026ndash;H bond, which is attributed to the attractive interaction between the π-electrons of the adsorbed molecules and the C\u0026ndash;H bond of 2-methylimidazole [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. One plausible attractive interaction is the CH-π interaction, which has been reported for various organic compounds, metal complexes, and molecular assemblies [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e]. It has been established that CH-π bonding induces stretching of the C\u0026ndash;H bonds, resulting in a redshift of the CH stretching band [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e, \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. Additionally, the protons of the methyl group may function as proton donors in CH-π interactions, thereby leading to an anticipated redshift of the CH symmetric stretching band in the methyl group of the 2-methylimidazolate moiety [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. However, the symmetric CH stretching band of the methyl group overlaps with the CH stretching bands of the adsorbate molecules. Deconvolution of the infrared bands potentially confirmed the redshift of the CH symmetric stretching band of the methyl group. However, the accuracy of the resultant redshift from this deconvolution was insufficient for a meaningful discussion of the minor wavenumber shifts observed in this study. The elongation of the C\u0026ndash;H bonds and the corresponding redshift of the CH symmetric stretching band were also qualitatively corroborated by the molecular assembly of 2-methylimidazole with benzene, 1,3-CHD, and 1,4-CHD, optimized using DFT calculations (Figure S9, Figure S10, and Table \u003cspan refid=\"MOESM3\" class=\"InternalRef\"\u003eS3\u003c/span\u003e). Generally, intermolecular CH-π interactions occur within 3 \u0026Aring; [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e] and exhibit a high degree of directionality, with the orientation of the C\u0026ndash;H bond aligned with that of the π-orbital. Consequently, the C\u0026ndash;H and π bonds tend to be oriented nearly perpendicular to each other. Following the formation of the CH-π bond, the adsorbed molecules were arranged such that their molecular planes were positioned at an angle that was nearly perpendicular to the 6MR aperture plane. This orientation facilitates the ingress of adsorbed molecules into the 6MR apertures. Thus, the CH-π interaction between the 2-methylimidazole ring and the adsorbed molecules is crucial for governing the molecular orientation necessary for the effective passage through the 6MR aperture for benzene and cyclic alkenes.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec16\" class=\"Section3\"\u003e\u003ch2\u003e3.3.2 \u003csup\u003e13\u003c/sup\u003eC CP/MAS NMR spectroscopy\u003c/h2\u003e\u003cp\u003eThe \u003csup\u003e13\u003c/sup\u003eC CP/MAS NMR spectra of ZIF-8, both in its pure form and with various adsorbates are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e7\u003c/span\u003e. Cross-polarization (CP) measurements are particularly effective for highlighting immobile molecular segments within solid matrices while attenuating signals from more mobile components. Consequently, the highly mobile adsorbed molecules appeared as minor signals in the spectra. The predominant peaks are attributed to the 2-methylimidazole moiety, suggesting that the ZIF-8 framework exhibited greater rigidity and reduced mobility than the adsorbate molecules. Specifically, the resonance lines at approximately 14, 124, and 151 ppm correspond to the methyl group, CH carbon of the imidazole ring, and quaternary carbon of the 2-methylimidazolate moiety, respectively. Figures\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e8\u003c/span\u003e(a), 8(b), and 8(c) provide enlarged views of these spectral components. Notably, the adsorption of benzene and cyclic alkenes resulted in downfield shifts of 0.5\u0026ndash;1.3 ppm in the resonance lines associated with the methyl group and CH carbon of the imidazole ring, whereas the adsorption of cyclohexane did not elicit such shifts. The resonance lines corresponding to quaternary carbons remained largely unchanged after adsorption. These observations suggest a significant interaction between the π-electrons of the adsorbed molecules and both the methyl group and the CH carbon of the 2-methylimidazolate moiety. Our previous studies indicated that both the methyl and CH protons of the imidazole ring can act as proton donors in CH-π bond. According to Pople's point magnetic dipole model [\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e, \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e], the shielding effect of carbon along a perpendicular line through the center of the molecular plane of the benzene ring is described by the equation Δ\u003cem\u003eσ\u003c/em\u003e\u003csub\u003eC\u003c/sub\u003e\u0026thinsp;=\u0026thinsp;51.96/\u003cem\u003eR\u003c/em\u003e\u003csup\u003e3\u003c/sup\u003e, where \u003cem\u003eR\u003c/em\u003e represents the average distance between the carbon nucleus of interest and the molecular plane in angstroms (\u0026Aring;). For the optimized structure of the 2-methylimidazole-benzene complex derived from DFT calculations (see Figure S9), \u003cem\u003eR\u003c/em\u003e was determined to be 3.27 \u0026Aring;, which yielded a shielding effect of approximately 1.5 ppm. In the case of CH-π bond, the C-H bond is anticipated to elongate owing to the attraction of the proton by the π electrons, resulting in an anti-shielding effect on the \u003csup\u003e13\u003c/sup\u003eC nucleus. Scheiner reported an anti-shielding effect of 200 ppm/\u0026Aring; for the tetrel bond of a methyl group [\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e]. To account for the observed downfield shift, an average elongation of approximately 0.01\u0026ndash;0.015 \u0026Aring; in the C-H bond is necessary, which is consistent with the expected bond-length increase associated with a strong CH-π bond. In summary, the downfield shifts in the resonance lines resulting from the adsorption of benzene and cyclic alkenes can be attributed to the CH-π bond between the CH protons of the methyl group and imidazole ring and π electrons of the adsorbed molecules.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Sec17\" class=\"Section2\"\u003e\u003ch2\u003e3.4 Local structure of molecular passage through 6MR\u003c/h2\u003e\u003cp\u003eThis study elucidates the intermolecular interactions that dictate molecular orientation during the adsorption of bulky molecules onto ZIF-8, particularly as these molecules navigate through 6MR apertures. The results suggest a specific configuration for benzene and six-membered-ring alicyclic hydrocarbons as they traverse the 6MR aperture (Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e9\u003c/span\u003e). The CH-π bond are pivotal in influencing the molecular orientation of benzene and cyclic alkenes during their passage through the aperture. As these molecules approach the CH bottleneck of the 6MR, they form CH-π bonds with the CH groups of the 2-methylimidazole ring at this bottleneck. This highly directional CH-π bonding results in an nearly perpendicular (T-shaped) alignment of the imidazole ring and adsorbed molecule. This arrangement ensured that the adsorbed molecules maintained a configuration that minimized their molecular cross-sectional areas (molecular planes) upon entering the 6MR aperture. After successfully passing through the aperture, the adsorbed molecule is stabilized by forming an additional CH-π bond with the methyl group on the methyl bottleneck side (Figure S10).\u003c/p\u003e\u003cp\u003eIn contrast, cyclic alkanes such as cyclohexane primarily engage in dispersion-force interactions. The optimal molecular orientation for enhancing these interactions occurs when the plane of the adsorbed molecule is parallel to the 2-methylimidazole ring, as this alignment maximizes the contact area between them. Consequently, a molecule approaching the CH bottleneck side of 6MR is captured by dispersion force interactions with the methyl group and quaternary carbon in the 2-methylimidazolate moieties. This results in the adsorbed molecule orienting its minimum molecular cross-section (molecular face) towards the 6MR aperture, thereby facilitating its entry. After passing through the 6MR aperture, the molecule stabilizes through dispersive force interactions with the three methyl groups on the methyl bottleneck side. Cyclohexene, being a cyclic alkene, also exhibits significant contributions from dispersion force interactions due to its substantial proportion of saturated hydrocarbon components, in addition to the CH-π interactions (see Figure S10).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eFinally, ZIF-8 exhibited remarkably different diffusion coefficients for 1.3-CHD and 1,4-CHD. This result is of particular interest because it highlights the potential of ZIF-8 as a molecular sieve for isomer separation. Further investigation is necessary to elucidate these detailed mechanisms. However, the isomer distinction may involve the symmetry of the phonon modes associated with the swing effect of the bridging ligands that constitute the 6MR, as well as the symmetry of the adsorbate molecule. Additionally, considering the aspect of configuration diffusion, this isomer distinction might be linked to the internal arrangement of π electrons through the CH-π bond, which controls the molecular orientation with respect to the 6MR aperture.\u003c/p\u003e\u003c/div\u003e"},{"header":"4 Conclusion","content":"\u003cp\u003eIn this study, we investigated the adsorption kinetics of benzene and six-membered-ring alicyclic hydrocarbons on ZIF-8, to clarify the interactions between the adsorbed molecules and linkers in relation to the adsorption behavior of larger molecules. Monte Carlo simulations were employed to elucidate the local structural dynamics of the molecules as they navigated through the six-membered ring aperture of ZIF-8. Additionally, the FT-IR and \u003csup\u003e13\u003c/sup\u003eC cross-polarization magic angle spinning nuclear magnetic resonance (CP/MAS NMR) were used to identify the adsorbate-adsorbate interactions that influence molecular orientation during transit through the aperture.\u003c/p\u003e\u003cp\u003eThe findings indicate that the diffusion of benzene and other six-membered ring hydrocarbons through the six-membered ring aperture is predominantly governed by a significant negative activation entropy ranging from \u0026minus;\u0026thinsp;172 to \u0026minus;\u0026thinsp;217 J K\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e mol\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. This observation suggests a pronounced effect of the molecular orientation on the transit of adsorbate molecules through the apertures. Monte Carlo simulations further elucidated this phenomenon, revealing that the aperture diameter significantly affected the activation energy required for the molecular passage. Order parameter analysis revealed that the molecular orientation was tightly constrained and aligned in a manner that minimizes the energy barrier for traversing the aperture. Spectroscopic analyses employing FT-IR and \u003csup\u003e13\u003c/sup\u003eC CP/MAS NMR techniques identified a specific CH-π interaction between the π electrons of the adsorbed molecules and the methyl group and imidazole ring CH of the 2-methylimidazolate linker. For unsaturated cyclic hydrocarbons possessing π electrons, the CH-π interaction predominantly governs the molecular orientation during aperture transit. Conversely, saturated cyclic hydrocarbons are primarily influenced by dispersion force interactions with the methyl groups of the bridging ligand. These orientation-controlling interactions are anticipated to play a critical role in governing molecular adsorption and diffusion behaviors of other metal-organic frameworks (MOFs), porous coordination polymers (PCPs), and porous organic materials (POMs) composed of organic materials.\u003c/p\u003e\u003cp\u003eFinally, the distinct diffusion coefficients of 1.3-CHD and 1,4-CHD are of particular interest. This indicated the potential of ZIF-8 as a molecular sieve for structural isomers.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors thank Dr. Naoya Inazumi and Dr. Yasuto Todokoro of the Analytical Instrument Facility, Graduate School of Science, Osaka University, for their helpful advice, guidance, and instructions regarding the solid-state NMR measurements. We would like to thank Editage (www.editage.jp) for English language editing.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor Contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTakahiro Ueda:\u003c/strong\u003e Conceptualization, methodology, validation, resources, software, formal analysis, visualization, writing\u0026ndash;original draft, review, editing, supervision, funding acquisition, and project administration. \u003cstrong\u003eYuta Yamada:\u003c/strong\u003e Methodology, investigation, data curation, and visualization. \u003cstrong\u003eKota Fujii:\u0026nbsp;\u003c/strong\u003eInvestigation, DFT calculations, data curation, and visualization.\u003cstrong\u003e\u0026nbsp;Ryota Mihara:\u0026nbsp;\u003c/strong\u003eInvestigation, Monte Carlo simulation, data curation, and visualization.\u003cstrong\u003e\u0026nbsp;Taku Iiyama:\u003c/strong\u003e Methodology, investigation, software, data curation, validation, visualization, review, and editing. \u003cstrong\u003eYasutaka Hamada:\u003c/strong\u003e Investigation, DFT calculations, data curation, and visualization. \u003cstrong\u003eMitsutaka Okumura:\u0026nbsp;\u003c/strong\u003eMethodology, software, and supervision.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis research was supported by JSPS KAKENHI Grant-in-Aid for Scientific Research (C) (Grant Numbers JP21K04979, JP24K8360). This work was the result of using research equipment shared in the MEXT Project to promote the public utilization of advanced research infrastructure (Program for supporting the construction of core facilities), Grant Numbers JPMXS0441200021, JPMXS0441200023, and JPMXS0441200024.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe data presented in this article will be available from the corresponding author on reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthical Approval\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConflicts of interest\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no conflicts of interest.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eFan, L., Tu, Z., Chan, S.H.: Recent development of hydrogen and fuel cell technologies: A review. Energy Rep. \u003cb\u003e7\u003c/b\u003e, 8421\u0026ndash;8446 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.egyr.2021.08.003\u003c/span\u003e\u003cspan address=\"10.1016/j.egyr.2021.08.003\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eCrandall, B.S., Brix, T., Weber, R.S., Jiao, F.: Techno-Economic Assessment of Green H\u003csub\u003e2\u003c/sub\u003e Carrier Supply Chains. Energy Fuels. \u003cb\u003e37\u003c/b\u003e, 1441\u0026ndash;1450 (2023). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acs.energyfuels.2c03616\u003c/span\u003e\u003cspan address=\"10.1021/acs.energyfuels.2c03616\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSalman, M.S., Rambhujun, N., Pratthana, C., Srivastava, K., Aguey-Zinsou, K.-F.: Catalysis in Liquid Organic Hydrogen Storage: Recent Advances, Challenges, and Perspectives. Ind. Eng. Chem. Res. \u003cb\u003e61\u003c/b\u003e, 6067\u0026ndash;6105 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acs.iecr.1c03970\u003c/span\u003e\u003cspan address=\"10.1021/acs.iecr.1c03970\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSwesi, Y., Kerleau, P., Pitault, I., Heurtaux, F., Ronze, D.: Purification of hydrogen from hydrocarbons by adsorption for vehicles application. Sep. Purif. Technol. \u003cb\u003e56\u003c/b\u003e, 25\u0026ndash;37 (2007). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.seppur.2006.09.023\u003c/span\u003e\u003cspan address=\"10.1016/j.seppur.2006.09.023\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eAcharya, D., Ng, D., Xie, Z.: Recent Advances in Catalysts and Membranes for MCH Dehydrogenation: A Mini Review. Membranes. \u003cb\u003e11\u003c/b\u003e, 955 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3390/membranes11120955\u003c/span\u003e\u003cspan address=\"10.3390/membranes11120955\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eMeng, J., Zhou, F., Ma, H., Yuan, X., Wang, Y., Zhang, J.: A Review of Catalysts for Methylcyclohexane Dehydrogenation. Top. Catal. \u003cb\u003e64\u003c/b\u003e, 509\u0026ndash;520 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s11244-021-01465-6\u003c/span\u003e\u003cspan address=\"10.1007/s11244-021-01465-6\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eChen, L., Verma, P., Hou, K., Qi, Z., Zhang, S., Liu, Y.-S., Guo, J., Stavila, V., Allendorf, M.D., Zheng, L., Salmeron, M., Prendergast, D., Somorjai, G.A., Su, J.: Reversible dehydrogenation and rehydrogenation of cyclohexane and methylcyclohexane by single-site platinum catalyst. Nat. Commun., 13, 1092 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.nature.com/articles/s41467-022-28607-y\u003c/span\u003e\u003cspan address=\"https://www.nature.com/articles/s41467-022-28607-y\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZhokh, A., Serebrii, T., Strizhak, P.: Diffusion-controlled adsorption\u0026ndash;desorption of C6 hydrocarbons by multiwall carbon nanotubes. Chem. Phys. Lett. \u003cb\u003e805\u003c/b\u003e, 139935 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.cplett.2022.139935\u003c/span\u003e\u003cspan address=\"10.1016/j.cplett.2022.139935\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eMukherjee, S., Sensharma, D., Qazvini, O.T., Dutta, S., Macreadie, L.K., Ghosh, S.K., Babarao, R.: Advances in adsorptive separation of benzene and cyclohexane by metal-organic framework adsorbents. Coord. Chem. Rev. \u003cb\u003e437\u003c/b\u003e, 213852 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.ccr.2021.213852\u003c/span\u003e\u003cspan address=\"10.1016/j.ccr.2021.213852\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eGraczov\u0026aacute;, E., Moln\u0026aacute;r, D., Steltenpohl, P.: Separation of Methylcyclohexane-Toluene Mixture Using Imidazolium Ionic Liquids. Chem. Eng. Trans. \u003cb\u003e94\u003c/b\u003e, 745\u0026ndash;750 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3303/CET2294124\u003c/span\u003e\u003cspan address=\"10.3303/CET2294124\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eValencia, A., Mu\u0026ntilde;iz-Valencia, R., Ceballos-Maga\u0026ntilde;a, S.G., Rojas-Mayorga, C.K., Bonilla-Petriciolet, A., Gonz\u0026aacute;lez, J., Aguayo-Villarreal, I.A.: Cyclohexane and benzene separation by fixed-bed adsorption on activated carbons prepared from coconut shell. Environ. Technol. Innov. \u003cb\u003e25\u003c/b\u003e, 102076 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.eti.2021.102076\u003c/span\u003e\u003cspan address=\"10.1016/j.eti.2021.102076\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eChen, C., Guan, H., Li, H., Zhou, Y., Huang, Y., Wei, W., Hong, M., Wu, M.: A Noncovalent π-Stacked Porous Organic Molecular Framework for Selective Separation of Aromatics and Cyclic Aliphatics. Angew Chem. Int. Ed. \u003cb\u003e61\u003c/b\u003e, e202201646 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/anie.202201646\u003c/span\u003e\u003cspan address=\"10.1002/anie.202201646\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003ePhan, A., Doonan, C.J., Uribe-Romo, F.J., Knobler, C.B., O\u0026rsquo;Keeffe, M., Yaghi, O.Y.: Synthesis, Structure, and Carbon Dioxide Capture Properties of Zeolitic Imidazolate Frameworks. Acc. Chem. Res. \u003cb\u003e43\u003c/b\u003e, 58\u0026ndash;67 (2010). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/ar900116g\u003c/span\u003e\u003cspan address=\"10.1021/ar900116g\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eBanerjee, R., Phan, A., Wang, B., Knobler, C., Furukawa, H., O\u0026rsquo;Keeffe, M., Yaghi, O.M.: High-Throughput Synthesis of Zeolitic Imidazolate Frameworks and Application to CO\u003csub\u003e2\u003c/sub\u003e Capture. Science. \u003cb\u003e319\u003c/b\u003e, 939\u0026ndash;943 (2008). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1126/science.1152516\u003c/span\u003e\u003cspan address=\"10.1126/science.1152516\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eCheng, P., Hu, Y.H.: H\u003csub\u003e2\u003c/sub\u003eO-Functionalized Zeolitic Zn(2-methylimidazole)\u003csub\u003e2\u003c/sub\u003e Framework (ZIF-8) for H\u003csub\u003e2\u003c/sub\u003e Storage. J. Phys. Chem. C. \u003cb\u003e118\u003c/b\u003e, 21866\u0026ndash;21872 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://pubs.acs.org/doi/\u003c/span\u003e\u003cspan address=\"https://pubs.acs.org/doi/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1021/jp507030g\u003c/span\u003e\u003cspan address=\"10.1021/jp507030g\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eB\u0026ouml;hme, U., Barth, B., Paula, C., Kuhnt, A., Schwieger, W., Mundstock, A., Caro, J., Hartmann, M.: Ethene/Ethane and Propene/Propane Separation via the Olefin and Paraffin Selective Metal\u0026ndash;Organic Framework Adsorbents CPO-27 and ZIF-8. Langmuir. \u003cb\u003e29\u003c/b\u003e, 8592\u0026ndash;8600 (2013). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/la401471g\u003c/span\u003e\u003cspan address=\"10.1021/la401471g\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZhang, C., Lively, R.P., Zhang, K., Johnson, J.R., Karvan, O., Koros, W.J.: Unexpected Molecular Sieving Properties of Zeolitic Imidazolate Framework-8. J. Phys. Chem. Lett. \u003cb\u003e3\u003c/b\u003e, 2130\u0026ndash;2134 (2012). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jz300855a\u003c/span\u003e\u003cspan address=\"10.1021/jz300855a\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZhang, K., Lively, R.P., Zhang, C., Chance, R.R., Koros, W.J., Sholl, D.S., Nair, S.: Exploring the Framework Hydrophobicity and Flexibility of ZIF-8: From Biofuel Recovery to Hydrocarbon Separations. J. Phys. Chem. Lett. \u003cb\u003e4\u003c/b\u003e, 3618\u0026ndash;3622 (2013). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jz402019d\u003c/span\u003e\u003cspan address=\"10.1021/jz402019d\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eFairen-Jimenez, D., Moggach, S.A., Wharmby, M.T., Wright, P.A., Parsons, S., D\u0026uuml;ren, T.: Opening the Gate: Framework Flexibility in ZIF-8 Explored by Experiments and Simulations. J. Am. Chem. Soc. \u003cb\u003e133\u003c/b\u003e, 8900\u0026ndash;8902 (2011). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/ja202154j\u003c/span\u003e\u003cspan address=\"10.1021/ja202154j\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTian, T., Wharmby, M.T., Parra, J.B., Ania, C.O., Fairen-Jimenez, D.: Role of crystal size on swing-effect and adsorption induced structure transition of ZIF-8. Dalton Trans. \u003cb\u003e45\u003c/b\u003e, 6893\u0026ndash;6900 (2016). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1039/C6DT00565A\u003c/span\u003e\u003cspan address=\"10.1039/C6DT00565A\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eCoudert, F.-X.: Molecular Mechanism of Swing Effect in Zeolitic Imidazolate Framework ZIF-8: Continuous Deformation upon Adsorption. ChemPhysChem. \u003cb\u003e18\u003c/b\u003e, 2732\u0026ndash;2738 (2017). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/cphc.201700463\u003c/span\u003e\u003cspan address=\"10.1002/cphc.201700463\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eMorris, W., Stevens, C.J., Taylor, R.E., Dybowski, C., Yaghi, O.M., Garcia-Garibay, M.A.: NMR and X-ray Study Revealing the Rigidity of Zeolitic Imidazolate Frameworks. J. Phys. Chem. C. \u003cb\u003e116\u003c/b\u003e, 13307\u0026ndash;13312 (2012). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jp303907p\u003c/span\u003e\u003cspan address=\"10.1021/jp303907p\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTanaka, H., Ohsaki, S., Hiraide, S., Yamamoto, D., Watanabe, S., Miyahara, M.T.: Adsorption-Induced Structural Transition of ZIF-8: A Combined Experimental and Simulation Study. J. Phys. Chem. C. \u003cb\u003e118\u003c/b\u003e, 8445\u0026ndash;8454 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jp500931g\u003c/span\u003e\u003cspan address=\"10.1021/jp500931g\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eAnia, C.O., Garc\u0026iacute;a-P\u0026eacute;rez, E., Haro, M., Guti\u0026eacute;rrez-Sevillano, J.J., Vald\u0026eacute;s-Sol\u0026iacute;s, T., Parra, J.B., Calero, S.: Understanding Gas-Induced Structural Deformation of ZIF-8. J. Phys. Chem. Lett. \u003cb\u003e3\u003c/b\u003e, 1159\u0026ndash;1164 (2012). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jz300292y\u003c/span\u003e\u003cspan address=\"10.1021/jz300292y\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eRyder, M.R., Civalleri, B., Bennett, T.D., Henke, S., Rudić, S., Cinque, G., Fernandez-Alonso, F., Tan, J.-C.: Identifying the Role of Terahertz Vibrations in Metal-Organic Frameworks: From Gate-Opening Phenomenon to Shear-Driven Structural Destabilization. Phys. Rev. Lett. \u003cb\u003e113\u003c/b\u003e, 215502 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/PhysRevLett.113.215502\u003c/span\u003e\u003cspan address=\"10.1103/PhysRevLett.113.215502\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eCasco, M.E., Cheng, Y.Q., Daemen, L.L., Fairen-Jimenez, D., Ramos-Fern\u0026aacute;ndez, E.V., Ramirez-Cuesta, A.J., Silvestre-Albero, J.: Gate-opening effect in ZIF-8: the first experimental proof using inelastic neutron scattering. Chem. Commun. \u003cb\u003e52\u003c/b\u003e, 3639\u0026ndash;3642 (2016). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1039/C5CC10222G\u003c/span\u003e\u003cspan address=\"10.1039/C5CC10222G\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTan, N.Y., Ruggiero, M.T., Orellana-Tavra, C., Tian, T., Bond, A.D., Korter, T.M., Fairen-Jimenez, D., Zeitler, J.A.: Investigation of the terahertz vibrational modes of ZIF-8 and ZIF-90 with terahertz time-domain spectroscopy. Chem. Commun. \u003cb\u003e51\u003c/b\u003e, 16037\u0026ndash;16040 (2015). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1039/C5CC06455D\u003c/span\u003e\u003cspan address=\"10.1039/C5CC06455D\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eKolokolov, D.I., Stepanov, A.G., Jobic, H.: Mobility of the 2-Methylimidazolate Linkers in ZIF-8 Probed by \u003csup\u003e2\u003c/sup\u003eH NMR: Saloon Doors for the Guests. J. Phys. Chem. C. \u003cb\u003e119\u003c/b\u003e, 27512\u0026ndash;27520 (2015). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acs.jpcc.5b09312\u003c/span\u003e\u003cspan address=\"10.1021/acs.jpcc.5b09312\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eUeda, T., Nakai, M., Yamatani, T.: A solid-state \u003csup\u003e1\u003c/sup\u003eH-NMR study of the dynamic structure of ZIF-8 and its role in the adsorption of bulky molecules. Adsorption. \u003cb\u003e23\u003c/b\u003e, 887\u0026ndash;901 (2017). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s10450-017-9902-5\u003c/span\u003e\u003cspan address=\"10.1007/s10450-017-9902-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eUeda, T., Yamatani, T., Okumura, M.: Dynamic Gate Opening of ZIF-8 for Bulky Molecule Adsorption as Studied by Vapor Adsorption Measurements and Computational Approach. J. Phys. Chem. C. \u003cb\u003e123\u003c/b\u003e, 27542\u0026ndash;27553 (2019). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acs.jpcc.9b07239\u003c/span\u003e\u003cspan address=\"10.1021/acs.jpcc.9b07239\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eYagi, R., Ueda, T.: Substitution (CH\u003csub\u003e3\u003c/sub\u003e, Cl, or Br) effects of the imidazolate linker on benzene adsorption kinetics for the zeolitic imidazolate framework (ZIF)-8. Phys. Chem. Chem. Phys. \u003cb\u003e25\u003c/b\u003e, 20585\u0026ndash;20596 (2023). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1039/D3CP01662E\u003c/span\u003e\u003cspan address=\"10.1039/D3CP01662E\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003ePoryvaev, A.S., Yazikova, A.A., Polyukhov, D.M., Fedin, M.V.: Ultrahigh selectivity of benzene/cyclohexane separation by ZIF-8 framework: Insights from spin-probe EPR spectroscopy. Micropor Mesopor Mater. \u003cb\u003e330\u003c/b\u003e, 111564 (2022). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.micromeso.2021.111564\u003c/span\u003e\u003cspan address=\"10.1016/j.micromeso.2021.111564\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eNishio, M.: CH/π hydrogen bonds in crystals. CrystEngComm. \u003cb\u003e6\u003c/b\u003e, 130\u0026ndash;158 (2004). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1039/B313104A\u003c/span\u003e\u003cspan address=\"10.1039/B313104A\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eScheiner, S.: Origins and properties of the tetrel bond. Phys. Chem. Chem. Phys. \u003cb\u003e23\u003c/b\u003e, 5702\u0026ndash;5717 (2021). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1039/D1CP00242B\u003c/span\u003e\u003cspan address=\"10.1039/D1CP00242B\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eParkes, M.V., Demir, H., Teich-McGoldrick, S.L., Sholl, D.S., Greathouse, J.A., Allendorf, M.D.: Molecular dynamics simulation of framework flexibility effects on noble gas diffusion in HKUST-1 and ZIF-8. Micropor Mesopor Mater. \u003cb\u003e194\u003c/b\u003e, 190\u0026ndash;199 (2014). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.micromeso.2014.03.027\u003c/span\u003e\u003cspan address=\"10.1016/j.micromeso.2014.03.027\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eZheng, B., Sant, M., Demontis, P., Suffritti, G.B.: Force Field for Molecular Dynamics Computations in Flexible ZIF-8 Framework. J. Phys. Chem. C. \u003cb\u003e116\u003c/b\u003e, 933\u0026ndash;938 (2012). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jp209463a\u003c/span\u003e\u003cspan address=\"10.1021/jp209463a\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eFrisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Scalmani, G., Barone, V., Petersson, G.A., Nakatsuji, H., Li, X., Caricato, M., Marenich, A.V., Bloino, J., Janesko, B.G., Gomperts, R., Mennucci, B., Hratchian, H.P., Ortiz, J.V., Izmaylov, A.F., Sonnenberg, J.L., Williams-Young, D., Ding, F., Lipparini, F., Egidi, F., Goings, J., Peng, B., Petrone, A., Henderson, T., Ranasinghe, D., Zakrzewski, V.G., Gao, J., Rega, N., Zheng, G., Liang, W., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Throssell, K., Montgomery, J.A. Jr., Peralta, J.E., Ogliaro, F., Bearpark, M.J., Heyd, J.J., Brothers, E.N., Kudin, K.N., Staroverov, V.N., Keith, T.A., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A.P., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Millam, J.M., Klene, M., Adamo, C., Cammi, R., Ochterski, J.W., Martin, R.L.: (2019). Morokuma, K., Farkas, O., Foresman, J. B., Fox, D. J.: Gaussian 16, Revision C.01, Gaussian, Inc., Wallingford CT\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eJorgensen, W.L., Tirado-Rives, J.: The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J. Am. Chem. Soc. \u003cb\u003e110\u003c/b\u003e, 1657\u0026ndash;1666 (1988). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://pubs.acs.org\u003c/span\u003e\u003cspan address=\"https://pubs.acs.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e/doi/10.1021/ja00214a001\u003c/span\u003e\u003cspan address=\"/doi/10.1021/ja00214a001\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eJorgensen, W.L., Maxwell, D.S., Tirado-Rives, J.: Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. \u003cb\u003e118\u003c/b\u003e, 11225\u0026ndash;11236 (1996). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/ja9621760\u003c/span\u003e\u003cspan address=\"10.1021/ja9621760\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eK\u0026auml;rger, J., Ruthven, D.M., Theodorou, D.N.: Diffusion in nanoporous materials. Wiley-VCH, Weinheim, Germany (2012)\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eBreck, D.W.: Zeolite Molecular Sieves, pp. 671\u0026ndash;688. John Willey \u0026amp; Sons, New York (1974)\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eGlasstone, S., Laidler, K.J., Eyring, H.: The Theory of Rate Processes: The Kinetics of Chemical Reactions, Viscosity, Diffusion and Electrochemical Phenomena. McGraw-Hill, New York, USA (1941)\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eWebster, C.E., Drago, R.S., Zerner, M.C.: Molecular dimensions for adsorptives. J. Am. Chem. Soc. \u003cb\u003e120\u003c/b\u003e, 5509\u0026ndash;5516 (1998). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/ja973906m\u003c/span\u003e\u003cspan address=\"10.1021/ja973906m\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eFujii, A., Morita, S., Miyazaki, M., Ebata, T., Mikami, N.: A Molecular Cluster Study on Activated CH/π Interactions: Infrared Spectroscopy of Aromatic Molecule\u0026thinsp;\u0026ndash;\u0026thinsp;Acetylene Clusters. J. Phys. Chem. A. \u003cb\u003e108\u003c/b\u003e(14), 2652\u0026ndash;2658 (2004). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jp049946b\u003c/span\u003e\u003cspan address=\"10.1021/jp049946b\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eTsuzuki, S., Fujii, A.: Nature and physical origin of CH/π interaction: significant difference from conventional hydrogen bonds. Phys. Chem. Chem. Phys. \u003cb\u003e10\u003c/b\u003e, 2584\u0026ndash;2594 (2008). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1039/B718656H\u003c/span\u003e\u003cspan address=\"10.1039/B718656H\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003ePople, J.A.: Proton Magnetic Resonance of Hydrocarbons. J. Chem. Phys. \u003cb\u003e24\u003c/b\u003e, 1111 (1956). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1063/1.1742701\u003c/span\u003e\u003cspan address=\"10.1063/1.1742701\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eSahakyan, A.B., Vendruscolo, M.: Analysis of the Contributions of Ring Current and Electric Field Effects to the Chemical Shifts of RNA Bases. J. Phys. Chem. B. \u003cb\u003e117\u003c/b\u003e, 1989\u0026ndash;1998 (2013). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/jp3057306\u003c/span\u003e\u003cspan address=\"10.1021/jp3057306\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003cli\u003e\u003cspan\u003eScheiner, S.: Dependence of NMR chemical shifts upon CH bond lengths of a methyl group involved in a tetrel bond. Chem. Phys. Lett. \u003cb\u003e714\u003c/b\u003e, 61\u0026ndash;64 (2019). \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.cplett.2018.10.069\u003c/span\u003e\u003cspan address=\"10.1016/j.cplett.2018.10.069\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"adsorption","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"adso","sideBox":"Learn more about [Adsorption](http://link.springer.com/journal/10450)","snPcode":"10450","submissionUrl":"https://submission.nature.com/new-submission/10450/3","title":"Adsorption","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Adsorption kinetics, Molecular sieving, Diffusion, Activation entropy, CH-π interaction","lastPublishedDoi":"10.21203/rs.3.rs-7276006/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-7276006/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThis study examined the adsorption kinetics of benzene and six-membered ring alicyclic hydrocarbons on Zeolitic imidazolate framework-8 (ZIF-8) to elucidate the interactions between the adsorbed molecules and linkers on the adsorption behavior of bulky molecules. The temperature dependence of the Fickian diffusion coefficient was analyzed to determine the activation entropy and energy of diffusion. The significant negative values observed for the activation entropy, from \u0026minus;\u0026thinsp;172 to \u0026minus;\u0026thinsp;217 J K\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e mol\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e, indicate a restriction on the degree of freedom of the adsorbate molecules within the diffusion transition state, thereby reducing the diffusion coefficients. The ZIF-8 differentiates between 1,3- and 1,4-cyclohexadiene based on kinetic mechanisms, indicating its potential as an isomeric molecular sieve. The molecular orientations of the adsorbate molecules traversing the ZIF-8 aperture were investigated using Monte Carlo simulations. Order parameter analysis revealed that the molecular orientation was constrained to minimize the energy barrier with passage through the aperture. Fourier-transform infrared (FT-IR) spectroscopy and \u003csup\u003e13\u003c/sup\u003eC cross-polarization magic angle spinning (CP/MAS) nuclear magnetic resonance (NMR) measurements supported the CH-π interactions between the adsorbate molecules and 2-methylimidazolate linkers. This interaction maintained the orientation of the adsorbate molecules, facilitating their passage through the ZIF-8 apertures. CH-π interactions influence molecular orientation control in the adsorption and diffusion of large molecules containing π electrons in porous materials with organic ligands.\u003c/p\u003e","manuscriptTitle":"Intermolecular interactions governing adsorption of bulky molecules on ZIF-8: Insight from adsorption kinetics of benzene and 6-membered ring alicyclic hydrocarbons","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-08-27 16:07:58","doi":"10.21203/rs.3.rs-7276006/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-09-15T00:12:59+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-09-04T10:05:54+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-09-01T13:03:45+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-08-27T15:10:15+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"148630918497182385164611962797482300246","date":"2025-08-21T01:54:39+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"101463799794143186919824091087694111850","date":"2025-08-19T13:36:24+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"320710133396725053944870047385331159145","date":"2025-08-19T12:38:06+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-08-19T11:06:53+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-08-04T12:13:16+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-08-04T11:37:12+00:00","index":"","fulltext":""},{"type":"submitted","content":"Adsorption","date":"2025-08-02T06:23:40+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"adsorption","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"adso","sideBox":"Learn more about [Adsorption](http://link.springer.com/journal/10450)","snPcode":"10450","submissionUrl":"https://submission.nature.com/new-submission/10450/3","title":"Adsorption","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"0d19d185-8524-419e-afab-e24bf20f71cd","owner":[],"postedDate":"August 27th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2026-03-23T16:02:44+00:00","versionOfRecord":{"articleIdentity":"rs-7276006","link":"https://doi.org/10.1007/s10450-026-00675-2","journal":{"identity":"adsorption","isVorOnly":false,"title":"Adsorption"},"publishedOn":"2026-03-19 15:58:35","publishedOnDateReadable":"March 19th, 2026"},"versionCreatedAt":"2025-08-27 16:07:58","video":"","vorDoi":"10.1007/s10450-026-00675-2","vorDoiUrl":"https://doi.org/10.1007/s10450-026-00675-2","workflowStages":[]},"version":"v1","identity":"rs-7276006","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-7276006","identity":"rs-7276006","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}