
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

38977694
49822
10.1038/s41467-024-49822-9
Matters Arising
Reply to: Diffusion anomaly in nanopores as a rich field for theorists and a challenge for experimentalists
http://orcid.org/0000-0002-7143-2658
Gao Mingbin 1
Yuan Jiamin 2
http://orcid.org/0000-0003-2872-0125
Liu Zhiqiang 2
http://orcid.org/0000-0002-7078-2402
Ye Mao maoye@dicp.ac.cn

1
http://orcid.org/0000-0001-7115-6510
Zheng Anmin zhenganm@wipm.ac.cn

2
1 grid.9227.e 0000000119573309 National Engineering Research Center of Lower-Carbon Catalysis Technology, Dalian National Laboratory for Clean Energy, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, People’s Republic of China
2 grid.9227.e 0000000119573309 State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, National Center for Magnetic Resonance in Wuhan, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan, People’s Republic of China
8 7 2024
8 7 2024
2024
15 57223 10 2023
18 6 2024
© The Author(s) 2024, corrected publication 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Subject terms

Heterogeneous catalysis
Chemical physics
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcreplying to S. Brandani et al. Nature Communications 10.1038/s41467-024-49821-w (2024)

Here we provide a comprehensive response to the Matter Arising by Brandani et al.1. They pointed out that the experimental data presented in Fig. 5e of our published article2 can be well-fitted by a first-order kinetic model (e.g., surface barriers3). Practically, the dual resistance model (DRM) can be used to fit the uptake process, reflecting the contributions of intracrystalline diffusion resistance and surface barriers with explicit physics meaning. In this reply, we illustrate the limitation of the dual-resistance model to decouple the surface barriers and intracrystalline diffusion based on the theorical analysis.

Introduction

Many scientists have previously proposed some interesting diffusion mechanisms such as “floating molecule”4, “levitation effect”5, and so on. However, little work has been done to consider the influence of molecular degrees of freedom on diffusion in confined channels, especially for long-chain molecules. In our work, we proposed a scheme to control the pore size of zeolite channels to adjust the degree of freedom of molecules to achieve ultrafast diffusion, in which we adopted both the molecular dynamics (MD) simulations (additional computational details in Supplementary Note 2 of Supplementary Information) and diffusion experiments as the research methods.

Results and discussion

In our original work2, we measured uptake curves of n-C12 and n-C4 over different zeolites and used DRM to obtain the surface permeability and intracrystalline diffusivity. In this reply, we would like to illustrate the limitation of DRM. In Fig. 1a, b, we plotted the theoretical uptake curves with different values of L (L = αl/D, where α is the surface permeability, D is the intracrystalline diffusivity and l is the half-length of the crystal). In Fig. 1a, the value of L is regulated by the different value of D while keeping the same value of α. In Fig. 1a, the uptake rate decreases as the increase in L value (L = αl/D), which is resulted from the decrease in intracrystalline diffusivity. In Fig. 1b, the uptake rate increases as the increase in L value, which is resulted from the increase in surface permeability. In Fig. 1a, for the range of L from 10−4 to 10−2, the uptake curves are unaffected by the intracrystalline diffusivity. Therefore, it can be concluded that when the value of L is below 10−1, the significant dominance of surface barriers on mass transfer makes it difficult to decouple the intracrystalline diffusivity from uptake curves by use of the dual-resistance model. In Fig. 1b, the value of L is regulated by the different value of α while keeping the same value of D. We show when the value of L is larger than 200 (for the sample interval ~1 s), it is difficult to determine the surface permeability by fitting the initial uptake curves.Fig. 1 Analysis based on dual-resistance model to determine the application range.

a Initial uptake rates by dual-resistance model for L = 1 × 10−4, 1 × 10−3, 1 × 10−2, and 1 × 10−1. b Initial uptake rates by dual-resistance model for L = 300, 240, and 120. The solid line is fitted by the equation (S1), and the discrete point is obtained by the equation (S2). Source data are provided as a Source Data file.

In our published paper2, we measured the diffusivities of n-C12 and n-C4 over TON, MTW, and AFI zeolites with small crystal sizes (≤ 2 μm). The results indicate that, in addition to intracrystalline diffusion, the presence of surface barriers can be observed. Therefore, we employed DRM to fit uptake curves, as DRM can be employed to decouple the intracrystalline diffusivity and surface permeability for a variety of guest molecules in different zeolites6–9. Remi et al.6 found the mass transport of methanol and butanol over individual H-SAPO-34 can be surface barriers-controlled based on the interference microscopy results in which they decoupled the intracrystalline diffusivity and surface permeability by DRM. Discrepant infrared signals in n-C12 presented in Supplementary Fig. 14 of our published paper2 show the different confinement effect imposed by intracrystalline-frameworks of TON, MTW and AFI zeolites. Therefore, we employed DRM to fit the intracrystalline diffusivity D and surface permeability α. Brandani et al. found it is not appropriate to use the value of intracrystalline diffusivity to present the experimental results due to the strong effect of surface barriers due to the small value of L. In this reply, based on our published data, we further presented the surface permeability in Fig. 2a, c. In Fig. 2b, We also calculated the inverse characteristic mass transport time (fitting by the first-order exponential model6 in Eq. 3) of n-C12 over TON-S, MTW-S and AFI-S samples is 1.18 × 10−2, 1.90 × 10−3, and 1.03 × 10−3  s−1, respectively. In Fig. 2d, the inverse characteristic mass transport time of n-C4 over TON-S, MTW-S and AFI-S samples is 4.60 × 10−5, 3.00 × 10−2 and 5.07 × 10−2  s−1, respectively. These results can obtain the same trend as found in MD simulations. In order to ensure the reliable fitting by DRM, zeolites with large MTW and AFI crystals (>100 μm) are adopted for experiments. In Fig. 3a, b, the uptake data over zeolites with large crystal sizes can be well-fitted by DRM, which shows the applicability of DRM for  the L range between 0.1-100. The effect of signal noise on the fitting errors, intracrystalline diffusivity and surface permeability has been examined. We found that the signal noise has significant effect on the fitting results (L value, intracrystalline diffusivity and surface permeability). In Fig. 3c, d, the average value and standard deviation of intracrystalline diffusivity and surface permeability are shown. We noticed that the uptake curves measured by microimaging techniques usually accompany with signal noise6,10. Therefore, how to determine the representative results under the signal noise effect is well-worth as an independent work. The average value of intracrystalline diffusivity and surface permeability can validate the trends obtained from MD simulations.Fig. 2 Uptake results of n-C12 and n-C4 in zeolites with small crystal size.

a The surface permeability and b inverse characteristic mass transport time (fitting by the first-order exponential model) of n-C12 over TON-S, MTW-S, and AFI-S samples at 298 K2. c The surface permeability and d, inverse characteristic diffusion mass transport time (fitting by the first-order exponential model) of n-C4 over TON-S, MTW-S, and AFI-S samples at 298 K2. The infrared microscopy (IRM) experimental conditions are as follows: the flowrate is 35 ml/min and the loading of zeolites in cell is ~1 mg. The intelligent gravimetric analyzer (IGA) experimental conditions are as follows: the loading of zeolites in IGA sample cell is 12 mg and the pressure of n-C4 changes from 0 to 4 mbar. Figure 2 were refitted from ref. 2. The value of bar graphs in (a, c) is obtained by fitting experimental data in ref. 2 by equation (S1). The value of bar graphs in (b, d) is obtained by fitting experimental data in ref. 2 by equation (S3). The error band is the standard error of fitting by equation (S1) and (S3). Source data are provided as a Source Data file.

Fig. 3 Uptake results of n-C12 and n-C4 in zeolites with large crystal size.

The uptake curves and fitting results of a n-C12 and b n-C4 over AFI-L and MTW-L samples. The infrared microscopy (IRM) experimental conditions are as follows: the flowrate is 35 ml/min and the loading of zeolites in cell is ~1 mg. c The intracrystalline diffusivity of n-C12 (Ln-C12,MTW = 0.77 ± 0.46 and Ln-C12,AFI = 1.45 ± 0.50) and n-C4 (Ln-C4,MTW = 3.44 ± 2.56 and Ln-C4,AFI = 1.61 ± 1.22) in AFI-L and MTW-L samples. d The surface permeability of n-C12 and n-C4 in AFI- L and MTW-L samples. a, b The solid line is fitted by the equation (S2), and discrete point is measured by the experiments. The error band is the standard error of experimental results of IGA measurements. Source data are provided as a Source Data file.

In the practical diffusion measurements over zeolite materials, the measured effective diffusivity can be affected by many factors besides zeolite frameworks, e.g., surface defects3, internal interfaces, and intergrowth structures11,12. At present, researchers usually compare the trends of experimental results and MD simulations to obtain a systematic understanding at both macro and micro levels13. So far, quantitative comparison between MD simulations and experimental diffusion measurements is a non-trivial task and remains a big challenge.

Methods

See Supplementary Information for details.

Supplementary information

Supplementary Information

Source data

source data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-49822-9.

Acknowledgements

This work was supported by the National Key R&D Program of China (No. 2022YFE0116000 (Z. Liu), and the National Natural Science Foundation of China (No. 22241801(Z. Liu), 22293021 (M. Ye), 22208337 (M. Gao)).

Author contributions

M.G. performed diffusion experiments and provided text for the corresponding section. J.Y. performed theoretical calculations and provided text for the corresponding section. Z.L., M.Y., and A.Z. supervised the work and revised the article.

Data availability

Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Mingbin Gao, Jiamin Yuan.

Change history

9/2/2024

A Correction to this paper has been published: 10.1038/s41467-024-52103-0
==== Refs
References

1. Brandani, S. et al. Diffusion anomaly in nanopores as a rich fieldfor theorists and a challenge for experimentalists. Nat. Commun.10.1038/s41467-024-49821-w (2024).
2. Yuan J Hyperloop-like diffusion of long-chain molecules under confinement Nat. Commun. 2023 14 1735 10.1038/s41467-023-37455-3 36977714
Yuan, J. et al. Hyperloop-like diffusion of long-chain molecules under confinement. Nat. Commun. 14, 1735 (2023).36977714 10.1038/s41467-023-37455-3
3. Kärger J Microimaging of transient guest profiles to monitor mass transfer in nanoporous materials Nat. Mater. 2014 13 333 343 10.1038/nmat3917 24651427
Kärger, J. et al. Microimaging of transient guest profiles to monitor mass transfer in nanoporous materials. Nat. Mater. 13, 333–343 (2014).24651427 10.1038/nmat3917
4. Derouane EG Andre JM Lucas AA Surface curvature effects in physisorption and catalysis by microporous solids and molecular sieves J. Catal. 1988 110 58 73 10.1016/0021-9517(88)90297-7
Derouane, E. G., Andre, J. M. & Lucas, A. A. Surface curvature effects in physisorption and catalysis by microporous solids and molecular sieves. J. Catal. 110, 58–73 (1988).10.1016/0021-9517(88)90297-7
5. Yashonath S Santikary P Diffusion of sorbates in zeolites Y and A: novel dependence on sorbate size and strength of sorbate-zeolite interaction J. Phys. Chem. 1994 98 6368 6376 10.1021/j100076a022
Yashonath, S. & Santikary, P. Diffusion of sorbates in zeolites Y and A: novel dependence on sorbate size and strength of sorbate-zeolite interaction. J. Phys. Chem. 98, 6368–6376 (1994).10.1021/j100076a022
6. Remi JCS The role of crystal diversity in understanding mass transfer in nanoporous materials Nat. Mater. 2016 15 401 406 10.1038/nmat4510 26689140
Remi, J. C. S. et al. The role of crystal diversity in understanding mass transfer in nanoporous materials. Nat. Mater. 15, 401–406 (2016).26689140 10.1038/nmat4510
7. Zhang B Boosting propane dehydrogenation over PtZn encapsulated in an epitaxial high-crystallized zeolite with a low surface barrier ACS Catal. 2022 12 1310 1314 10.1021/acscatal.1c04092
Zhang, B. et al. Boosting propane dehydrogenation over PtZn encapsulated in an epitaxial high-crystallized zeolite with a low surface barrier. ACS Catal. 12, 1310–1314 (2022).10.1021/acscatal.1c04092
8. Liu X A diffusion anisotropy descriptor links morphology effects of H-ZSM-5 zeolites to their catalytic cracking performance Commun. Chem. 2021 4 107 10.1038/s42004-021-00543-w 36697565
Liu, X. et al. A diffusion anisotropy descriptor links morphology effects of H-ZSM-5 zeolites to their catalytic cracking performance. Commun. Chem. 4, 107 (2021).36697565 10.1038/s42004-021-00543-w
9. Le TT Elemental zoning enhances mass transport in zeolite catalysts for methanol to hydrocarbons Nat. Catal. 2023 6 254 265 10.1038/s41929-023-00927-2
Le, T. T. et al. Elemental zoning enhances mass transport in zeolite catalysts for methanol to hydrocarbons. Nat. Catal. 6, 254–265 (2023).10.1038/s41929-023-00927-2
10. Smirnova O Precise control over gas-transporting channels in zeolitic imidazolate framework glasses Nat. Mater. 2024 23 262 270 10.1038/s41563-023-01738-3 38123813
Smirnova, O. et al. Precise control over gas-transporting channels in zeolitic imidazolate framework glasses. Nat. Mater. 23, 262–270 (2024).38123813 10.1038/s41563-023-01738-3
11. Karwacki L Morphology-dependent zeolite intergrowth structures leading to distinct internal and outer-surface molecular diffusion barriers Nat. Mater. 2009 8 959 965 10.1038/nmat2530 19767739
Karwacki, L. et al. Morphology-dependent zeolite intergrowth structures leading to distinct internal and outer-surface molecular diffusion barriers. Nat. Mater. 8, 959–965 (2009).19767739 10.1038/nmat2530
12. Liu L Imaging defects and their evolution in a metal–organic framework at sub-unit-cell resolution Nat. Chem. 2019 11 622 628 10.1038/s41557-019-0263-4 31086300
Liu, L. et al. Imaging defects and their evolution in a metal–organic framework at sub-unit-cell resolution. Nat. Chem. 11, 622–628 (2019).31086300 10.1038/s41557-019-0263-4
13. Nagumo R Takaba H Nakao S Prediction of ideal permeability of hydrocarbons through an MFI-type zeolite membrane by a combined method using molecular simulation techniques and permeation theory J. Phys. Chem. B 2003 107 14422 14428 10.1021/jp034401n
Nagumo, R., Takaba, H. & Nakao, S. Prediction of ideal permeability of hydrocarbons through an MFI-type zeolite membrane by a combined method using molecular simulation techniques and permeation theory. J. Phys. Chem. B 107, 14422–14428 (2003).10.1021/jp034401n
