Could Mesophases Play a Role in the Nucleation and Polymorph Selection of Zeolites?

The Current Understanding of Mechanistic Pathways in Zeolite Crystallization

Zeolite catalysts and adsorbents have been an integral part of many commercial processes and are projected to play a significant role in emerging technologies to address the changing energy and environmental landscapes. The ability to rationally design zeolites with tailored properties relies on a fundamental understanding of crystallization pathways to strategically manipulate processes of nucleation and growth. The complexity of zeolite growth media engenders a diversity of crystallization mechanisms that can manifest at different synthesis stages. In this review, we discuss the current understanding of classical and nonclassical pathways associated with the formation of (alumino)silicate zeolites. We begin with a brief overview of zeolite history and seminal advancements, followed by a comprehensive discussion of different classes of zeolite precursors with respect to their methods of assembly and physicochemical properties. The following two sections provide detailed discussions of nucleation and growth pathways wherein we emphasize general trends and highlight specific observations for select zeolite framework types. We then close with conclusions and future outlook to summarize key hypotheses, current knowledge gaps, and potential opportunities to guide zeolite synthesis toward a more exact science.

Ranking the synthesizability of hypothetical zeolites with the sorting hat

Zeolites are nanoporous alumino-silicate frameworks widely used as catalysts and adsorbents. Even though millions of siliceous networks can be generated by computer-aided searches, no new hypothetical framework has yet been synthesized. The needle-in-a-haystack problem of finding promising candidates among large databases of predicted structures has intrigued materials scientists for decades; yet, most work to date on the zeolite problem has been limited to intuitive structural descriptors. Here, we tackle this problem through a rigorous data science scheme-the "Zeolite Sorting Hat"-that exploits interatomic correlations to discriminate between real and hypothetical zeolites and to partition real zeolites into compositional classes that guide synthetic strategies for a given hypothetical framework. We find that, regardless of the structural descriptor used by the Zeolite Sorting Hat, there remain hypothetical frameworks that are incorrectly classified as real ones, suggesting that they might be good candidates for synthesis. We seek to minimize the number of such misclassified frameworks by using as complete a structural descriptor as possible, thus focusing on truly viable synthetic targets, while discovering structural features that distinguish real and hypothetical frameworks as an output of the Zeolite Sorting Hat. Further ranking of the candidates can be achieved based on thermodynamic stability and/or their suitability for the desired applications. Based on this workflow, we propose three hypothetical frameworks differing in their molar volume range as the top targets for synthesis, each with a composition suggested by the Zeolite Sorting Hat. Finally, we analyze the behavior of the Zeolite Sorting Hat with a hierarchy of structural descriptors including intuitive descriptors reported in previous studies, finding that intuitive descriptors produce significantly more misclassified hypothetical frameworks, and that more rigorous interatomic correlations point to second-neighbor Si-O distances around 3.2-3.4 Å as the key discriminatory factor.

The Promoter Role of Amines in the Condensation of Silicic Acid: A First-Principles Investigation

Though well-recognized, the molecular-level understanding of the multifunctional roles of amines in the condensation of polysilicic acids, which is one of the key processes in hydrothermal synthesis of zeolites, is still limited. Taking ethylamine as a prototype, we investigated the mechanism of polysilicic acid condensation in the existence of organic amines in aqueous solution with extensive first-principles-based calculations. Because of the high proton affinity, ethylamine exists as amine silicates and alters the subsequent condensation mechanisms from a 1-step lateral attack mechanism accompanied with simultaneous intermolecular proton transfer in neutral aqueous solution to a 2-step S_N2-like mechanism. Specifically, the 5-coordinated Si species that were not observed on pathways of condensation in neutral solution are effectively stabilized by the ethylamine cations as intermediates, and the barriers for condensation of ortho-silicic acid are significantly reduced from 133 kJ/mol in neutral solution to 58 and 63 kJ/mol for formation of the 5-coordinated Si intermediate and proton transfer for water release, respectively. Similar variations of mechanisms and barriers for condensation were also observed in the formation of cyclic trimers as well as linear and cyclic tetramers of ortho-silicic acids. Based on these, it was proposed that apart from acting as structure-directing agents, pore fillers, and pH adjusters, organic amines can also function as promoters in the condensation of polysilicic acids.

Formation of porous ice frameworks at room temperature

Bulk crystalline ices with ultralow densities have been demonstrated to be thermodynamically metastable at negative pressures. However, the direct formation of these bulk porous ices from liquid water at negative pressures is extremely challenging. Inspired by approaches toward porous media based on host-guest chemistry, such as metal-organic frameworks and covalent organic frameworks, we herein demonstrate via molecular dynamics simulations that a class of ultralow-density porous ices with upright channels can be formed spontaneously from liquid water at 300 K with the assistance of carbon nanotube arrays. We refer to these porous ice polymorphs as water oxygen-vertex frameworks (WOFs). Notably, our simulations revealed that the liquid-WOF phase transition is first-order and occurs at room temperature. All the WOFs exhibited the unique structural feature that they can be regarded as assemblies of nanoribbons of hexagonal bilayer ice (2D ice I) at their armchair or zigzag edges. Based on density functional theory calculations, a comprehensive phase diagram of the WOFs was constructed considering both the thermodynamic and thermal stabilities of the porous ices at negative pressures. Like other types of porous media, these WOFs may be applicable to gas storage, purification, and separation. Moreover, these biocompatible porous ice networks may be exploited as medical-related carriers.

A dramatic conformational effect of multifunctional zwitterions on zeolite crystallization

Carnitine functions as a mesoporogen in LTA zeolite synthesis whereas its structural analogue acetylcarnitine acts as a crystal growth modifier. An array of experimental and theoretical studies reveal a remarkable effect of molecular conformation on the actual roles of organic functional groups during zeolite crystallization.

Oxidative desulfurization of model oil over Ta-Beta zeolite synthesized via structural reconstruction

As to metallosilicate zeolites, ions with larger size such as Ta⁵⁺ in the gels greatly retarded their crystallization during the hydrothermal synthesis, affording long-winded synthesis periods, up-limited framework-substituted metal contents, or even frustrated outcome. An efficient hydrothermal synthesis strategy for metallosilicate, in this case of Ta framework-substituted *BEA zeolite, via structural reconstruction was proposed to stride the gap. The Ta content in our developed Ta-Beta-Re-50 zeolite achieved up to 5.48 % (Si/Ta = 52), breaking through the limitation of Ta contents for conventional method (Si/Ta > 100). Additionally, this Ta-Beta-Re zeolite possessed nanosized crystals (20-40 nm) and short crystallization time (8 h), significantly improving space-time yields of practical zeolite production. Through spectroscopic study, it was confirmed that the existence of zeolite structural units intensively facilitated the formation of nucleation and crystal growth. This innovative Ta-Beta zeolite demonstrated high catalytic performances for oxidation desulfurization, far outperforming traditional fluoride-mediated Ta-Beta-F, which was ascribed to its excellent diffusion properties and incredible high isolated Ta contents. Additionally, the catalytic performance of Ta-Beta-Re could be regenerated after simple calcination and the deactivation may be caused by pore blocking of organics. This work provides a new method for rationally design and construction of metallosilicate materials with high activity for catalytic oxidation applications, which can bridge the conceptual and technical gap between periodic trends and zeolite material synthesis.

An ultralow-density porous ice with the largest internal cavity identified in the water phase diagram

The recent back-to-back findings of low-density porous ice XVI and XVII have rekindled the century-old field of the solid-state physics and chemistry of water. Experimentally, both ice XVI and XVII crystals can be produced by extracting guest atoms or molecules enclosed in the cavities of preformed ice clathrate hydrates. Herein, we examine more than 200 hypothetical low-density porous ices whose structures were generated according to a database of zeolite structures. Hitherto unreported porous EMT ice, named according to zeolite nomenclature, is identified to have an extremely low density of 0.5 g/cm3 and the largest internal cavity (7.88 Å in average radius). The EMT ice can be viewed as dumbbell-shaped motifs in a hexagonal close-packed structure. Our first-principles computations and molecular dynamics simulations confirm that the EMT ice is stable under negative pressures and exhibits higher thermal stability than other ultralow-density ices. If all cavities are fully occupied by hydrogen molecules, the EMT ice hydrate can easily outperform the record hydrogen storage capacity of 5.3 wt % achieved with sII hydrogen hydrate. Most importantly, in the reconstructed temperature-pressure (T-P) phase diagram of water, the EMT ice is located at deeply negative pressure regions below ice XVI and at higher temperature regions next to FAU. Last, the phonon spectra of empty-sII, FAU, EMT, and other zeolite-like ice structures are computed by using the dispersion corrected vdW-DF2 functional. Compared with those of ice XI (0.93 g/cm3), both the bending and stretching vibrational modes of the EMT ice are blue-shifted due to their weaker hydrogen bonds.

Following the nucleation pathway from disordered liquid to gyroid mesophase

Mesophases have order intermediate between liquids and crystals and arise in systems with frustration, such as surfactants, block copolymers, and Janus nanoparticles. The gyroid mesophase contains two interpenetrated, nonintersecting chiral networks that give it properties useful for photonics. It is challenging to nucleate a gyroid from the liquid. Elucidating the reaction coordinate for gyroid nucleation could assist in designing additives that facilitate the formation of the mesophase. However, the complexity of the gyroid structure and the extreme weakness of the first-order liquid to gyroid transition make this a challenging quest. Here, we investigate the pathway and transition states for the nucleation of a gyroid from the liquid in molecular simulations with a mesogenic binary mixture. We find that the gyroid nuclei at the transition states have a large degree of positional disorder and are not compact, consistent with the low surface free energy of the liquid-gyroid interface. A combination of bond-order parameters for the minor component is best to describe the passage from liquid to gyroid, among those we consider. The committor analyses, however, show that this best coordinate is not perfect and suggests that accounting for the relative ordering of the two interpenetrated networks in infant nuclei, as well as for signatures of ordering in the major component of the mesophase, would improve the accuracy of the reaction coordinate for gyroid formation and its use to evaluate nucleation barriers. To our knowledge, this study is the first to investigate the reaction coordinate and critical nuclei for the formation of any mesophase from an amorphous phase.

Does Local Structure Bias How a Crystal Nucleus Evolves?

The broad scientific and technological importance of crystallization has led to significant research probing and rationalizing crystal nucleation processes. Previous work has generally neglected the possibility of the molecular-level dynamics of individual crystal nuclei coupling to local structures. However, recent experimental work has conjectured that this can occur. Therefore, to address a deficiency in scientific understanding of crystallization, we have probed the nucleation of prototypical single and multicomponent crystals (specifically, ice and mixed gas hydrates). We establish that local structures can bias the evolution of nascent crystal phases on a nanosecond time scale by, for example, promoting the appearance or disappearance of specific crystal motifs and thus reveal a new facet of crystallization behavior. Moreover, we demonstrate structural biases are likely present during crystallization processes beyond ice and gas hydrate formation. Structurally biased dynamics are a lens for understanding existing computational and experimental results while pointing to future opportunities.