1
|
Vervust W, Zhang DT, Ghysels A, Roet S, van Erp TS, Riccardi E. PyRETIS 3: Conquering rare and slow events without boundaries. J Comput Chem 2024; 45:1224-1234. [PMID: 38345082 DOI: 10.1002/jcc.27319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Revised: 01/16/2024] [Accepted: 01/18/2024] [Indexed: 04/19/2024]
Abstract
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.
Collapse
Affiliation(s)
- Wouter Vervust
- IBiTech-BioMMedA Group, Ghent University, Ghent, Belgium
| | - Daniel T Zhang
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| | - An Ghysels
- IBiTech-BioMMedA Group, Ghent University, Ghent, Belgium
| | - Sander Roet
- Department of Chemistry, Utrecht University, Utrecht, The Netherlands
| | - Titus S van Erp
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| | - Enrico Riccardi
- Department of Energy Resources, University of Stavanger, Stavanger, Norway
| |
Collapse
|
2
|
Matinfar M, Nychka JA. A review of sodium silicate solutions: Structure, gelation, and syneresis. Adv Colloid Interface Sci 2023; 322:103036. [PMID: 37952363 DOI: 10.1016/j.cis.2023.103036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Revised: 10/05/2023] [Accepted: 10/28/2023] [Indexed: 11/14/2023]
Abstract
Sodium silicate solutions, also known as waterglass, have been found to have remarkable utility in a variety of applications. The cumulative weight of evidence from 70 years of varied analysis indicates that silicate solutions consist of a wide range of species, from monomers through oligomers, up to colloids. Moreover, the structure and distribution of these species are greatly dependent upon many parameters, such as solute concentrations, silica to alkali ratio, pH, and temperature. The most interesting and characteristic property of silicate solutions is their ability to form silica gels. Overall, despite extensive research using different spectroscopic and scattering techniques, many questions related to sodium silicate's dynamic structure, stability, polymerization, and gelation remain difficult to answer. The multitude of simultaneous reactions which restructure the silicate species at the atomic scale in response to variation in solution and environmental parameters, makes it difficult to investigate the individual events using only experimental data. Molecular modelling provides an alternative way to study the unknown areas in the aqueous silicate and silica gel systems, generating key insights into the chemical reactions at microscopic length scales. However, sufficient sampling remains a challenge for the practical use of molecular simulation for these systems. Based on both experimental and modelling studies, this review provides a detailed discussion over the structure and speciation of sodium silicate solutions, their gelation mechanism and kinetics, and the syneresis phenomenon. The goal is not only to review the current level of understanding of sodium silicate solutions, silica gels and characterization techniques suitable for studying them, but also to identify the gaps in the literature and open up opportunities for advancing knowledge about these complex systems. We believe that the future direction of research should be toward correlating atomistic, molecular, and meso-scale level details of interactions and reactions in silicate solution and establishing a fundamental understanding of its gelation mechanism and kinetics. We believe that this knowledge could eliminate the "trial and error" approach in manufacturing, and improve structural control in the synthesis of important materials derived from these solutions, such as silica gels and zeolites.
Collapse
Affiliation(s)
- Marzieh Matinfar
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 1H9, Canada.
| | - John A Nychka
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 1H9, Canada
| |
Collapse
|
3
|
Ho TH, Do TH, Tong HD, Meijer EJ, Trinh TT. The Role of Chloride ion in the Silicate Condensation Reaction from ab Initio Molecular Dynamics Simulations. J Phys Chem B 2023; 127:7748-7757. [PMID: 37647302 PMCID: PMC10510376 DOI: 10.1021/acs.jpcb.3c04256] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2023] [Revised: 08/24/2023] [Indexed: 09/01/2023]
Abstract
The comprehension of silicate oligomer formation during the initial stage of zeolite synthesis is of significant importance. In this study, we investigated the effect of chloride ions (Cl-) on silicate oligomerization using ab initio molecular dynamics simulations with explicit water molecules. The results show that the presence of Cl- increases the free energy barriers of all reactions compared to the case without the anion. The formation of the 4-ring structure has the lowest free energy barrier (73 kJ/mol), while the formation of the 3-ring structure has the highest barrier (98 kJ/mol) in the presence of Cl-. These findings suggest that Cl- suppresses the formation of 3-rings and favors the formation of larger oligomers in the process of zeolite synthesis. Our study provides important insights into the directing role of Cl- in silicate oligomerization by regulating thermodynamic and kinetic parameters. An important point to consider is the impact of the anion on aqueous reactions, particularly in altering the hydrogen bond network around reactive species. These results also provide a basis for further studies of the formations of larger silicate oligomers in solution.
Collapse
Affiliation(s)
- Thi H. Ho
- Laboratory
for Computational Physics Institute for Computational Science and
Artificial Intelligence, Van Lang University, Ho Chi Minh City 700000, Vietnam
- Faculty
of Mechanical - Electrical and Computer Engineering School of Technology, Van Lang University, Ho Chi Minh City 700000, Vietnam
| | - Tuong Ha Do
- Faculty
of Applied Sciences, Ton Duc Thang University, 19 Nguyen Huu Tho, Tan Phong ward
District 7, Ho Chi Minh City 700000, Vietnam
| | - Hien Duy Tong
- Faculty
of Engineering, Vietnamese-German University
(VGU), Thu Dau
Mot City, Binh Duong Province 75000, Vietnam
| | - Evert Jan Meijer
- Van
’t Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam 1012 WX, The Netherlands
| | - Thuat T. Trinh
- Porelab,
Department of Chemistry, Norwegian University
of Science and Technology, NO-7491 Trondheim, Norway
| |
Collapse
|
4
|
Do TH, Tong HD, Tran KQ, Meijer EJ, Trinh TT. Insight into the role of excess hydroxide ions in silicate condensation reactions. Phys Chem Chem Phys 2023; 25:12723-12733. [PMID: 37114395 DOI: 10.1039/d3cp00475a] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
The formation of silicate oligomers in the early stages is key to zeolite synthesis. The pH and the presence of hydroxide ions are important in regulating the reaction rate and the dominant species in solutions. This paper describes the formation of silicate species, from dimers to 4-membered rings, using ab initio molecular dynamics simulations in explicit water molecules with an excess hydroxide ion. The thermodynamic integration method was used to calculate the free energy profile of the condensation reactions. The hydroxide group's role is not only to control the pH of the environment, but also to actively participate in the condensation reaction. The results show that the most favorable reactions are linear-tetramer and 4-membered-ring formation, with overall barriers of 71 kJ mol-1 and 73 kJ mol-1, respectively. The formation of trimeric silicate, with the largest free-energy barrier of 102 kJ mol-1, is the rate-limiting step under these conditions. The excess hydroxide ion aids in the stabilization of the 4-membered-ring structure over the 3-membered-ring structure. Due to a relatively high free-energy barrier, the 4-membered ring is the most difficult of the small silicate structures to dissolve in the backward reaction. This study is consistent with the experimental observation that silicate growth in zeolite synthesis is slower in a very-high-pH environment.
Collapse
Affiliation(s)
- Tuong Ha Do
- Faculty of Applied Sciences, Ton Duc Thang University, 19 Nguyen Huu Tho, Tan Phong ward, District 7, Ho Chi Minh City, Vietnam.
| | - Hien Duy Tong
- Faculty of Engineering, Vietnamese-German University (VGU), Thu Dau Mot City, Binh Duong Province, Vietnam
| | - Khanh-Quang Tran
- Department of Energy and Process Engineering, Norwegian University of Science and Technology, Kolbjørn Hejes vei 1B, NO- 7491 Trondheim, Norway
| | - Evert Jan Meijer
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands.
| | - Thuat T Trinh
- Porelab, Department of Chemistry, Norwegian University of Science and Technology, Høgskoleringen 5, 7491-Trondheim, Norway.
| |
Collapse
|
5
|
Hong S, Mallette AJ, Neeway JJ, Motkuri RK, Rimer JD, Mpourmpakis G. Understanding formation thermodynamics of structurally diverse zeolite oligomers with first principles calculations. Dalton Trans 2023; 52:1301-1315. [PMID: 36625388 DOI: 10.1039/d2dt02764j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The mechanisms of many zeolitic processes, including nucleation and interzeolite transformation, are not fully understood owing to complex growth mixtures that obfuscate in situ monitoring of molecular events. In this work, we provide insights into zeolite chemistry by investigating the formation thermodynamics of small zeolitic species using first principles calculations. We systematically study how formation energies of pure-silicate and aluminosilicate species differ by structure type and size, temperature, and the presence of alkali or alkaline earth metal cations (Na+, K+, and Ca2+). Highly condensed (cage-like) species are found to be strongly preferred to simple rings in the pure-silicate system, and this thermodynamic preference increases with temperature. Introducing aluminum leads to more favorable formation thermodynamics for all species. Moreover, for species with a low Si/Al ratio (≤2), a thermodynamic preference does not exist among structure types; instead, a pool of diverse aluminosilicate structures compete in formation. Metal cation effects strongly depend on the presence of aluminum, cage size, cation type, and location, since each of these factors can alter electrostatic interactions between cations and zeolitic species. We reveal that confined metal cations may destabilize pure-silicate cages due to localized interactions; conversely, they stabilize aluminosilicates due to strong cation-framework attractions in sufficiently large cages. Importantly, this work rationalizes a series of experimental observations and can potentially guide efforts for controlling zeolite nucleation/crystallization processes.
Collapse
Affiliation(s)
- Sungil Hong
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, PA, 15261, USA.
| | - Adam J Mallette
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, TX, 77204, USA
| | - James J Neeway
- Energy and Environment Directorate, Pacific Northwest National Laboratory, Richland, WA, 99352, USA
| | - Radha Kishan Motkuri
- Energy and Environment Directorate, Pacific Northwest National Laboratory, Richland, WA, 99352, USA
| | - Jeffrey D Rimer
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, TX, 77204, USA
| | - Giannis Mpourmpakis
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, PA, 15261, USA.
| |
Collapse
|
6
|
Zhang DT, Riccardi E, van Erp TS. Enhanced path sampling using subtrajectory Monte Carlo moves. J Chem Phys 2023; 158:024113. [PMID: 36641412 DOI: 10.1063/5.0127249] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Path sampling allows the study of rare events, such as chemical reactions, nucleation, and protein folding, via a Monte Carlo (MC) exploration in path space. Instead of configuration points, this method samples short molecular dynamics (MD) trajectories with specific start- and end-conditions. As in configuration MC, its efficiency highly depends on the types of MC moves. Since the last two decades, the central MC move for path sampling has been the so-called shooting move in which a perturbed phase point of the old path is propagated backward and forward in time to generate a new path. Recently, we proposed the subtrajectory moves, stone-skipping (SS) and web-throwing, that are demonstrably more efficient. However, the one-step crossing requirement makes them somewhat more difficult to implement in combination with external MD programs or when the order parameter determination is expensive. In this article, we present strategies to address the issue. The most generic solution is a new member of subtrajectory moves, wire fencing (WF), that is less thrifty than the SS but more versatile. This makes it easier to link path sampling codes with external MD packages and provides a practical solution for cases where the calculation of the order parameter is expensive or not a simple function of geometry. We demonstrate the WF move in a double-well Langevin model, a thin film breaking transition based on classical force fields, and a smaller ruthenium redox reaction at the ab initio level in which the order parameter explicitly depends on the electron density.
Collapse
Affiliation(s)
- Daniel T Zhang
- Norwegian University of Science and Technology, Department of Chemistry, NO-7491 Trondheim, Norway
| | - Enrico Riccardi
- Department of Informatics, UiO, Gaustadalléen 23B, 0373 Oslo, Norway
| | - Titus S van Erp
- Norwegian University of Science and Technology, Department of Chemistry, NO-7491 Trondheim, Norway
| |
Collapse
|
7
|
Hare SR, Pfaendtner J. Elucidating the role of catalytic amino acid residues in the peptide-mediated silica oligomerization reaction mechanism. Phys Chem Chem Phys 2022; 24:3664-3674. [PMID: 35080220 DOI: 10.1039/d1cp04542c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Understanding the detailed mechanism by which the proteins of marine diatoms such as silaffins are able to control the morphology of silica oligomers has eluded synthetic chemists and materials scientists for decades. In this study, we use DFT calculations to determine how individual amino acid residues of silaffin catalyze silica dimerization. The reaction network for formation of a silica dimer was explored using several different small molecules, including water, guanidinium ions, and methylammonium ions, the latter two molecules representing analogs of arginine and lysine, both of which are known to play critical roles in enabling the catalytic function of naturally occurring protein and synthetic analogs of silaffin. It was found that the lysine analog selectively lowers the energy of a direct water removal pathway for silicate dimerization. Comparing the energy landscapes and mechanisms for various catalysts for this reaction provides direct evidence for the role of lysine side chains of silaffins in the oligmerization of silica.
Collapse
Affiliation(s)
- Stephanie R Hare
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, USA.
| | - Jim Pfaendtner
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, USA. .,Physical Science Division, Pacific Northwest National Laboratory, Richmond, Washington 99354, USA
| |
Collapse
|
8
|
Maghrebi K, Gam S, Hammami B, Alsadiri A, Abderrabba M, Messaoudi S. Exploration of the Mechanism of the Dimerization of Hydroxymethylsilanetriol Using Electronic Structure Methods. ACS OMEGA 2022; 7:2661-2670. [PMID: 35097264 PMCID: PMC8792941 DOI: 10.1021/acsomega.1c05027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/11/2021] [Accepted: 12/24/2021] [Indexed: 06/14/2023]
Abstract
Hydroxymethylsilanetriol undergoes condensation reactions to form new structures with an organic part in the formed bridges. As a first step to explore the formation of these bridges, we studied the corresponding mechanisms using simple models and theoretical methods. Three mechanisms were studied for the formation of dimers of hydroxymethylsilanetriol with bridges: Si-O-C-Si, Si-O-Si, and Si-C-O-C-Si. Energies are calculated using M06/6-311+G(d,p) single-point calculations on B3LYP-optimized geometries in solution and including B3LYP thermodynamic corrections. The first mechanism for the formation of the Si-O-C-Si bridge consists of one step. The second mechanism for the formation of the Si-O-Si bridge consists of two steps. The barrier for the last mechanism for the formation of the Si-C-O-C-Si bridge is too high and cannot occur at room temperature. The energy barriers are 31.8, 27.6, and 65.9 kcal mol-1 for the first, second, and third mechanisms, respectively. When adding one explicit water molecule, these energies are 25.9, 22.9, and 80.3 kcal mol-1, respectively. The first and second mechanisms can occur at room temperature, which is in agreement with the experimental results.
Collapse
Affiliation(s)
- Khaoula Maghrebi
- Faculty
of Sciences of Bizerte FSB, University of
Carthage, 7021 Jarzouna, Tunisia
- Laboratory
of Materials, Molecules and Applications, IPEST, University of Carthage, Sidi Bou Said Road, B.P. 51 2070 La Marsa, Tunisia
| | - Safa Gam
- Faculty
of Sciences of Bizerte FSB, University of
Carthage, 7021 Jarzouna, Tunisia
- Laboratory
of Materials, Molecules and Applications, IPEST, University of Carthage, Sidi Bou Said Road, B.P. 51 2070 La Marsa, Tunisia
| | - Bechir Hammami
- Faculty
of Sciences of Bizerte FSB, University of
Carthage, 7021 Jarzouna, Tunisia
- Department
of Chemistry, College of Science, Qassim
University, Buraidah 51452, Saudi Arabia
| | - Atheer Alsadiri
- Department
of Chemistry, College of Science, Qassim
University, Buraidah 51452, Saudi Arabia
| | - Manef Abderrabba
- Laboratory
of Materials, Molecules and Applications, IPEST, University of Carthage, Sidi Bou Said Road, B.P. 51 2070 La Marsa, Tunisia
| | - Sabri Messaoudi
- Faculty
of Sciences of Bizerte FSB, University of
Carthage, 7021 Jarzouna, Tunisia
- Laboratory
of Materials, Molecules and Applications, IPEST, University of Carthage, Sidi Bou Said Road, B.P. 51 2070 La Marsa, Tunisia
- Department
of Chemistry, College of Science, Qassim
University, Buraidah 51452, Saudi Arabia
| |
Collapse
|
9
|
Xia H, Fan X, Zhang J, He H, Guo Q. Theoretical Studies on the Oligomerization of Silicate Species in Basic Solution. J Phys Chem A 2021; 125:8827-8835. [PMID: 34587739 DOI: 10.1021/acs.jpca.1c06203] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The possible reaction pathways of silicate species to linear- and ring-structure oligomers up to silicate hexamers in the basic medium have been studied using the density functional theory. The calculations were performed at the ωB97XD/6-31+G(d,p) level, and the integral equation formalism polarizable continuum model was employed to simulate the solvent effects, and it was found that they are appropriate in exploring the reaction mechanism of silicate species condensation. There are two steps in the anionic silicate condensation reactions: the SiO-Si bond formation step and the dehydration step. Moreover, the latter is the rate-limiting step for most of the reaction pathways except for the cyclization reaction of the linear pentamer to the 5-ring. The short linear oligomers would be likely formed from the reaction between monomers and oligomers, while the longer ones are easily formed through the reactions between short oligomers. The 4-ring and branched 5-ring oligomers are found to be formed very favorable both in kinetic and thermodynamic and could have great influences on the initial stage of zeolite synthesis. The intramolecular and intermolecular hydrogen bond effect of silicate species is an important factor affecting the reaction mechanism.
Collapse
Affiliation(s)
- HongQiang Xia
- State Key Laboratory of High-efficiency Utilization of Coal and Green Chemical Engineering, Ningxia University, Yinchuan 750021, P. R. China
| | - Xiyu Fan
- State Key Laboratory of High-efficiency Utilization of Coal and Green Chemical Engineering, Ningxia University, Yinchuan 750021, P. R. China
| | - Jianli Zhang
- State Key Laboratory of High-efficiency Utilization of Coal and Green Chemical Engineering, Ningxia University, Yinchuan 750021, P. R. China
| | - Haiming He
- State Key Laboratory of High-efficiency Utilization of Coal and Green Chemical Engineering, Ningxia University, Yinchuan 750021, P. R. China
| | - Qingjie Guo
- State Key Laboratory of High-efficiency Utilization of Coal and Green Chemical Engineering, Ningxia University, Yinchuan 750021, P. R. China
| |
Collapse
|
10
|
Roet S, Daub CD, Riccardi E. Chemistrees: Data-Driven Identification of Reaction Pathways via Machine Learning. J Chem Theory Comput 2021; 17:6193-6202. [PMID: 34555907 PMCID: PMC8515787 DOI: 10.1021/acs.jctc.1c00458] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
![]()
We propose to analyze
molecular dynamics (MD) output via a supervised machine
learning (ML) algorithm, the decision tree.
The approach aims to identify the predominant geometric features which
correlate with trajectories that transition between two arbitrarily
defined states. The data-driven algorithm aims to identify these features
without the bias of human “chemical intuition”. We demonstrate
the method by analyzing the proton exchange reactions in formic acid
solvated in small water clusters. The simulations were performed with ab initio MD combined with a method to efficiently sample
the rare event, path sampling. Our ML analysis identified relevant
geometric variables involved in the proton transfer reaction and how
they may change as the number of solvating water molecules changes.
Collapse
Affiliation(s)
- Sander Roet
- Department of Chemistry, Norwegian University of Science and Technology, Høgskoleringen 5, 7491 Trondheim, Norway
| | - Christopher D Daub
- Department of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| | - Enrico Riccardi
- Department of Informatics, UiO, Gaustadalléen 23B, 0373 Oslo, Norway
| |
Collapse
|
11
|
Liu X, Liu C, Feng Z, Meng C. The Promoter Role of Amines in the Condensation of Silicic Acid: A First-Principles Investigation. ACS OMEGA 2021; 6:22811-22819. [PMID: 34514252 PMCID: PMC8427787 DOI: 10.1021/acsomega.1c03235] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/21/2021] [Accepted: 08/18/2021] [Indexed: 06/13/2023]
Abstract
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 SN2-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.
Collapse
Affiliation(s)
- Xin Liu
- State Key Laboratory of Fine
Chemicals and Department of Chemistry, Dalian
University of Technology, No. 2 Linggong Road, Dalian 116024, P. R. China
| | - Cai Liu
- State Key Laboratory of Fine
Chemicals and Department of Chemistry, Dalian
University of Technology, No. 2 Linggong Road, Dalian 116024, P. R. China
| | - Zhe Feng
- State Key Laboratory of Fine
Chemicals and Department of Chemistry, Dalian
University of Technology, No. 2 Linggong Road, Dalian 116024, P. R. China
| | - Changgong Meng
- State Key Laboratory of Fine
Chemicals and Department of Chemistry, Dalian
University of Technology, No. 2 Linggong Road, Dalian 116024, P. R. China
| |
Collapse
|
12
|
Ciantar M, Trinh TT, Michel C, Sautet P, Mellot-Draznieks C, Nieto-Draghi C. Impact of Organic Templates on the Selective Formation of Zeolite Oligomers. Angew Chem Int Ed Engl 2021; 60:7111-7116. [PMID: 33237634 DOI: 10.1002/anie.202014027] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Indexed: 11/07/2022]
Abstract
Zeolites are essential materials to industry due to their adsorption and catalytic properties. The best current approach to prepare a targeted zeolite still relies on trial and error's synthetic procedures since a rational understanding of the impact of synthesis variables on the final structures is still missing. To discern the role of a variety of organic templates, we perform simulations of the early stages of condensation of silica oligomers by combining DFT, Brønsted-Evans-Polanyi relationships and kinetic Monte Carlo simulations. We investigate an extended reaction path mechanism including 258 equilibrium reactions and 242 chemical species up to silica octamers, comparing the computed concentrations of Si oligomers with 29 SI NMR experimental data. The effect of the templating agent is linked to the modification of the intramolecular H-bond network in the growing oligomer, which produces higher concentration of 4-membered ring intermediates, precursors of the key double-four ring building blocks present on more than 39 known zeolite topologies.
Collapse
Affiliation(s)
- Marine Ciantar
- IFP Energies nouvelles, 1-4 Avenue de Bois-Préau, 92852, Rueil-Malmaison, France
| | - Thuat T Trinh
- Univ Lyon, Ens de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratoire de Chimie, 69342, Lyon, France
| | - Carine Michel
- Univ Lyon, Ens de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratoire de Chimie, 69342, Lyon, France
| | - Philippe Sautet
- Univ Lyon, Ens de Lyon, CNRS UMR 5182, Université Claude Bernard Lyon 1, Laboratoire de Chimie, 69342, Lyon, France.,Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Los Angeles, CA, 90095, USA.,Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, CA, 90095, USA
| | - Caroline Mellot-Draznieks
- Laboratoire de Chimie des Processus Biologiques. UMR CNRS 8229, Collège de France, Sorbonne Université, PSL Research University, 11 Place Marcelin Berthelot, 75231, Paris, Cedex 05, France
| | - Carlos Nieto-Draghi
- IFP Energies nouvelles, 1-4 Avenue de Bois-Préau, 92852, Rueil-Malmaison, France
| |
Collapse
|
13
|
Bolhuis PG, Swenson DWH. Transition Path Sampling as Markov Chain Monte Carlo of Trajectories: Recent Algorithms, Software, Applications, and Future Outlook. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202000237] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Peter G. Bolhuis
- Amsterdam Center for Multiscale Modeling van 't Hoff Institute for Molecular Sciences University of Amsterdam PO Box 94157 1090 GD Amsterdam The Netherlands
| | - David W. H. Swenson
- Centre Blaise Pascal Ecole Normale Superieure 46, allée d'Italie 69364 Lyon Cedex 07 France
| |
Collapse
|
14
|
Ciantar M, Trinh TT, Michel C, Sautet P, Mellot‐Draznieks C, Nieto‐Draghi C. Impact of Organic Templates on the Selective Formation of Zeolite Oligomers. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202014027] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Marine Ciantar
- IFP Energies nouvelles 1–4 Avenue de Bois-Préau 92852 Rueil-Malmaison France
| | - Thuat T. Trinh
- Univ Lyon Ens de Lyon, CNRS UMR 5182 Université Claude Bernard Lyon 1 Laboratoire de Chimie 69342 Lyon France
| | - Carine Michel
- Univ Lyon Ens de Lyon, CNRS UMR 5182 Université Claude Bernard Lyon 1 Laboratoire de Chimie 69342 Lyon France
| | - Philippe Sautet
- Univ Lyon Ens de Lyon, CNRS UMR 5182 Université Claude Bernard Lyon 1 Laboratoire de Chimie 69342 Lyon France
- Department of Chemical and Biomolecular Engineering University of California, Los Angeles Los Angeles CA 90095 USA
- Department of Chemistry and Biochemistry University of California, Los Angeles Los Angeles CA 90095 USA
| | - Caroline Mellot‐Draznieks
- Laboratoire de Chimie des Processus Biologiques. UMR CNRS 8229 Collège de France Sorbonne Université PSL Research University 11 Place Marcelin Berthelot 75231 Paris, Cedex 05 France
| | - Carlos Nieto‐Draghi
- IFP Energies nouvelles 1–4 Avenue de Bois-Préau 92852 Rueil-Malmaison France
| |
Collapse
|
15
|
Aarøen O, Kiær H, Riccardi E. PyVisA
: Visualization and Analysis of path sampling trajectories. J Comput Chem 2020; 42:435-446. [DOI: 10.1002/jcc.26467] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2020] [Revised: 11/03/2020] [Accepted: 11/26/2020] [Indexed: 12/20/2022]
Affiliation(s)
- Ola Aarøen
- Department of Biotechnology and Food Science Norwegian University of Science and Technology Trondheim Norway
| | - Henrik Kiær
- Department of Chemistry Norwegian University of Science and Technology Trondheim Norway
| | - Enrico Riccardi
- Department of Chemistry Norwegian University of Science and Technology Trondheim Norway
| |
Collapse
|
16
|
Mai NL, Do HT, Hoang NH, Nguyen AH, Tran KQ, Meijer EJ, Trinh TT. Elucidating the Role of Tetraethylammonium in the Silicate Condensation Reaction from Ab Initio Molecular Dynamics Simulations. J Phys Chem B 2020; 124:10210-10218. [PMID: 33119320 PMCID: PMC7735729 DOI: 10.1021/acs.jpcb.0c06607] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The understanding of the formation of silicate oligomers in the initial stage of zeolite synthesis is important. The use of organic structure-directing agents (OSDAs) is known to be a key factor in the formation of different silicate species and the final zeolite structure. For example, tetraethylammonium ion (TEA+) is a commonly used organic template for zeolite synthesis. In this study, ab initio molecular dynamics (AIMD) simulation is used to provide an understanding of the role of TEA+ in the formation of various silicate oligomers, ranging from dimer to 4-ring. Calculated free-energy profiles of the reaction pathways show that the formation of a 4-ring structure has the highest energy barrier (97 kJ/mol). The formation of smaller oligomers such as dimer, trimer, and 3-ring has lower activation barriers. The TEA+ ion plays an important role in regulating the predominant species in solution via its coordination with silicate structures during the condensation process. The kinetics and thermodynamics of the oligomerization reaction indicate a more favorable formation of the 3-ring over the 4-ring structure. The results from AIMD simulations are in line with the experimental observation that TEA+ favors the 3-ring and double 3-ring in solution. The results of this study imply that the role of OSDAs is not only important for the host-guest interaction but also crucial for controlling the reactivity of different silicate oligomers during the initial stage of zeolite formation.
Collapse
Affiliation(s)
- Ngoc Lan Mai
- Faculty of Applied Sciences, Ton Duc Thang University, 19 Nguyen Huu Tho Str., Tan Phong Ward, District 7, Ho Chi Minh City, Vietnam
| | - Ha T Do
- Faculty of Applied Sciences, Ton Duc Thang University, 19 Nguyen Huu Tho Str., Tan Phong Ward, District 7, Ho Chi Minh City, Vietnam
| | - Nguyen Hieu Hoang
- Department of Materials and Nanotechnology, SINTEF Industry, 7034 Trondheim, Norway
| | - Anh H Nguyen
- Electrical Engineering and Computer Sciences, University of California Irvine, Irvine, California 92697, United States
| | - Khanh-Quang Tran
- Department of Energy and Process Engineering, Norwegian University of Science and Technology, Kolbjørn Hejes vei 1B, 7491 Trondheim, Norway
| | - Evert Jan Meijer
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam 1012 WX, The Netherlands
| | - Thuat T Trinh
- Department of Civil and Environmental Engineering, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| |
Collapse
|
17
|
Freeman EE, Neeway JJ, Motkuri RK, Rimer JD, Mpourmpakis G. Understanding initial zeolite oligomerization steps with first principles calculations. AIChE J 2020. [DOI: 10.1002/aic.17107] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Emily E. Freeman
- Department of Chemical and Petroleum Engineering University of Pittsburgh Pittsburgh Pennsylvania USA
| | - James J. Neeway
- Energy and Environment Directorate Pacific Northwest National Laboratory Richland Washington USA
| | - Radha Kishan Motkuri
- Energy and Environment Directorate Pacific Northwest National Laboratory Richland Washington USA
| | - Jeffrey D. Rimer
- Department of Chemical and Biomolecular Engineering University of Houston Houston Texas USA
| | - Giannis Mpourmpakis
- Department of Chemical and Petroleum Engineering University of Pittsburgh Pittsburgh Pennsylvania USA
| |
Collapse
|
18
|
|
19
|
Daub CD, Riccardi E, Hänninen V, Halonen L. Path sampling for atmospheric reactions: formic acid catalysed conversion of SO3 + H2O to H2SO4. PEERJ PHYSICAL CHEMISTRY 2020. [DOI: 10.7717/peerj-pchem.7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Atmospheric reactions, hitherto studied computationally mainly with static computations in conjunction with transition state theories, can be further described via path sampling calculations. Here we report on an exploratory study of the formic acid catalysed hydrolysis of SO3 to produce H2SO4. We demonstrate that precise measurements of rate coefficients can be obtained in principle for such reactions with an acceptable expenditure of computational resources, and that new insights into the reaction can be obtained by the analysis of the path space explored via path sampling techniques.
Collapse
Affiliation(s)
| | - Enrico Riccardi
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim, Norway
| | - Vesa Hänninen
- Department of Chemistry, University of Helsinki, Helsinki, Finland
| | - Lauri Halonen
- Department of Chemistry, University of Helsinki, Helsinki, Finland
| |
Collapse
|
20
|
Riccardi E, Lervik A, Roet S, Aarøen O, Erp TS. PyRETIS 2: An improbability drive for rare events. J Comput Chem 2019; 41:370-377. [DOI: 10.1002/jcc.26112] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Revised: 10/25/2019] [Accepted: 10/29/2019] [Indexed: 01/27/2023]
Affiliation(s)
- Enrico Riccardi
- Department of ChemistryNorwegian University of Science and Technology Høgskoleringen 5 7491 Trondheim Norway
| | - Anders Lervik
- Department of ChemistryNorwegian University of Science and Technology Høgskoleringen 5 7491 Trondheim Norway
| | - Sander Roet
- Department of ChemistryNorwegian University of Science and Technology Høgskoleringen 5 7491 Trondheim Norway
| | - Ola Aarøen
- Department of Biotechnology and Food ScienceNorwegian University of Science and Technology Høgskoleringen 5 7491 Trondheim Norway
| | - Titus S. Erp
- Department of ChemistryNorwegian University of Science and Technology Høgskoleringen 5 7491 Trondheim Norway
- Center for Molecular Modeling (CMM)Ghent University Technologiepark 903 9052 Zwijnaarde Belgium
| |
Collapse
|
21
|
Affiliation(s)
- Zeke A. Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Oluwaseun O. Mesele
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H. Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
22
|
Cabriolu R, Skjelbred Refsnes KM, Bolhuis PG, van Erp TS. Foundations and latest advances in replica exchange transition interface sampling. J Chem Phys 2017; 147:152722. [DOI: 10.1063/1.4989844] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Raffaela Cabriolu
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, 7491 Trondheim, Norway
| | - Kristin M. Skjelbred Refsnes
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, 7491 Trondheim, Norway
| | - Peter G. Bolhuis
- Van ’t Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| | - Titus S. van Erp
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, 7491 Trondheim, Norway
| |
Collapse
|
23
|
Riccardi E, Dahlen O, van Erp TS. Fast Decorrelating Monte Carlo Moves for Efficient Path Sampling. J Phys Chem Lett 2017; 8:4456-4460. [PMID: 28857565 DOI: 10.1021/acs.jpclett.7b01617] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Many relevant processes in chemistry, physics, and biology are rare events from a computational perspective as they take place beyond the accessible time scale of molecular dynamics (MD). Examples are chemical reactions, nucleation, and conformational changes of biomolecules. Path sampling is an approach to break this time scale limit via a Monte Carlo (MC) sampling of MD trajectories. Still, many trajectories are needed for accurately predicting rate constants. To improve the speed of convergence, we propose two new MC moves, stone skipping and web throwing. In these moves, trajectories are constructed via a sequence of subpaths obeying superdetailed balance. By a reweighting procedure, almost all paths can be accepted. Whereas the generation of a single trajectory becomes more expensive, the reduced correlation results in a significant speedup. For a study on DNA denaturation, the increase was found to be a factor 12.
Collapse
Affiliation(s)
- Enrico Riccardi
- Department of Chemistry, Norwegian University of Science and Technology , NO-7491 Trondheim, Norway
| | - Oda Dahlen
- Department of Chemistry, Norwegian University of Science and Technology , NO-7491 Trondheim, Norway
| | - Titus S van Erp
- Department of Chemistry, Norwegian University of Science and Technology , NO-7491 Trondheim, Norway
| |
Collapse
|