Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

One-dimensional intergrowths in two-dimensional zeolite nanosheets and their effect on ultra-selective transport


Zeolite MFI is a widely used catalyst and adsorbent that also holds promise as a thin-film membrane. The discovery of nanometre-thick two-dimensional (2D) MFI nanosheets has facilitated methods for thin-film zeolite fabrication that open new horizons for membrane science and engineering. However, the crystal structure of 2D-MFI nanosheets and their relationship to separation performance remain elusive. Using transmission electron microscopy, we find that one- to few-unit-cell-wide intergrowths of zeolite MEL exist within 2D-MFI. We identify the planar distribution of these 1D or near-1D-MEL domains, and show that a fraction of nanosheets have high (~25% by volume) MEL content while the majority of nanosheets are MEL-free. Atomistic simulations show that commensurate knitting of 1D-MEL within 2D-MFI creates more rigid and highly selective pores compared to pristine MFI nanosheets, and permeation experiments show a separation factor of 60 using an industrially relevant (undiluted 1 bar xylene mixture) feed. Confined growth in graphite is shown to increase the MEL content in MFI nanosheets. Our observation of these intergrowths suggests strategies for the development of ultra-selective zeolite membranes.

This is a preview of subscription content

Access options

Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: Structure of the nanosheet.
Fig. 2: Identification of MFI and MEL heterostructures in nanosheets.
Fig. 3: Electron diffraction analysis for quantification of MEL in MFI nanosheets.
Fig. 4: Effect of MEL domains on separation and mechanical properties of the MFI nanosheet.

Data availability

The main data supporting the findings of this study are available within this Article and its Supplementary Information. Further data are available from the corresponding authors upon request.

Code availability

The template matching code used to identify MFI and MEL is available from K.A.M. and P.K. upon reasonable request.


  1. Rangnekar, N., Mittal, N., Elyassi, B., Caro, J. & Tsapatsis, M. Zeolite membranes—a review and comparison with MOFs. Chem. Soc. Rev. 44, 7128–7154 (2015).

    Article  CAS  Google Scholar 

  2. Daramola, M. O., Aransiola, E. F. & Ojumu, T. V. Potential applications of zeolite membranes in reaction coupling separation processes. Materials (Basel) 5, 2101–2136 (2012).

    Article  CAS  Google Scholar 

  3. Kokotailo, G., Lawton, S. & Olson, D. Structure of synthetic zeolite ZSM-5. Nature 272, 437–438 (1978).

    Article  CAS  Google Scholar 

  4. Flanigen, E. M. et al. Silicalite, a new hydrophobic crystalline silica molecular sieve. Nature 271, 512–516 (1978).

    Article  CAS  Google Scholar 

  5. Degnan, T. F., Chitnis, G. K. & Schipper, P. H. History of ZSM-5 fluid catalytic cracking additive development at mobil. Microporous Mesoporous Mater. 35–36, 245–252 (2000).

    Article  Google Scholar 

  6. Clough, M. et al. Nanoporous materials forge a path forward to enable sustainable growth: technology advancements in fluid catalytic cracking. Microporous Mesoporous Mater. 254, 45–58 (2017).

    Article  CAS  Google Scholar 

  7. Lai, Z. et al. Microstructural optimization of a zeolite membrane for organic vapor separation. Science 300, 456–460 (2003).

    CAS  Google Scholar 

  8. Hedlund, J. et al. High-flux MFI membranes. Microporous Mesoporous Mater. 52, 179–189 (2002).

    Article  CAS  Google Scholar 

  9. Jeon, M. Y. et al. Ultra-selective high-flux membranes from directly synthesized zeolite nanosheets. Nature 543, 690–694 (2017).

    Article  CAS  Google Scholar 

  10. Choi, M. et al. Stable single-unit-cell nanosheets of zeolite MFI as active and long-lived catalysts. Nature 461, 246–249 (2009).

    Article  CAS  Google Scholar 

  11. Varoon, K. et al. Dispersible exfoliated zeolite nanosheets and their application as a selective membrane. Science 334, 72–75 (2011).

    Article  CAS  Google Scholar 

  12. Kokotailo, G., Chu, P. & Lawton, S. Synthesis and structure of synthetic zeolite ZSM-11. Nature 275, 119–120 (1978).

    Article  CAS  Google Scholar 

  13. Ohsuna, T. & Terasaki, O. Electron microscopic study of intergrowth of MFI and MEL: crystal faults in B-MEL. J. Phys. Chem. B 101, 9881–9885 (1997).

    Article  CAS  Google Scholar 

  14. Pan, M. High resolution electron microscopy of zeolites. Micron 27, 219–238 (1996).

    Article  CAS  Google Scholar 

  15. Millward, G. R., Ramdas, S., Thomas, J. M. & Barlow, M. T. Evidence for semi-regularly ordered sequences of mirror and inversion symmetry planes in ZSM-5/ZSM-11 shape-selective zeolitic catalysts. J. Chem. Soc. Faraday Trans. 2 79, 1075–1082 (1983).

    Article  CAS  Google Scholar 

  16. Zhang, H. et al. Open-pore two-dimensional MFI zeolite nanosheets for the fabrication of hydrocarbon-isomer-selective membranes on porous polymer supports. Angew. Chem. Int. Ed. 55, 7184–7187 (2016).

    Article  CAS  Google Scholar 

  17. Aramburo, L. R. et al. The porosity, acidity, and reactivity of dealuminated zeolite ZSM-5 at the single particle level: the influence of the zeolite architecture. Chemistry 17, 13773–13781 (2011).

    Article  CAS  Google Scholar 

  18. Anderson, M. W. & Klinowski, J. Direct observation of shape selectivity in zeolite ZSM-5 by magic-angle-spinning NMR. Nature 339, 200–203 (1989).

    Article  CAS  Google Scholar 

  19. John, N. S., Stevens, S. M., Terasaki, O. & Anderson, M. W. Evolution of surface morphology with introduction of stacking faults in zeolites. Chemistry 16, 2220–2230 (2010).

    Article  CAS  Google Scholar 

  20. Ugurlu, O. et al. Radiolysis to knock-on damage transition in zeolites under electron beam irradiation. Phys. Rev. B 83, 113408 (2011).

    Article  CAS  Google Scholar 

  21. Mkhoyan, K. A., Maccagnano-Zacher, S. E., Kirkland, E. J. & Silcox, J. Effects of amorphous layers on ADF-STEM imaging. Ultramicroscopy 108, 791–803 (2008).

    Article  CAS  Google Scholar 

  22. Szostak, R., Pan, M. & Lillerud, K. P. High-resolution TEM imaging of extreme faulting in natural zeolite tschernichite. J. Phys. Chem. 99, 2104–2109 (1995).

    Article  CAS  Google Scholar 

  23. Terasaki, O. & Ohsuna, T. What can we observe in zeolite related materials by HRTEM? Catal. Today 23, 201–218 (1995).

    Article  CAS  Google Scholar 

  24. Pan, M. & Crozier, P. A. Quantitative imaging and diffraction of zeolites using a slow-scan CCD camera. Ultramicroscopy 52, 487–498 (1993).

    Article  CAS  Google Scholar 

  25. Bursill, L. A., Thomas, J. M. & Rao, K. J. Stability of zeolites under electron irradiation and imaging of heavy cations in silicates. Nature 289, 157–158 (1981).

    Article  CAS  Google Scholar 

  26. Treacy, M. M. J. & Newsam, J. M. Electron beam sensitivity of zeolite L. Ultramicroscopy 23, 411–419 (1987).

    Article  CAS  Google Scholar 

  27. Yoshida, K. & Sasaki, Y. Optimal accelerating voltage for HRTEM imaging of zeolite. J. Electron Microsc. 62, 369–375 (2013).

    CAS  Google Scholar 

  28. Kirkland, E. J. Advanced Computing in Electron Microscopy 2nd edn (Springer, 2010).

  29. Treacy, M. M. J. Z dependence of electron scattering by single atoms into annular dark-field detectors. Microsc. Microanal. 17, 847–858 (2011).

    Article  CAS  Google Scholar 

  30. Gonzalez, R. C. & Woods, R. E. Digital Image Processing 3rd edn (Pearson, 2008).

  31. Senftle, T. P. et al. The ReaxFF reactive force-field: development, applications and future directions. npj Comput. Mater. 2, 15011 (2016).

    Article  CAS  Google Scholar 

  32. Agrawal, K. V. et al. Oriented MFI membranes by gel-less secondary growth of sub-100 nm MFI-nanosheet seed layers. Adv. Mater. 27, 3243–3249 (2015).

    Article  CAS  Google Scholar 

  33. Yuan, W., Lin, Y. S. & Yang, W. Molecular sieving MFI-type zeolite membranes for pervaporation separation of xylene isomers. J. Am. Chem. Soc. 126, 4776–4777 (2004).

    Article  CAS  Google Scholar 

  34. Zhu, X. et al. Establishing hierarchy: the chain of events leading to the formation of silicalite-1 nanosheets. Chem. Sci. 7, 6506–6513 (2016).

    Article  CAS  Google Scholar 

  35. Rohling, R. Y., Szyja, B. M. & Hensen, E. J. M. Insight into the formation of nanostructured MFI sheets and MEL needles driven by molecular recognition. J. Phys. Chem. C 123, 5326–5335 (2019).

    Article  CAS  Google Scholar 

  36. Rangnekar, N. et al. 2D zeolite coatings: Langmuir–Schaefer deposition of 3 nm thick MFI zeolite nanosheets. Angew. Chem. Int. Ed. 54, 6571–6575 (2015).

    Article  CAS  Google Scholar 

  37. Kilaas, R. Optimal and near-optimal filters in high-resolution electron microscopy. J. Microsc. 190, 45–51 (1998).

    Article  Google Scholar 

  38. Frank, J. in Computer Processing of Electron Microscope Images (ed. Hawkes, P. W.) 187–222 (Springer, 1980).

  39. Cowley, J. M. & Moodie, A. F. The scattering of electrons by atoms and crystals. I. A new theoretical approach. Acta Crystallogr. 10, 609–619 (1957).

    Article  CAS  Google Scholar 

  40. Treacy, M. M. J., Newsam, J. M. & Deem, M. W. A general recursion method for calculating diffracted intensities from crystals containing planar stacking faults. Proc. R. Soc. Lond. A 433, 499–520 (1991).

    Article  Google Scholar 

  41. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).

    Article  CAS  Google Scholar 

  42. Newsome, D. A., Sengupta, D., Foroutan, H., Russo, M. F. & van Duin, A. C. T. Oxidation of silicon carbide by O2 and H2O: a ReaxFF reactive molecular dynamics study, part I. J. Phys. Chem. C 116, 16111–16121 (2012).

    Article  CAS  Google Scholar 

  43. Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558–561 (1993).

    Article  CAS  Google Scholar 

  44. Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996).

    Article  CAS  Google Scholar 

  45. Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994).

    Article  Google Scholar 

  46. Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758–1775 (1999).

    Google Scholar 

  47. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).

    CAS  Google Scholar 

  48. Grossfield, A. WHAM: the weighted histogram analysis method, version 2.0.9 (2013).

  49. Hutter, J., Iannuzzi, M., Schiffmann, F. & Vandevondele, J. cp2k: atomistic simulations of condensed matter systems. WIREs Comput. Mol. Sci. 4, 15–25 (2014).

    Article  CAS  Google Scholar 

  50. Hartwigsen, C., Goedecker, S. & Hutter, J. Relativistic separable dual-space Gaussian pseudopotentials from H to Rn. Phys. Rev. B 58, 3641–3662 (1998).

    Article  CAS  Google Scholar 

  51. VandeVondele, J. & Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 127, 114105 (2007).

    Article  CAS  Google Scholar 

  52. Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104 (2010).

    Article  CAS  Google Scholar 

  53. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 81, 511–519 (1984).

    Article  Google Scholar 

  54. Hoover, W. G. Canonical dynamics: equilibrium phase-space distributions. Phys. Rev. A 31, 1695–1697 (1985).

    Article  CAS  Google Scholar 

  55. Martyna, G. J., Klein, M. L. & Tuckerman, M. Nosé–Hoover chains: the canonical ensemble via continuous dynamics. J. Chem. Phys. 97, 2635–2643 (1992).

    Article  Google Scholar 

Download references


This work was primarily supported by the National Science Foundation (CBET-1705687). XRD was performed at Argonne National Laboratory supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, under contract no. DE-AC02-06CH11357. Parts of this work were carried out in the Characterization Facility, University of Minnesota, which receives partial support from the NSF through the MRSEC and NNIN programs (DMR-1420013). The FPMD simulations used resources of the Argonne Leadership Computing Facility, which is a DOE Office of Science User Facility supported under contract no. DE-AC02-06CH11357. Additional computer resources were provided by the Minnesota Supercomputing Institute. P.K. and E.O.F. acknowledge support from a Doctoral Dissertation Fellowship received from the Graduate School at the University of Minnesota. H.X. and T.D. acknowledge support from NSF 1332228. Q.X. acknowledges support from the National Natural Science Foundation of China (21471131).

Author information

Authors and Affiliations



P.K., along with K.A.M. and M.T., conceived the project. P.K. performed TEM data collection and analysis and developed the template matching algorithm with input from M.T. and K.A.M. H.Z., Q.X. and N.R. synthesized the 2D nanosheets. D.W.K. prepared the MEL-containing directly synthesized nanosheets (in the presence of graphite) with input from M.T., and S.G. performed TEM characterization. N.R. and M.S. performed XRD and AFM measurements. N.R. synthesized membranes and, with B.M., performed xylene membrane permeation measurements. H.X. and T.D. performed DFT simulations and, along with K.A.M. and P.K., interpreted the results on the mechanical properties of the nanosheets. E.O.F. and J.I.S. performed and analysed the FPMD simulations for the separation properties of the nanosheets. P.K., M.T. and K.A.M. prepared the manuscript with written contributions from all co-authors. K.A.M. and M.T. co-directed all aspects of the project.

Corresponding authors

Correspondence to Prashant Kumar, Michael Tsapatsis or K. Andre Mkhoyan.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

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

Extended data

Extended Data Fig. 1 MEL content quantification through XRD.

a, Indexed in-plane XRD data after background removal is labeled with corresponding peak positions. b, Bulk powder XRD patterns for MFI with 0–10% MEL are simulated using DiffaX computer code. The (102) peak decreases in intensity with increasing MEL content (bottom panel) while the (101) peak remains unchanged. c, Calibration curve from simulated data is plotted for % MEL vs I102/I101. Using the fitted curve, the MEL content of the entire sample is calculated: the experimentally determined I102/I101 ratio is 0.07, which corresponds to 4.3% MEL content.

Extended Data Fig. 2 ADF-STEM imaging of MFI nanosheets.

a, (left panel) Low magnification ADF-STEM image of an MFI nanosheet deposited on a carbon support. (middle and right panel) As-acquired raw ADF-STEM images of the sections shown in main text (Fig. 2a,b). b, ADF-STEM image (top) and the corresponding FFT (bottom) at increasing electron doses (left to right). Image acquired at the optimal dose is highlighted in blue. At less than 1.6x109 e, pores are barely visible (missing FFT spots as compared to FFT on right), while at higher doses the nanosheet is amorphized as seen by the disappeared spots in FFT (when compared to FFT on left). Cumulative electrons (e-), electron dose (e-2) and the total time of exposure (t) for image acquisition are listed at the top right corner of the as-acquired images. Maximum detectable periodicity is indicated by the arrows in FFT. c, Filtering steps (Left to right) for as-acquired raw image are shown. The same filtering steps are implemented as step 1 in template matching algorithm described in Extended data Fig. 4. Top row shows the ADF-STEM images. Middle row shows a smaller section from the image on top. Bottom column is the fast-fourier transform (FFT) of the image.

Extended Data Fig. 3 ADF-STEM defocus series of a 1.5 u.c. thick (along b-axis) MFI and MEL unit cell.

The simulation parameters were: Cs = 0 mm, Vo = 200 kV, Δf = -20 to 200 Å, 35 and 317 mrad for ADF detector inter- and outer-angles, 256 x 256 pix2 image size. The frozen phonons were not included in these calculations. Δf = 100 Å simulated image (enclosed by dotted lines) matched the experimental defocus conditions and is used as starting template (shown at the bottom after Gaussian blurring) for the template matching algorithm described in Extended data Fig. 4. Each simulated image is scaled in intensity individually from 0 to 255.

Extended Data Fig. 4 Algorithm for finding MFI and MEL distribution across ADF-STEM.

Demonstration of the cross-correlation based template matching code (written in MATLAB) to find multislice simulated MFI and MEL template in TEM images. a, Flowchart describing the four steps (bound by grey boxes) for executing the algorithm. b, Representative outputs of the code from the locations marked as 1* - 7* in the flowchart shown in a.

Extended Data Fig. 5 Results of template matching algorithm on different ADF-STEM images.

As-acquired raw images from the microscope are shown on the left. Dotted lines represent the section which is aligned with a- and c-axes of the crystal as vertical and horizontal lines respectively. Filtered images processed using radial Wiener, Band-pass, and moving 5x5 average filters of the section shown in the first column are shown in the second column. Results of the template matching algorithm are shown in the third and fourth column with shades of red and yellow circles representing two different orientations of MEL and MFI, respectively (same color scheme as in the Extended Data Fig. 4b and main text Fig. 2a). Third column images are obtained by cross-correlation of simulated unit-cell TEM images of MFI and MEL, with the filtered images in column 2, while the fourth column is obtained by cross-correlating the (cross-correlated and averaged) templates created from the matched units of each color in the third column with the corresponding images in column 1. Fourth column shows more complete filling since the scanned templates are created from the experimental image, thereby incorporating the noise levels and astigmatism present during collection of experimental data.

Extended Data Fig. 6 Effect of different stages of filtering an ADF-STEM image on template matching.

a, ADF-STEM at different stages of filtering. b, Effect of each stage of filtering on the ADF-STEM images. Difference image between Stage 2 and Stage 1 shows a low-frequency intensity variation arising due to the carbon residues on the nanosheet. Histogram of the image at Stage 3 shows a decrease in FWHM by 14% as compared to histogram of the image at Stage 2. c, Results of the template matching algorithm (algorithm is shown in Extended Data Fig. 4) at each stage of filtering show that the aforementioned reduction in FWHM correspond to a decrease in the noise in image.

Extended Data Fig. 7 Template matching results for BF-TEM images.

Filtered BF-TEM images (first column), template matching results using simulated MFI and MEL unit cell images (second column) and results using templates created from the matched units in run 1 (third column) for a, nominally pristine MFI nanosheet b, MFI-MEL nanosheet. Simulated and cross-correlated experimental units extracted from c, panels a for MFI nanosheet and d, panels b for MFI-MEL nanosheet. First row indicates the simulated BF-TEM templates used for run 1. Second and third rows indicate the cross-correlated and averaged units identified after run 1 and 2, respectively.

Extended Data Fig. 8 Effect of domain size on the reciprocal space of MFI nanosheet.

a, 3D model of a MFI nanosheet with lengths along a-, b-, and c-direction as Na, Nb, and Nc unit cells (u.c.). b, Changes in shape of the structure factor iso-surface at a reciprocal lattice point are shown, as the length of the MFI supercell is changed along a-direction from 1-4 u.c. The width of the reciprocal lattice point along a*- and b*- direction is inversely proportional to length along a- and b-direction (Na, Nb). c, MFI supercells (Na = 1-4 u.c., Nb = 1.5 u.c., Nc = 11 u.c.) along with different projections of the corresponding structure factor iso-surfaces in the 3D reciprocal space are shown. It is seen that as the length of the MFI supercell increases in a-direction, the elongation of the spot in a* - direction decreases. This corroborates with the streaking of spots along a* direction in acquired electron diffraction patterns.

Extended Data Fig. 9 Simulation of ED patterns for various nanosheet heterostructures.

a, ED patterns simulated using multislice code for nanosheet heterostructures shown in bottom right of each panel. Total length of the heterostructure along a-, b-, and c-direction is 11 u.c., 1.5 u.c. and 20 u.c. b, Diffraction pattern section indicated by dotted line in panel a is shown in a magnified view (top row) for models 1-10. Iso-surfaces for structure factor, |F|2 are plotted for the same models (bottom row). c, Line scans from |F|2 taken across (101) and (102) spot along a*-direction show changes in intensity of (102) spot while (101) spot remains unchanged as the MEL content and domain size varies. d, Intensity ratio vs % MEL content in the heterostructures (1 to 10) shows linear variation (line fit has R2=0.98). Intensities of individual spots are calculated by adding up the pixel values (or area under the 2D gaussian fit) for each spot after background subtraction. e, MFI domain width along a-direction (d) vs FWHM of (102) peak over FWHM of (101) peak (R) for heterostructures (1 to 10) shows exponential variation (fit has R2=0.97). FWHM is calculated along a*- and c*-direction by fitting a 2D-Gaussian function.

Extended Data Fig. 10 Analysis of experimental SAED patterns.

a, Representative SAED patterns showing variation in diffraction spots across different nanosheets. (101) spot is highlighted in red and (102) spot is highlighted in white. b, (101) and (102) spots from SAED patterns in panel a are magnified and shown. Streaking and splitting of (102) spots is seen while the corresponding (101) spot shows minimal variations in shape and intensity. c, Line-scans across diffraction spots shown in panel b are plotted. Linescans 2 and 3 show splitting of spots with distance along a* between the peaks being 0.0027 Å-1 and 0.0036 Å-1. This splitting indicates the presence of a repeating heterostructure in the nanosheet which has a periodicity of 37 nm and 28 nm respectively. d, Averaged (101) and (102) spots are shown from SAED patterns of 50 MFI nanosheets and 50 MFI-MEL nanosheets. e, Intensities of diffraction spots shown in panel d are listed in the table, with MFI nanosheets having 0% MEL and MFI-MEL nanosheets having 25% MEL. Percentage of MEL is calculated using the calibration chart (I102/I101 vs %MEL) shown in Extended data Fig. 9d.

Supplementary information

Supplementary Information

Membrane fabrication and testing discussion, Supplementary Figs. 1–7, Tables 1–3 and references 1–5.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kumar, P., Kim, D.W., Rangnekar, N. et al. One-dimensional intergrowths in two-dimensional zeolite nanosheets and their effect on ultra-selective transport. Nat. Mater. 19, 443–449 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI:

Further reading


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing