Tunable CHA/AEI Zeolite Intergrowths with A Priori Biselective Organic Structure-Directing Agents: Controlling Enrichment and Implications for Selective Catalytic Reduction of NOx.

A novel ab initio methodology based on high-throughput simulations has permitted designing unique biselective organic structure-directing agents (OSDAs) that allow the efficient synthesis of CHA/AEI zeolite intergrowth materials with controlled phase compositions. Distinctive local crystallographic ordering of the CHA/AEI intergrowths was revealed at the nanoscale level using integrated differential phase contrast scanning transmission electron microscopy (iDPC STEM). These novel CHA/AEI materials have been tested for the selective catalytic reduction (SCR) of NOx, presenting an outstanding catalytic performance and hydrothermal stability, even surpassing the performance of the well-established commercial CHA-type catalyst. This methodology opens the possibility for synthetizing new zeolite intergrowths with more complex structures and unique catalytic properties.

Benchmarking binding energy calculations for organic structure-directing agents in pure-silica zeolites.

Molecular modeling plays an important role in the discovery of organic structure-directing agents (OSDAs) for zeolites. By quantifying the intensity of host-guest interactions, it is possible to select cost-effective molecules that maximize binding toward a given zeolite framework. Over the last few decades, a variety of methods and levels of theory have been used to calculate these binding energies. Nevertheless, there is no consensus on the best calculation strategy for high-throughput virtual screening undertakings. In this work, we compare binding affinities from density functional theory (DFT) and Dreiding force field calculations for 272 zeolite-OSDA pairs obtained from static and time-averaged simulations. Enabled by automation software, we show that Dreiding binding energies from the frozen pose method correlate best with DFT energies. They are also less sensitive to the choice of initial lattice parameters and optimization algorithms, as well as less computationally expensive than their time-averaged counterparts. Furthermore, we demonstrate that a broader exploration of the conformation space from molecular dynamics simulations does not provide significant improvements in binding energy trends over the frozen pose method despite being orders of magnitude more expensive. The code and benchmark data are open-sourced and provide robust and computationally efficient guidelines to calculating binding energies in zeolite-OSDA pairs.

Graph similarity drives zeolite diffusionless transformations and intergrowth.

Predicting and directing polymorphic transformations is a critical challenge in zeolite synthesis1-3. Interzeolite transformations enable selective crystallization4-7, but are often too complex to be designed by comparing crystal structures. Here, computational and theoretical tools are combined to both exhaustively data mine polymorphic transformations reported in the literature and analyse and explain interzeolite relations. It was found that crystallographic building units are weak predictors of topology interconversion and insufficient to explain intergrowth. By introducing a supercell-invariant metric that compares crystal structures using graph theory, we show that diffusionless (topotactic and reconstructive) transformations occur only between graph-similar pairs. Furthermore, all the known instances of intergrowth occur between either structurally similar or graph similar frameworks. We identify promising pairs to realize diffusionless transformations and intergrowth, with hundreds of low-distance pairs identified among known zeolites, and thousands of hypothetical frameworks connected to known zeolite counterparts. The theory may enable the understanding and control of zeolite polymorphism.

Temperature-transferable coarse-graining of ionic liquids with dual graph convolutional neural networks.

Computer simulations can provide mechanistic insight into ionic liquids (ILs) and predict the properties of experimentally unrealized ion combinations. However, ILs suffer from a particularly large disparity in the time scales of atomistic and ensemble motion. Coarse-grained models are therefore used in place of costly all-atom simulations, accessing longer time scales and larger systems. Nevertheless, constructing the many-body potential of mean force that defines the structure and dynamics of a coarse-grained system can be complicated and computationally intensive. Machine learning shows great promise for the linked challenges of dimensionality reduction and learning the potential of mean force. To improve the coarse-graining of ILs, we present a neural network model trained on all-atom classical molecular dynamics simulations. The potential of mean force is expressed as two jointly trained neural network interatomic potentials that learn the coupled short-range and many-body long range molecular interactions. These interatomic potentials treat temperature as an explicit input variable to capture its influence on the potential of mean force. The model reproduces structural quantities with high fidelity, outperforms the temperature-independent baseline at capturing dynamics, generalizes to unseen temperatures, and incurs low simulation cost.

One-Pot Synthesis of CHA/ERI-Type Zeolite Intergrowth from a Single Multiselective Organic Structure-Directing Agent.

We report the one-pot synthesis of a chabazite (CHA)/erionite (ERI)-type zeolite intergrowth structure characterized by adjustable extents of intergrowth enrichment and Si/Al molar ratios. This method utilizes readily synthesizable 6-azaspiro[5.6]dodecan-6-ium as the exclusive organic structure-directing agent (OSDA) within a potassium-dominant environment. High-throughput simulations were used to accurately determine the templating energy and molecular shape, facilitating the selection of an optimally biselective OSDA from among thousands of prospective candidates. The coexistence of the crystal phases, forming a distinct structure comprising disk-like CHA regions bridged by ERI-rich pillars, was corroborated via rigorous powder X-ray diffraction and integrated differential-phase contrast scanning transmission electron microscopy (iDPC S/TEM) analyses. iDPC S/TEM imaging further revealed the presence of single offretite layers dispersed within the ERI phase. The ratio of crystal phases between CHA and ERI in this type of intergrowth could be varied systematically by changing both the OSDA/Si and K/Si ratios. Two intergrown zeolite samples with different Si/Al molar ratios were tested for the selective catalytic reduction (SCR) of NOx with NH3, showing competitive catalytic performance and hydrothermal stability compared to that of the industry-standard commercial NH3-SCR catalyst, Cu-SSZ-13, prevalent in automotive applications. Collectively, this work underscores the potential of our approach for the synthesis and optimization of adjustable intergrown zeolite structures, offering competitive alternatives for key industrial processes.

Approaching enzymatic catalysis with zeolites or how to select one reaction mechanism competing with others.

Approaching the level of molecular recognition of enzymes with solid catalysts is a challenging goal, achieved in this work for the competing transalkylation and disproportionation of diethylbenzene catalyzed by acid zeolites. The key diaryl intermediates for the two competing reactions only differ in the number of ethyl substituents in the aromatic rings, and therefore finding a selective zeolite able to recognize this subtle difference requires an accurate balance of the stabilization of reaction intermediates and transition states inside the zeolite microporous voids. In this work we present a computational methodology that, by combining a fast high-throughput screeening of all zeolite structures able to stabilize the key intermediates with a more computationally demanding mechanistic study only on the most promising candidates, guides the selection of the zeolite structures to be synthesized. The methodology presented is validated experimentally and allows to go beyond the conventional criteria of zeolite shape-selectivity.

Mechanistic Insights on Permeation of Water over Iron Cations in Nanoporous Silicon Oxide Films for Selective H2 and O2 Evolution.

Electrocatalysts encapsulated by an ultrathin and semipermeable oxide layer offer a promising avenue for efficient, selective, and cost-effective production of hydrogen through photoelectrochemical water splitting. This architecture is especially attractive for Z-scheme water splitting, for which a nanoporous oxide film can be leveraged to mitigate undesired, yet kinetically facile, reactions involving redox shuttles, such as aqueous iron cations, by limiting transport of these species to catalytically active sites. In this work, molecular dynamics simulations were combined with electrochemical measurements to provide a mechanistic understanding of permeation of water and Fe(III)/Fe(II) redox shuttles through nanoporous SiO2 films. It is shown that even for SiO2 pores with a width as small as 0.8 nm, water does not experience any energy barrier for permeating into the pores due to a favorable interaction with hydrophilic silanol groups on the oxide surface. In contrast, permeation of Fe(III) and Fe(II) into microporous SiO2 pores is limited due to high energy barriers, which stem from a combination of distortion and dehydration of the second and third ion solvation shells. Our simulations and experimental results show that SiO2 coatings can effectively mitigate undesired Fe(III)/Fe(II) redox reactions at underlying electrodes by attenuating permeation of iron cations, while allowing water to permeate and thus participate in water splitting reactions. In a broader context, our study demonstrates that selectivity of solvated cations can be manipulated by controlling the pore size and surface chemistry of oxide films.

Differentiable sampling of molecular geometries with uncertainty-based adversarial attacks.

Neural network (NN) interatomic potentials provide fast prediction of potential energy surfaces, closely matching the accuracy of the electronic structure methods used to produce the training data. However, NN predictions are only reliable within well-learned training domains, and show volatile behavior when extrapolating. Uncertainty quantification methods can flag atomic configurations for which prediction confidence is low, but arriving at such uncertain regions requires expensive sampling of the NN phase space, often using atomistic simulations. Here, we exploit automatic differentiation to drive atomistic systems towards high-likelihood, high-uncertainty configurations without the need for molecular dynamics simulations. By performing adversarial attacks on an uncertainty metric, informative geometries that expand the training domain of NNs are sampled. When combined with an active learning loop, this approach bootstraps and improves NN potentials while decreasing the number of calls to the ground truth method. This efficiency is demonstrated on sampling of kinetic barriers, collective variables in molecules, and supramolecular chemistry in zeolite-molecule interactions, and can be extended to any NN potential architecture and materials system.

Data-Driven Design of Biselective Templates for Intergrowth Zeolites.

Zeolites are inorganic materials with wide industrial applications due to their topological diversity. Tailoring confinement effects in zeolite pores, for instance by crystallizing intergrown frameworks, can improve their catalytic and transport properties, but controlling zeolite crystallization often relies on heuristics. In this work, we use computational simulations and data mining to design organic structure-directing agents (OSDAs) to favor the synthesis of intergrown zeolites. First, we propose design principles to identify OSDAs which are selective toward both end members of the disordered structure. Then, we mine a database of hundreds of thousands of zeolite-OSDA pairs and downselect OSDA candidates to synthesize known intergrowth zeolites such as CHA/AFX, MTT/TON, and BEC/ISV. The computationally designed OSDAs balance phase competition metrics and shape selectivity toward the frameworks, thus bypassing expensive dual-OSDA approaches typically used in the synthesis of intergrowths. Finally, we propose potential OSDAs to obtain hypothesized disordered frameworks such as AEI/SAV. This work may accelerate zeolite discovery through data-driven synthesis optimization and design.

Discovering Relationships between OSDAs and Zeolites through Data Mining and Generative Neural Networks.

Organic structure directing agents (OSDAs) play a crucial role in the synthesis of micro- and mesoporous materials especially in the case of zeolites. Despite the wide use of OSDAs, their interaction with zeolite frameworks is poorly understood, with researchers relying on synthesis heuristics or computationally expensive techniques to predict whether an organic molecule can act as an OSDA for a certain zeolite. In this paper, we undertake a data-driven approach to unearth generalized OSDA-zeolite relationships using a comprehensive database comprising of 5,663 synthesis routes for porous materials. To generate this comprehensive database, we use natural language processing and text mining techniques to extract OSDAs, zeolite phases, and gel chemistry from the scientific literature published between 1966 and 2020. Through structural featurization of the OSDAs using weighted holistic invariant molecular (WHIM) descriptors, we relate OSDAs described in the literature to different types of cage-based, small-pore zeolites. Lastly, we adapt a generative neural network capable of suggesting new molecules as potential OSDAs for a given zeolite structure and gel chemistry. We apply this model to CHA and SFW zeolites generating several alternative OSDA candidates to those currently used in practice. These molecules are further vetted with molecular mechanics simulations to show the model generates physically meaningful predictions. Our model can automatically explore the OSDA space, reducing the amount of simulation or experimentation needed to find new OSDA candidates.

A priori control of zeolite phase competition and intergrowth with high-throughput simulations.

Zeolites are versatile catalysts and molecular sieves with large topological diversity, but managing phase competition in zeolite synthesis is an empirical, labor-intensive task. In this work, we controlled phase selectivity in templated zeolite synthesis from first principles by combining high-throughput atomistic simulations, literature mining, human-computer interaction, synthesis, and characterization. Proposed binding metrics distilled from more than 586,000 zeolite-molecule simulations reproduced the extracted literature and rationalized framework competition in the design of organic structure-directing agents. Energetic, geometric, and electrostatic descriptors of template molecules were found to regulate synthetic accessibility windows and aluminum distributions in pure-phase zeolites. Furthermore, these parameters allowed us to realize an intergrowth zeolite through a single bi-selective template. The computation-first approach enables control of both zeolite synthesis and structure composition using a priori theoretical descriptors.