Insight into the role of excess hydroxide ions in silicate condensation reactions.

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.

Impact of Organic Templates on the Selective Formation of Zeolite Oligomers.

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.

Thiosquaramide-Based Supramolecular Polymers: Aromaticity Gain in a Switched Mode of Self-Assembly.

Despite a growing understanding of factors that drive monomer self-assembly to form supramolecular polymers, the effects of aromaticity gain have been largely ignored. Herein, we document the aromaticity gain in two different self-assembly modes of squaramide-based bolaamphiphiles. Importantly, O → S substitution in squaramide synthons resulted in supramolecular polymers with increased fiber flexibility and lower degrees of polymerization. Computations and spectroscopic experiments suggest that the oxo- and thiosquaramide bolaamphiphiles self-assemble into "head-to-tail" versus "stacked" arrangements, respectively. Computed energetic and magnetic criteria of aromaticity reveal that both modes of self-assembly increase the aromatic character of the squaramide synthons, giving rise to stronger intermolecular interactions in the resultant supramolecular polymer structures. These examples suggest that both hydrogen-bonding and stacking interactions can result in increased aromaticity upon self-assembly, highlighting its relevance in monomer design.

Diffusion of gas mixtures in the sI hydrate structure.

Replacing methane with carbon dioxide in gas hydrates has been suggested as a way of harvesting methane, while at the same time storing carbon dioxide. Experimental evidence suggests that this process is facilitated if gas mixtures are used instead of pure carbon dioxide. We studied the free energy barriers for diffusion of methane, carbon dioxide, nitrogen, and hydrogen in the sI hydrate structure using molecular simulation techniques. Cage hops between neighboring cages were considered with and without a water vacancy and with a potential inclusion of an additional gas molecule in either the initial or final cage. Our results give little evidence for enhanced methane and carbon dioxide diffusion if nitrogen is present as well. However, the inclusion of hydrogen seems to have a substantial effect as it diffuses rapidly and can easily enter occupied cages, which reduces the barriers of diffusion for the gas molecules that co-occupy a cage with hydrogen.

Geometrical flexibility of platinum nanoclusters: impacts on catalytic decomposition of ethylene glycol.

Catalytic decomposition of ethylene glycol on the Pt13 cluster was studied as a model system for hydrogen production from a lignocellulosic material. Ethylene glycol was chosen as a starting material because of two reasons, it is the smallest oxygenate with a 1 : 1 carbon to oxygen ratio and it contains the C-H, O-H, C-C, and C-O bonds also present in biomass. Density functional theory calculations were employed for predictions of reaction pathways for C-H, O-H, C-C and C-O cleavages, and Brønsted-Evans-Polanyi relationships were established between the final state and the transition state for all mechanisms. The results show that Pt13 catalyzes the cleavage reactions of ethylene glycol more favourably than a Pt surface. The flexibility of Pt13 clusters during the reactions is the key factor in reducing the activation barrier. Overall, the results demonstrate that ethylene glycol and thus biomass can be efficiently converted into hydrogen using platinum nanoclusters as catalysts.

Rare event simulations reveal subtle key steps in aqueous silicate condensation.

A replica exchange transition interface sampling (RETIS) study combined with Born-Oppenheimer molecular dynamics (BOMD) is used to investigate the dynamics, thermodynamics and the mechanism of the early stages of the silicate condensation process. In this process, two silicate monomers, of which one is an anionic species, form a negatively charged five-coordinated silicate dimer. In a second stage, this dimer can fall apart again, forming the original monomers, or release a water molecule into the solution. We studied the association and dissociation reaction in the gas phase, and the dissociation and water removal step in the aqueous phase. The results on the aqueous phase dissociation suggest two possible mechanisms. The breakage of the bond between the intermediate oxygen and the five-coordinated silicon is sometimes accompanied by a proton transfer. After dissociation into silicate monomers, the anionic monomer is either the previously four-coordinated silicon or the previously five-coordinated silicon depending on whether the hydrogen transfer occurs or not. Our results show that the mechanism of proton transfer is highly predominant. Water removal simulations also show two possible mechanisms distinguished by the proton transfer reaction path. Proton transfer can occur either via a direct or via a water mediated reaction step. The calculations reveal that although both mechanisms contribute to the water removal process, the direct proton transfer is slightly favorable and occurs roughly in six out of ten occasions.

Chemically accurate energy barriers of small gas molecules moving through hexagonal water rings.

We present chemically accurate potential energy curves of CH4, CO2 and H2 moving through hexagonal water rings, calculated by CCSD(T)/aug-cc-pVTZ with counterpoise correction. The barriers are extracted from a potential energy surface obtained by allowing the water ring to expand while the gas molecule diffuses through. State-of-the-art XC-functionals are evaluated against the CCSD(T) potential energy surface.

The mechanism of the initial step of germanosilicate formation in solution: a first-principles molecular dynamics study.

The condensation reactions between Ge(OH)4 and Si(OH)4 units in solution are studied to understand the mechanism and stable species during the initial steps of the formation process of Ge containing zeolites under basic conditions. The free energy of formation of (OH)3Ge-O-Ge-(OH)2O(-), (OH)3Si-O-Si-(OH)2O(-), (OH)3Ge-O-Si-(OH)2O(-) and (OH)3Si-O-Ge-(OH)2O(-) dimers is calculated with ab initio molecular dynamics and thermodynamic integration, including an explicit description of the water solvent molecules. Calculations show that the attack of the conjugated base (Ge(OH)3O(-) and Si(OH)3O(-)) proceeds with a smaller barrier at the Ge center. In addition, the formation of the pure germanate dimer is more favorable than that of the germano-silicate structure. These results explain the experimental observation of Ge-Ge and Si-Ge dimer species in solutions, with a few Si-Si ones.

Heat transport through a solid-solid junction: the interface as an autonomous thermodynamic system.

We perform computational experiments using nonequilibrium molecular dynamics simulations, showing that the interface between two solid materials can be described as an autonomous thermodynamic system. We verify the local equilibrium and give support to the Gibbs description of the interface also away from the global equilibrium. In doing so, we reconcile the common formulation of the thermal boundary resistance as the ratio between the temperature discontinuity at the interface and the heat flux with a more rigorous derivation from nonequilibrium thermodynamics. We also show that thermal boundary resistance of a junction between two pure solid materials can be regarded as an interface property, depending solely on the interface temperature, as implicitly assumed in some widely used continuum models, such as the acoustic mismatch model. Thermal rectification can be understood on the basis of different interface temperatures for the two flow directions.

Note: A new truncation correction for the configurational temperature extends its applicability to interaction potentials with a discontinuous force.

We present a simple truncation correction for the configurational temperature which, unlike previous corrections, works even at low truncation values for the shifted and truncated Lennard-Jones potential. The success of the new correction suggests that the expression for the configurational temperature is valid also for interaction potentials with a discontinuous force, given that the discontinuity is properly accounted for.

Heat and mass transfer across interfaces in complex nanogeometries.

Heat and mass transfer in nanodevices depends much on the geometry due to the strong influence of curvature on interfacial properties, such as the Kapitza resistance. We present a method which combines nonequilibrium square gradient theory and nonequilibrium molecular dynamics simulations to obtain the coefficients in a curvature expansion of the interface transfer coefficients. The expansion can be used directly to describe heat and mass transfer in complex nanogeometries. As examples of complex nanogeometries, we consider an oblate spheroidal droplet, a prolate spheroidal bubble, and a toroidal bubble. Depending on the sign and magnitude of the curvature, transfer is enhanced or reduced significantly. The presented method is applicable to many types of interfaces and substances, and we expect it to contribute to the understanding and design of future nanodevices.

Low barriers for hydrogen diffusion in sII clathrate.

The transport of gas molecules in hydrates is presently poorly understood. In sII structured hydrates with hydrogen guests there is, for instance, a mismatch between experimental and computed values for diffusion constants. We provide an explanation for the experimentally observed diffusion rates, using DFT-based molecular dynamics simulations at 100 K. By considering the effect of cage occupancy, as well as the flexibility of the water lattice, we show that barriers for hydrogen diffusing between cages, can approach values as low as 5 kJ mol(-1), which is very close to experimental values.

The role of a structure directing agent tetramethylammonium template in the initial steps of silicate oligomerization in aqueous solution.

The understanding of the formation of silicate oligomers in the initial stage of zeolite synthesis is of fundamental scientific and technological importance. The use of different organic structure directing agents is known to be a key factor in the formation of different silicate species, and the final zeolite structure. Tetramethylammonium (TMA(+)), for example, is indispensable for the formation of the LTA zeolite type. However, the role of a TMA(+) template has not yet been elucidated at the molecular level. In this study, ab initio molecular dynamic simulations were combined with thermodynamic integration to arrive at an understanding of the role of TMA(+) in the formation of various silicate species, ranging from dimer to 4-ring. Free energy profiles show that trimer and 3-ring silicate are less favourable than other oligomers such as linear tetramer, branched tetramer and 4-ring structures. TMA(+) exhibits an important role in controlling the predominant species in solution via its coordination with silicate structures during the reaction process. This can explain that formation of D4R·8TMA crystals, as observed in experiment, is controlled by the single 4-ring formation step.

Density Functional Theory Study on the Interactions of Metal Ions with Long Chain Deprotonated Carboxylic Acids.

In this work, interactions between carboxylate ions and calcium or sodium ions are investigated via density functional theory (DFT). Despite the ubiquitous presence of these interactions in natural and industrial chemical processes, few DFT studies on these systems exist in the literature. Special focus has been placed on determining the influence of the multibody interactions (with up to 4 carboxylates and one metal ion) on an effective pair-interaction potential, such as those used in molecular mechanics (MM). Specifically, DFT calculations are employed to quantify an effective pair-potential that implicitly includes multibody interactions to construct potential energy curves for carboxylate-metal ion pairs. The DFT calculated potential curves are compared to a widely used molecular mechanics force field (OPLS-AA). The calculations indicate that multibody effects do influence the energetic behavior of these ionic pairs and the extent of this influence is determined by a balance between (a) charge transfer from the carboxylate to the metal ions which stabilizes the complex and (b) repulsion between carboxylates, which destabilizes the complex. Additionally, the potential curves of the complexes with 1 and 2 carboxylates and one counterion have been examined to higher separation distance (20 Å) by the use of relaxed scan optimization and constrained density functional theory (CDFT). The results from the relaxed scan optimization indicate that near the equilibrium distance, the charge transfer between the metal ion and the deprotonated carboxylic acid group is significant and leads to non-negligible differences between the DFT and MM potential curves, especially for calcium. However, at longer separation distances the MM calculated interaction potential functions converge to those calculated with CDFT, effectively indicating the approximate domain of the separation distance coordinate where charge transfer between the ions is occurring.

Finite-size and truncation effects for microscopic expressions for the temperature at equilibrium and nonequilibrium.

Several expressions have been proposed for the temperature in molecular simulations, where some of them have configurational contributions. We investigate how their accuracy is influenced by the number of particles in the simulation and the discontinuity in the derivatives of the interaction potential introduced by truncation. For equilibrium molecular dynamics with fixed total volume and fixed average total energy per particle, all the evaluated expressions including that for the kinetic temperature give a dependence on the total number of particles in the simulation. However, in a partitioned simulation volume under the same conditions, the mean temperature of each bin is independent of the number of bins. This finding is important for consistently defining a local temperature for use in nonequilibrium simulations. We identify the configurational temperature expressions which agree most with the kinetic temperature and find that they give close to identical results in nonequilibrium molecular dynamics (NEMD) simulations with a temperature gradient, for high and low density bulk-systems (both for transient and steady-state conditions), and across vapor-liquid interfaces, both at equilibrium and during NEMD simulations. The work shows that the configurational temperature is equivalent to the kinetic temperature in steady-state molecular dynamics simulations if the discontinuity in the derivatives of the interaction potential is handled properly, by using a sufficiently long truncation-distance or tail-corrections.

A test on reactive force fields for the study of silica dimerization reactions.

We studied silica dimerization reactions in the gas and aqueous phase by density functional theory (DFT) and reactive force fields based on two parameterizations of ReaxFF. For each method (both ReaxFF force fields and DFT), we performed constrained geometry optimizations, which were subsequently evaluated in single point energy calculations using the other two methods. Standard fitting procedures typically compare the force field energies and geometries with those from quantum mechanical data after a geometry optimization. The initial configurations for the force field optimization are usually the minimum energy structures of the ab initio database. Hence, the ab initio method dictates which structures are being examined and force field parameters are being adjusted in order to minimize the differences with the ab initio data. As a result, this approach will not exclude the possibility that the force field predicts stable geometries or low transition states which are realistically very high in energy and, therefore, never considered by the ab initio method. Our analysis reveals the existence of such unphysical geometries even at unreactive conditions where the distance between the reactants is large. To test the effect of these discrepancies, we launched molecular dynamics simulations using DFT and ReaxFF and observed spurious reactions for both ReaxFF force fields. Our results suggest that the standard procedures for parameter fitting need to be improved by a mutual comparative method.

Aromatic Gain in a Supramolecular Polymer.

The synergy of aromatic gain and hydrogen bonding in a supramolecular polymer is explored. Partially aromatic bis(squaramide) bolaamphiphiles were designed to self-assemble through a combination of hydrophobic, hydrogen-bonding, and aromatic effects into stiff, high-aspect-ratio fibers. UV and IR spectroscopy show electron delocalization and geometric changes within the squaramide ring indicative of strong hydrogen bonding and aromatic gain of the monomer units. The aromatic contribution to the interaction energy was further supported computationally by nucleus-independent chemical shift (NICS) and harmonic oscillator model of aromaticity (HOMA) indices, demonstrating greater aromatic character upon polymerization: at least 30% in a pentamer. The aromatic gain-hydrogen bonding synergy results in a significant increase in thermodynamic stability and a striking difference in aggregate morphology of the bis(squaramide) bolamphiphile compared to isosteres that cannot engage in this effect.

Thermal conductivity of carbon dioxide from non-equilibrium molecular dynamics: a systematic study of several common force fields.

We report a systematic investigation of the thermal conductivity of various three-site models of carbon dioxide (CO2) using nonequilibrium molecular dynamics in the temperature range 300-1000 K and for pressures up to 200 MPa. A direct comparison with experimental data is made. Three popular CO2 force fields (MSM, EPM2, and TraPPE) and two flexible models (based on EPM2) were investigated. All rigid force fields accurately predict the equation of state for carbon dioxide for the given range of variables. They can also reproduce the thermal conductivity of CO2 at room temperature and predict a decrease of the thermal conductivity with increasing temperature. At high temperatures, the rigid models underestimate the thermal conductivity.

Molecular insights into the interactions between PEG carriers and drug molecules from Celastrus hindsii: a multi-scale simulation study.

Efficient drug delivery is crucial for the creation of effective pharmaceutical treatments, and polyethylene glycol (PEG) carriers have been emerged as promising candidates for this purpose due to their bio-compatibility, enhancement of drug solubility, and stability. In this study, we utilized molecular simulations to examine the interactions between PEG carriers and selected drug molecules extracted from Celastrus hindsii: Hindsiilactone A, Hindsiiquinoflavan B, Maytenfolone A, and Celasdin B. The simulations provided detailed insights into the binding affinity, stability, and structural properties of these drug molecules when complexed with PEG carriers. A multi-scale approach combining density functional theory (DFT), extended tight-binding (xTB), and molecular dynamics (MD) simulations was conducted to investigate both unbound and bound states of PEG/drug systems. The results from DFT and xTB calculations revealed that the unbound complex has an unfavorable binding free energy, primarily due to negative contributions of delta solvation free energy and entropy. The MD simulations provided more detailed insights into the interactions between PEG and drug molecules in water solutions. By integrating the findings from the multi-scale simulations, a comprehensive picture of the unbound and bound states of PEG and drug systems were obtained. This information is valuable for understanding the molecular mechanisms governing the binding of drugs in PEG-based delivery platforms, and it contributes to the rational design and optimization of these systems.

Water adsorption on surfaces of calcium aluminosilicate crystal phase of stone wool: a DFT study.

Stone wool is widely used as an efficient thermal insulator within the construction industry; however, its performance can be significantly impacted by the presence of water vapor. By altering the material's characteristics and effective thermo-physical properties, water vapor can reduce overall efficacy in various environmental conditions. Therefore, understanding water adsorption on stone wool surfaces is crucial for optimizing insulation properties. Through the investigation of interaction between water molecules and calcium aluminosilicate (CAS) phase surfaces within stone wool using density functional theory (DFT), we can gain insight into underlying mechanisms governing water adsorption in these materials. This research aims to elucidate the molecular-level interaction between water molecules and CAS surfaces, which is essential for understanding fundamental properties that govern their adsorption process. Both dissociative and molecular adsorptions were investigated in this study. For molecular adsorption, the adsorption energy ranged from -  84 to -  113 kJ mol - 1 depending on surface orientation. A wider range of adsorption energy ( -  132 to -  236 kJ mol - 1 ) was observed for dissociative adsorption. Molecular adsorption was energetically favored on (010) surfaces while dissociative adsorption was most favorable on (111) surfaces. This DFT study provides valuable insights into the water adsorption behavior on low index surfaces of CAS phase in stone wool, which can be useful for designing effective strategies to manage moisture-related issues in construction materials. Based on these findings, additional research on the dynamics and kinetics of water adsorption and desorption processes of this thermal isolation material is suggested.