Prospects of transition interface sampling simulations for the theoretical study of zeolite synthesis

PyRETIS 3: Conquering rare and slow events without boundaries

We present and discuss the advancements made in PyRETIS 3, the third instalment of our Python library for an efficient and user-friendly rare event simulation, focused to execute molecular simulations with replica exchange transition interface sampling (RETIS) and its variations. Apart from a general rewiring of the internal code towards a more modular structure, several recently developed sampling strategies have been implemented. These include recently developed Monte Carlo moves to increase path decorrelation and convergence rate, and new ensemble definitions to handle the challenges of long-lived metastable states and transitions with unbounded reactant and product states. Additionally, the post-analysis software PyVisa is now embedded in the main code, allowing fast use of machine-learning algorithms for clustering and visualising collective variables in the simulation data.

Use of Multiscale Data-Driven Surrogate Models for Flowsheet Simulation of an Industrial Zeolite Production Process

The production of catalysts such as zeolites is a complex multiscale and multi-step process. Various material properties, such as particle size or moisture content, as well as operating parameters—e.g., temperature or amount and composition of input material flows—significantly affect the outcome of each process step, and hence determine the properties of the final product. Therefore, the design and optimization of such processes is a complex task, which can be greatly facilitated with the help of numerical simulations. This contribution presents a modeling framework for the dynamic flowsheet simulation of a zeolite production sequence consisting of four stages: precipitation in a batch reactor; concentration and washing in a block of centrifuges; formation of droplets and drying in a spray dryer; and burning organic residues in a chain of rotary kilns. Various techniques and methods were used to develop the applied models. For the synthesis in the reactor, a multistage strategy was used, comprising discrete element method simulations, data-driven surrogate modeling, and population balance modeling. The concentration and washing stage consisted of several multicompartment decanter centrifuges alternating with water mixers. The drying is described by a co–current spray dryer model developed by applying a two-dimensional population balance approach. For the rotary kilns, a multi-compartment model was used, which describes the gas–solid reaction in the counter–current solids and gas flows.

Identifying a Feasible Transition Pathway between Two Conformational States for a Protein

Proteins usually need to transit between different conformational states to fulfill their biological functions. In the mechanistic study of such transition processes by molecular dynamics simulations, identification of the minimum free energy path (MFEP) can substantially reduce the sampling space, thus enabling rigorous thermodynamic evaluation of the process. Conventionally, the MFEP is derived by iterative local optimization from an initial path, which is typically generated by simple brute force techniques like the targeted molecular dynamics (tMD). Therefore, the quality of the initial path determines the successfulness of MFEP estimation. In this work, we propose a method to improve derivation of the initial path. Through iterative relaxation-biasing simulations in a bidirectional manner, this method can construct a feasible transition pathway connecting two known states for a protein. Evaluation on small, fast-folding proteins against long equilibrium trajectories supports the good sampling efficiency of our method. When applied to larger proteins including the catalytic domain of human c-Src kinase as well as the converter domain of myosin VI, the paths generated by our method deviate significantly from those computed with the generic tMD approach. More importantly, free energy profiles and intermediate states obtained from our paths exhibit remarkable improvements over those from tMD paths with respect to both physical rationality and consistency with a priori knowledge.

Advances in theory and their application within the field of zeolite chemistry

Zeolites are versatile and fascinating materials which are vital for a wide range of industries, due to their unique structural and chemical properties, which are the basis of applications in gas separation, ion exchange and catalysis. Given their economic impact, there is a powerful incentive for smart design of new materials with enhanced functionalities to obtain the best material for a given application. Over the last decades, theoretical modeling has matured to a level that model guided design has become within reach. Major hurdles have been overcome to reach this point and almost all contemporary methods in computational materials chemistry are actively used in the field of modeling zeolite chemistry and applications. Integration of complementary modeling approaches is necessary to obtain reliable predictions and rationalizations from theory. A close synergy between experimentalists and theoreticians has led to a deep understanding of the complexity of the system at hand, but also allowed the identification of shortcomings in current theoretical approaches. Inspired by the importance of zeolite characterization which can now be performed at the single atom and single molecule level from experiment, computational spectroscopy has grown in importance in the last decade. In this review most of the currently available modeling tools are introduced and illustrated on the most challenging problems in zeolite science. Directions for future model developments will be given.