Targeted synthesis of zeolites from calculated interaction between zeolite structure and organic template

Abstract Zeolites, a class of silica-based porous materials, have been widely employed in the chemical industry for uses such as sorption, separation, catalysis and ion exchange. Normally, the synthesis of zeolites is performed in the presence of organic templates via a trial-and-error route, which is labor-intensive and empirical. In recent years, theoretical simulation from the interaction between a zeolite structure and an organic template has been used to guide the synthesis of zeolites, which is time-saving. In this review, recent progress in the targeted synthesis of zeolites from interaction between a zeolite structure and an organic template are briefly outlined including the design of new templates for zeolite synthesis, preparation of zeolites with new composition, development of novel routes for zeolite synthesis, synthesis of intergrowth zeolites, generation of novel zeolite structures, control of zeolite morphology and modulation of aluminum distribution in zeolites. These targeted syntheses reveal that the minimum energy principle from the theoretical simulation is key for guiding zeolite crystallization. This review will be important for zeolite researchers for rationally synthesizing zeolites and effectively designing new zeolite structures.


INTRODUCTION
Zeolites consist of corner-sharing TO 4 (Al, Si and P) tetrahedra as the primary building units, having well-defined micropores. In the initial stage, the artificial synthesis of zeolites such as A and X reported by Barrer and Milton was always carried out in inorganic systems. Later, the addition of organic quaternary ammonium cations and organic amines as organic templates in zeolite synthesis generates many novel zeolite structures. So far, >250 zeolite frameworks have been accepted by the Structure Commission of the International Zeolite Association and each of them has a special three-letter code [1]. Currently, zeolites have been widely used in many industrial processes such as sorption, separation, catalysis and ion-exchange processes [2][3][4][5][6][7][8][9][10].
Normally, the discovery of novel structure zeolites in the presence of organic templates is performed via the trial-and-error route, which is labor-intensive and empirical. To overcome this issue, theoretical simulation from the interaction between zeolite frameworks and organic templates to guide the experimental zeolite synthesis has been employed. For example, many organic templates have been theoretically designed for the rational synthesis of zeolites with new composition and structure, desirable synthesis of controllable zeolite morphology, construction of intergrowth zeolites and prediction of aluminum distribution in zeolite crystals.
In this review, we will briefly summarize recent progress in the targeted synthesis of zeolites from interaction between zeolite frameworks and organic templates. After a simple introduction, typical examples for the targeted synthesis of zeolites guided by theoretical simulation are shown. Finally, a conclusion and perspectives are given. We believe that this review will be helpful for zeolite researchers for rationally synthesizing zeolites and easily understanding zeolite crystallization.

TARGETED SYNTHESIS OF ZEOLITES
In the past decade, great efforts have been made in synthesizing zeolites guided by energies simulated from the interaction between zeolite frameworks and organic templates. The targeted synthesis of zeolites mainly includes the design of new templates for zeolite synthesis, preparation of zeolites with new compositions, development of novel routes for zeolite synthesis, synthesis of intergrowth zeolites, generation of novel zeolite structures, control of zeolite morphology and modulation of aluminum distribution in zeolites, as shown in Fig. 1.

Design of new organic templates for zeolite synthesis
Generally, the structure of organic templates is critical for the successful synthesis of specific zeolites. In many cases, organic templates with complex structures are costly in zeolite synthesis, which makes them very difficult to use for practical application in industrial processes. One of solutions is to discover low-cost organic templates to synthesize these zeolites assisted by theoretical simulation. This topic has been very successful in recent years and typical examples for the synthesis of zeolites such as BOG, STW, SFW, RTH and AEI structures are given in this subsection.
Boggsite as a natural zeolite is difficult to prepare artificially [11]. To solve this problem, Corma et al. calculated the stabilization energy between the BOG structure and organic templates, showing that the commercial tert-butyl-iminotris(dimethylamino)phosphorane has the best stabilization between this organic template and Boggsite. As a result, this organic template has been successfully directed for the synthesis of Boggsite zeolite (ITQ-47, BOG structure) from a high-throughput strategy [12]. This structure was confirmed using the X-ray diffraction technique, and 13 C and 31 P MAS NMR techniques showed that the organic molecules were indeed positioned in the zeolite pore after the crystallization. In addition, elemental analysis of the as-synthesized product displayed similar N/C and P/N ratios to those of pure organic molecules. These results demonstrate that the tert-butyl-iminotris(dimethylamino)phosphorane is an efficient organic template for directing the ITQ-47 zeolite.
HPM-1 zeolite with STW structure is the first pure silica chiral zeolite, which is potentially important for applications [13,14]. However, the synthesis of HPM-1 zeolite is very difficult, which was only performed in the presence of 2-ethyl-1,3,4trimethylimidazolium under a very narrow phase The stabilization energies (kJ (mol Si) -1 ) calculated between organic templates and a pure silica STW framework. Adapted with permission from ref. [15].
diagram [13]. Davis et al. calculated the stabilization interaction energies between a pure silica STW zeolite structure with a series of imidazole-based organic molecules, finding that pentamethylimidazolium has the same stabilization energy with the conventional organic template [15], as shown in Fig. 2. Notably, the pentamethylimidazolium has no rotational degrees of freedom compared with the conventional organic template that has many different conformations in aqueous solution. Therefore, pentamethylimidazolium might be more favorable to synthesize pure silica STW zeolite. As a result, it is directly synthesized pure silica STW zeolite in the presence of pentamethylimidazolium as an organic template. SSZ-52 zeolite with a SFW structure has good performance in the selective catalytic reduction of NO x with NH 3 (NH 3 -SCR) reaction [16,17]. However, the conventional template of polycyclic quarternary ammonium [N,N-diethyl-5,8-dimethylazonium bicyclo(3.2.2.)nonane] for the synthesis of SSZ-52 zeolite is very complex and costly, which strongly hinders its wide application. Davis et al. investigated the interaction energies between the SSZ-52 framework with organic templates that were inexpensive and easily synthesized from commercially available reagents. They found that N-ethyl-N- (2,4,4trimethylcyclopentyl)pyrrolidinium and N-ethyl-N-(3,3,5-trimethylcyclohexyl)pyrrolidinium with large stabilization energies (more negative) were suitable templates for the rational synthesis of SSZ-52 zeolite [18]. In order to synthesize SSZ-52 zeolite economically, the chosen organic templates should (i) be easily available raw materials, (ii) require no more than three steps for the preparation and (iii) have a high yield of organic molecules.
RTH zeolite exhibits excellent performance in methanol-to-olefin (MTO) and NH 3 -SCR reactions [19][20][21]. In the past decade, great progress has been made in developing different routes to synthesize RTH zeolite [19][20][21][22]. However, the crystallization rate of RTH zeolite is still slow, giving low synthesis efficiency. Xu et al. found that the interaction energy of the novel organic template 2,6methyl-N-methylpyridinium cation with the RTH zeolite framework is much larger (more negative) than those of conventional templates for the synthesis of RTH zeolites such as N-ethyl-N-methyl-5,7,7trimethyl-azonium bicyclo[4.1.1]octane cations and 1,2,3-trimethylimidazolium cations. These results mean that the 2,6-methyl-N-methylpyridinium cation as an organic template might be more appropriate for directing the RTH zeolite structure [23]. In the following experiments, they found that the activation energy of RTH zeolite directed by this novel organic template was much lower than that of the conventional one. As a result, compared with the conventional synthesis (72 h at 130 • C), the crystallization of RTH zeolite in this work took a very short time (50 min at 240 • C and 12 h at 130 • C). In addition, RTH zeolite exhibited high selectivity for propylene and ethylene of >75% in the MTO reaction, which would be of great importance for the selective production of light olefins in practical applications. Moreover, the Cu 2+ -exchanged RTH zeolite showed comparable performance to that of the commercial Cu-SSZ-13 zeolite, suggesting a promising zeolite for industrial applications in the NH 3 -SCR reaction.
SSZ-39 zeolite with an AEI structure is also a very important candidate for the NH 3 -SCR reaction [24][25][26]. However, Si/Al ratios in the conventional synthesis of AEI zeolite were <10, giving relatively low stability [24]. To enhance the stability, Schmidt et al. simulated the interaction energies between the zeolite framework with organic templates, predicting that a bulky organic template of N-ethyl-Nmethyl-2,2,6,6-tetramethylpiperidinium has larger (more negative) stabilization energy than that of conventional organic templates for the formation of a high-silica SSZ-39 zeolite. Accordingly, they have successfully used this template for the synthesis of SSZ-39 with a Si/Al ratio of 13. As expected, it significantly enhanced the stability of the SSZ-39 zeolite, thus facilitating its application [27]. Considering the complexity and high cost of alkyl-substituted piperidinium cations in zeolite synthesis, Schwalbe-Koda et al. performed theoretical simulation by phase competition and shape analysis, reporting that a facile tris(dimethylamino)(methyl)phosphonium might be a novel candidate for synthesizing SSZ-39, which was confirmed experimentally. The obtained zeolite products displayed similar crystallinity and textural properties to those prepared with the conventional organic template, but the crystallization conditions are much milder [28].
In addition, Corma et al. showed a new method for the synthesis of zeolites by the employment of organic templates that mimic the transition state of pre-established reactions such as toluene disproportionation, ethylbenzene isomerization and adamantane synthesis [29]. Theoretical calculation displayed that the above organic templates are very suitable for synthesizing zeolites (ITQ-27, ITQ-64 and MIT-1). Later, they reported that small-pore zeolites with eight-membered rings (AEI, RTH and CHA) could also be successfully synthesized by using organic templates that mimic the key molecular species in the process of the MTO reaction [30].

Preparation of zeolites with new composition
Although >250 zeolite frameworks have been discovered so far, there are only ∼20 zeolites being applied in industrial processes [31][32][33]. Among them, most of the zeolites are aluminosilicate compositions, which have both suitable acidic sites and excellent thermal and hydrothermal stabilities. Aluminogermanosilicate zeolites such as IWR and ITH, one of typical zeolite families, also exhibit excellent performance in the fields of sorption and catalysis [34][35][36][37][38], but the existence of germanium species in zeolite frameworks not only increase the zeolite cost, but also reduces their stabilities, which strongly limits their practical application [39,40]. Therefore, it is highly desirable to directly synthesize the aluminosilicate zeolites with the same structures as aluminogermanosilicate zeolites. Hong et al. reported that aluminosilicate IWR zeolite could be directly synthesized using a designed organic template from a computational simulation screening by a strong interaction between a zeolite framework and pyrrolidine-based cations [41]. Based on minimum energy principle, it is suggested that larger stabilization energies (more negative) mean higher suitability for the zeolite synthesis. Balancing the interaction stabilization energy and cost of organic templates, a suitable organic template has been selected for directing the synthesis of targeted IWR zeolite, as shown in Fig. 3. As expected, the aluminosilicate IWR zeolite with a fourcoordinated Al species (Al-IWR) has been successfully prepared from the designed organic template, in good agreement with the theoretical simulation. The Al-IWR zeolite could give an Si/Al ratio as low as 15. In methanol-to-propylene tests, the Al-IWR zeolite showed higher selectivity for propylene than that of commercial ZSM-5 zeolite, which might be of importance for the selective production of propylene from methanol-a platform compound from coal and biomass.
Later, Lei et al. reported the direct synthesis of aluminosilicate ITH zeolite by employing a designed cationic oligomer as an efficient organic template [42]. First, the theoretical simulation was performed to search for the suitable cationic oligomers. As observed from the interaction stabilization energy between the ITH framework and the structure of the cationic oligomer, it was suggested that the cationic oligomer was suitable for directing the crystallization of ITH zeolite. Second, it was found that the cationic oligomer had stronger complexation ability with aluminum species than that of the conventional organic template. Combining both advantages, aluminosilicate ITH zeolite was successfully synthesized. Moreover, aluminosilicate ITH zeolite exhibited excellent selectivity for propylene and a long lifetime in the MTP reaction. In addition, aluminosilicate ITH zeolite was a good additive for enhancing selectivities for butylene and propylene in the fluid catalytic cracking process. The employment of a cationic oligomer as an efficient organic template guided by theoretical simulation might open the door to directly preparing other aluminosilicate zeolites that could not be synthesized up to now.

Development of novel routes for zeolite synthesis
The modern synthesis of zeolites normally requires the presence of organic templates and sometimes costly raw materials are necessary, which are environmentally unfriendly and have low efficiency [43]. Thus, developing novel routes for zeolite synthesis in a sustainable manner is highly desirable. Despite great efforts having been made in past decades for synthesizing zeolites [44][45][46][47][48][49], there are still many challenges in the field of zeolites synthesis with environmentally unfriendly features. For example, pure silica zeolites cannot be synthesized in the absence of organic templates and the use of organic templates both increases the zeolite cost and produces harmful gases. The synthesis of SSZ-39 zeolite is always performed in the presence of high-silica Y zeolite as a raw material and the use of this kind of highsilica Y zeolite is costly. To address these issues, both theoretical simulations and experimental work have been performed, and the typical examples are the organotemplate-free synthesis of pure silica MFI zeolite and the transformation of MFI and Beta zeolites into SSZ-39 zeolite.
It is well known that pure silica zeolites have been widely used as catalytic supports and adsorbents in industrial processes [50][51][52], but their synthesis must use organic templates. Wu et al. for the first time showed that pure silica zeolites could be successfully synthesized in the presence of zeolite seeds and ethanol but in the absence of organic templates [53]. In this case, the zeolite seeds were used for directing zeolite growth, while the ethanol was employed for filling the zeolite pores. The filling role of ethanol is well confirmed by theoretical simulation, as shown in Fig. 4a. In the process of   SSZ-39 zeolite with an AEI structure is one of the most promising candidates instead of SSZ-13 zeolite used for the selective catalytic reduction of NO x with NH 3 (NH 3 -SCR). In general, SSZ-39 zeolite is synthesized by interzeolite transformation from highsilica Y zeolite, which is very costly [24]. Therefore, it is hopeful to use low-cost zeolites such as ZSM-5 and Beta as raw materials for the synthesis of SSZ-39 zeolite. However, it seems difficult because the interzeolite transformation from ZSM-5 and Beta zeolites into SSZ-39 zeolite would challenge the rule from a low framework density into a high framework density of zeolites. Notably, in published work, organic templates are always present for the interzeolite transformation but the contribution of organic templates for this transformation is completely ignored. If the interaction between organic templates with zeolite frameworks is considered, it is simulated that the stabilization energies between organic templates and zeolite structures (FAU, * BEA and MFI) are higher than that between organic templates and AEI zeolite structures, which suggests that it is possible to synthesize SSZ-39 zeolite by interzeolite transformations from ZSM-5 and Beta zeolite, as shown in Fig. 5 [54]. Later experiments showed that the interzeolite transformation of SSZ-39 zeolite from ZSM-5 and Beta zeolites was successful. After being exchanged with Cu 2+ , the SSZ-39 zeolite exhibited extraordinary performance in the NH 3 -SCR reaction, which was fully comparable with that of SSZ-39 zeolite exchanged with Cu 2+ from the conventional transformation synthesis.
Because raw materials such as colloidal silica are much cheaper than MFI and Beta zeolites, it is reasonable to directly synthesize SSZ-39 zeolite from  Figure 6. (a) The comparison between the binding energies of organic templates in AEI and CHA zeolites; (b) the relationship between the shapes of organic templates and the corresponding binding energies of AEI and CHA zeolites. Adapted with permission from ref. [28]. Structure C 3 C 4 C 5 C 6 C 3 C 4 C 5 C 6 C 3 C 4 C 5 C 6 C 3 C 4 C 5 C 6 C 3 C 4 C 5 C 6 C 3 C 4 C 5 C 6 C 3 C 4 C 5 C 6 C 3 C 4 C 5 C 6 colloidal silica. Initially, the as-synthesized products always form MOR zeolite rather than SSZ-39 zeolite. After the theoretical simulation of stabilization energies between organic templates and the zeolite structures (MOR and AEI), it is found that both stabilization energies are very low and similar, indicating that both zeolites could be formed for the synthesis of SSZ-39 zeolite. By adjusting the ratios of Si/Al in the starting materials for the avoidance of the MOR phase, it is successful for the direct synthesis of SSZ-39 zeolite using colloidal silica as a source [55], which further reduces the cost of SSZ-39 zeolite for practical applications.

Synthesizing intergrowth zeolites
Compared with a sole-zeolite structure, intergrowth zeolites such as CHA/AEI, MFI/MEL, FAU/EMT, MTT/TON and SBS/SBT always give distinctive properties leading to excellent performance for practical applications [56][57][58][59][60]. In general, the synthesis of intergrowth zeolites was performed in the presence of dual organic templates, which required fine adjustment of the synthetic parameters to avoid the formation of separated zeolite phases. Thus, it is expected to design a bi-selective organic template for directing the targeted intergrowth zeolites. Schwalbe-Koda et al. for the first time reported the formation of aluminosilicate CHA/AEI intergrowth zeolite from N-ethyl-N-isopropyl-Nmethylpropan-2-aminium as an organic template through theoretical calculations. This template has similar binding energies to both pure CHA and AEI phases as well as its shape positioned on the phase boundary between the CHA and AEI structures, as shown in Fig. 6. Based on theoretical calculation, the aluminosilicate CHA/AEI zeolite was successfully synthesized under suitable conditions. The intergrowth property of CHA/AEI zeolite was confirmed by the SEM, XRD and HRTEM techniques. In addition, the result of DIFFaX simulation confirmed the CHA/AEI intergrowth ratio at ∼1.0 [28]. Later, they employed computational simulations and data mining to design bi-selective organic templates used for directing formation of other intergrown zeolites (BEC/ISV, MTT/TON and CHA/AFX). Furthermore, a bi-selective organic template for the potential formation of a hypothesized intergrowth zeolite structure (AEI/SAV) was also proposed [61].

Generation of novel zeolite structures
The generation of novel zeolite structures is always a hot topic in the field of zeolite research [62][63][64][65][66][67][68][69][70]. Currently, most novel zeolite structures are made from organic templates via a trial-and-error route, which is time-consuming and has low efficiency. One efficient method is to use theoretical simulation.
A typical example reported by Hong et al. is the successful synthesis of a novel zeolite structure of PST-30 guided by calculating the stabilization energies between organic templates and hypothetical zeolite structures, as shown in Fig. 7 [71]. At first, the stabilization energies of the reported organic templates for PST-21 and PST-22 zeolites were calculated [72], showing that 1,2,3-trimethylimidazolium had a stabilization energy of -13.7 kJ (mol Si) -1 in the PST-21 zeolite structure, while 1,3,4-trimethylimidazolium and 1,2,3,4-tetramethylimidazolium gave the stabilization energies of approximately -13.4 and -13.8 kJ (mol Si) -1 in the PST-22 zeolite structure. Therefore, the benchmark energy for considering whether organic templates would direct the formation of novel targeted zeolite structures was set at -13.0 kJ (mol Si) -1 . After calculating the stabilization energies by interaction between pyrazoliumbased and diazole-based organic templates with hypothetical zeolite structures, it was found that [1,1 -(1,4-butane-diyl)bis(2,5-dimethyl-1H-pyrazol-2ium)] might be a suitable candidate for directing BRE L1 1 (PST-30) zeolite due to its large stabilization energy (more negative) by computational simulation. After the introduction of this template into the synthesis system, PST-30 zeolite with good crystallinity could be successfully synthesized under suitable crystallization conditions, as confirmed by the XRD technique, in good agreement with the simulated XRD pattern. The structure analysis of PST-30 zeolite showed that it had a 2D micropore consisting of 8-MRs (4.2 × 4.4Å) and 10-MRs (4.7 × 6.7Å). This special micropore system of PST-30 zeolite displayed excellent catalytic performance in the skeletal isomerization of the 1-butene to isobutene reaction. The targeted synthesis of the novel zeolite structure chosen from the database of hypothetical zeolites might open the door for synthesizing other novel zeolite structures guided by computational calculations of organic templates.
Recently, great efforts have been witnessed for this kind of work [83,84] and one line of progress has been to control zeolite morphology assisted with theoretical simulation.
Xu et al. reported that ultra-thin nanosheets of FER zeolite (N-FER) with a thickness of ∼8 nm were successfully synthesized by the employment of a small organic ammonium salt (N,N-diethylcis-2,6-dimethylpiperidinium, DMP) [85]. Theoretical simulations were used for understanding the N-FER zeolite with four-unit-cells along with the [100] direction. Figure 8a shows the calculated adsorption energy associated with the number of DMP molecules on FER [100] surfaces, showing that the model with four DMP molecules had the largest adsorption energy (the most negative). Then, the dependence of the surface energy of FER zeolite on layers of its structure was calculated, showing that three-to five-layered zeolite structure is stable, as shown in Fig. 8b. The above results indicated that the DMP molecules would adsorb on the [100] plane of FER zeolite, which limited the FER layers' growth and thus formed the ultra-thin nanosheets FER zeolite with a thickness of ∼8 nm. Furthermore, theoretical calculation also predicted the position of the DMP molecules in the FER zeolite structure, showing that two DMP molecules per unit cell had the largest stabilization energy (the most negative), as shown in Fig. 8c. This result suggests that the DMP is an efficient organic template for directing the formation of a FER zeolite structure. In addition, the N-FER zeolite exhibited higher conversion and isomer selectivity than the conventional FER zeolite due to the feature of ultra-thin nanosheets.
Later, Peng et al. showed that a small, cheap and simple imidazole-based molecule could not only direct the formation of FER zeolite structure, but also limit the growth of the FER zeolite nanosheet by π -π stacking interaction [86].  T1  T2  T3  T4   T1  T2  T3  T4   39  result, ultra-thin FER zeolite nanosheets (ECNU-17) with a thickness of only 2.9 nm could be successfully synthesized. Theoretical simulation showed that the ECNU-17 nanosheets formed easily when two layers of organic template molecules were positioned on the external surface of the layers. Due to a significant increase in accessible acidic sites, the ultra-thin FER zeolite exhibited superior catalytic performance in the catalytic cracking of 1,3,5-triisopropylbenzene, a model compound for bulky molecules.

Modulation of aluminum distribution in zeolites
Aluminosilicate zeolites as acidic catalysts have been widely applied for petrochemical processes. Notably, the same aluminosilicate zeolites sometimes show quite different performances, which are related to the distinguished aluminum distribution governed by the different synthetic parameters [87][88][89][90][91][92][93][94]. Among these parameters, the role of organic templates is critical for the modulation of aluminum distribution in zeolites, where the development of computational simulation techniques offers a good chance for the design of organic templates. Muraoka et al. reported the preparation of IFR zeolites with controllable aluminum distribution at different T sites using three different organic templates. Theoretical simulation showed that the T1 site was the most stable position for aluminum for any of the organic templates, which was in good agreement with the result from the MAS NMR technique. At the same time, short N + · · ·O(Al) distances would be observed for three organic templates and the T1 site, which confirms the highest occupancy of the T1 site in three products. On the contrary, theoretical simulation showed the highest energy at the T3 site confirmed by the long N + · · ·O(Al) distances, which led to the low occupancy in the T3 site of products, as shown in Fig. 9 [95].
Chu et al. demonstrated the constrained aluminum sites in aluminosilicate FER zeolite. Theoretical simulation was performed by interaction between the FER zeolite structure and protonated pyridines. The results showed that the T1 and T3 sites of FER zeolite were the preferable occupancy, in good agreement with the results of Rietveld refinement. As a result, the above FER zeolite displayed limited catalytic performance in the carbonylation of dimethyl ether since the acidic sites of T1 and T3 in the open channels cannot act as the active center for this reaction [96].
More recently, Schwalbe-Koda et al. showed that the aluminum distribution of SSZ-13 zeolite was associated with the spatial charge distribution of organic templates through theoretical simulation. The nitrogen center descriptor showed that the specific organic template with the N center was close to the center of the CHA cage, which could make it difficult to form more paired aluminum sites [28]. The above principle might be a good pathway to the modulation of aluminum distribution in SSZ-13 zeolite, then improving the performance of SSZ-13 zeolite in the NH 3 -SCR reaction.

CONCLUSION AND PERSPECTIVE
In summary, the targeted synthesis of zeolites from the calculated interaction between zeolite structures and organic templates has been briefly concluded. The targeted syntheses of zeolites mainly include the design of new templates for zeolite synthesis, preparation of zeolites with new composition, development of novel routes for zeolite synthesis, synthesis of intergrowth zeolites, generation of novel zeolite structures, control of zeolite morphology and modulation of aluminum distribution in the zeolites.
In these examples, it is revealed that the minimum energy principle for the calculated interaction energies between zeolite frameworks and organic templates is critical for the successful experimental synthesis of targeted zeolites. Of course, a minimum energy calculation is capturing the local minimum in the potential energy surface. Following this way, it is believable that more work is still on the way for the targeted synthesis of zeolites guided by theoretical simulation.
Despite great efforts having been made for the targeted synthesis of zeolites, there are still challenges to be addressed. For example, many zeolite structures show excellent catalytic performances but their syntheses are costly due to the use of complex organic templates; many zeolite structures are easily synthesized in the lab but they are difficult for scale-up preparation for industrial production. Possibly, theoretical simulation offers an alternative way to solve these kinds of problems in the future.
Considering the costly organic templates in many zeolite syntheses, it is strongly desirable to use inorganic templates such K + cation instead of organic templates for synthesizing zeolites, which can be guided from the calculated interaction energies between zeolite frameworks and inorganic cations. For example, Han et al. showed the successful preparation of a new high-silica KFI zeolite with a Si/Al ratio of >5.0 in the presence of K + cation but in the absence of any organic templates, which is effectively inspired by interaction energies between the zeolite structure and the K + cation, as shown in Fig. 10 [97]. The simulated interaction energies between zeolite frameworks and inorganic cations might help to experimentally synthesize more zeolite structures in the future.
However, the targeted synthesis of zeolite guided by computational simulation also has technical challenges and limitation. For example, the effect of synthetic conditions such as having water as a solvent are not often considered in computational simulations; the defects in targeted zeolite frameworks are usually ignored in theoretical calculations. In addition, it is worth noting that the zeolites are metastable in terms of thermodynamics, which are often influenced by many synthetic parameters [98,99].
Possibly, the targeted synthesis of zeolites guided by theoretical simulation would open a new door for synthesizing zeolites, which is should be efficient, low-cost and time-saving. Therefore, more efforts should be made for this kind of research in the future.