Absorption and reactive properties of fluids in porous media are key to the design and improvement of numerous energy related applications. Molecular simulations of these systems require accurate force fields that capture the involved chemical reactions and have the ability to describe the vapor-liquid equilibrium (VLE). Two new reactive force fields (ReaxFF) for CO2 and H2O are developed, which are capable of not only modeling bond breaking and formation in reactive environments but also predicting their VLEs at saturation conditions. These new force fields include extra terms (ReaxFF-lg) to improve the long-range interactions between the molecules. For validation, we have developed a new Gibbs ensemble Monte Carlo (GEMC-ReaxFF) approach to predict the VLE. Computed VLE data show good agreement with National Institute of Standards and Technology reference data as well as existing nonreactive force fields. This validation proves the applicability of the GEMC-ReaxFF method to test new reactive force fields, and simultaneously it proves the applicability to extend newly developed ReaxFF force fields to other more complex reactive systems.

Thermochemical heat-storage applications, based on the reversible endo-/exothermic hydration reaction of salts, are intensively investigated to search for compact heat-storage devices. To achieve a truly valuable storage system, progressively complex salts are investigated. For these salts, the equilibrium temperature and pressure conditions are not always easy to predict. However, these conditions are crucial for the design of thermochemical heat-storage systems. A biased grand-canonical Monte Carlo (GCMC) tool is developed, enabling the study of equilibrium conditions at the molecular level. The GCMC algorithm is combined with reactive force field molecular dynamics (ReaxFF), which allows bond formation within the simulation. The Weeks-Chandler-Andersen (WCA) potential is used to scan multiple trial positions for the GCMC algorithm at a small cost. The most promising trial positions can be selected for recomputation with the more expensive ReaxFF. The developed WCA-ReaxFF-GCMC tool was used to study the hydration of MgCl2·nH2O. The simulation results show a good agreement with experimental and thermodynamic equilibriums for multiple hydration levels. The hydration shows that water, present at the surface of crystalline salt, deforms the surface layers and promotes further hydration of these deformed layers. Additionally, the WCA-ReaxFF-GCMC algorithm can be used to study other, non-TCM-related, reactive sorption processes.

The Diels-Alder cycloaddition (DAC) is a powerful tool to construct C-C bonds. The DAC reaction can be accelerated in several ways, one of which is reactant confinement as observed in supramolecular complexes and Diels-Alderases. Another method is altering the frontier molecular orbitals (FMOs) of the reactants by using homogeneous transition-metal complexes whose active sites exhibit d-orbitals suitable for net-bonding orbital interactions with the substrates. Both features can be combined in first row d-block (TM) exchanged faujasite catalysts where the zeolite framework acts as a stabilizing ligand for the active site while confining the reactants. Herein, we report on a mechanistic and periodic DFT study on TM-(Cu(I), Cu(II), Zn(II), Ni(II), Cr(III), Sc(III), V(V))exchanged faujasites to elucidate the effect of d-shell filling on the DAC reaction between 2,5-dimethylfuran and ethylene. Two pathways were found: one being the concerted one-step and the other being the stepwise two-step pathway. A decrease in d-shell filling results in a concomitant increase in reactant activation as evidenced by increasingly narrow energy gaps and lower activation barriers. For models holding relatively small d-block cations, the zeolite framework was found to bias the DAC reaction toward an asynchronous one-step pathway instead of the two-step pathway. This work is an example of how the active site properties and the surrounding chemical environment influence the reaction mechanism of chemical transformations.

The Diels-Alder cycloaddition (DAC) reaction is a commonly employed reaction for the formation of C-C bonds. DAC catalysis can be achieved by using Lewis acids and via reactant confinement in aqueous nanocages. Low-silica alkali-exchanged faujasite catalysts combine these two factors in one material. They can be used in the tandem DAC/dehydration reaction of biomass-derived 2,5-dimethylfuran (DMF) with ethylene toward p-xylene, in which the DAC reaction step initiates the overall reaction cycle. In this work, we performed periodic density functional theory (DFT) calculations on the DAC reaction between DMF and C₂H₄ in low-silica alkali(M)-exchanged faujasites (MY; Si/Al = 2.4; M = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺). The aim was to investigate how confinement of reactants in MY catalysts changed their electronic structure and the DAC-reactivity trend among the evaluated MY zeolites. The conventional high-silica alkali-exchanged isolated site model (MFAU; Si/Al = 47) served as a reference. The results show that confinement leads to initial-state (IS) destabilization and transition-state (TS) stabilization. Among the tested MY, most significant IS destabilization is found in RbY. Only antibonding orbital interactions between the reactants/reactive complex and cations were found, indicating that TS stabilization arises from ionic interactions. Additionally, in RbY the geometry of the transition state is geometrically most similar to that of the initial and final state. RbY also exhibits an optimal combination of the confinement-effects, resulting in having the lowest computed DAC-activation energy. The overall effect is a DAC-reactivity trend inversion in MY as compared to the trend found in MFAU where the activation energy correlates with the Lewis acidity of the exchangeable cations.

The one-pot Diels-Alder cycloaddition (DAC)/dehydration (D) tandem reaction between 2,5-dimethylfuran and ethylene is a potent pathway toward biomass-derived p-xylene. In this work, we present a cheap and active low-silica potassium-exchanged faujasite (KY, Si/Al = 2.6) catalyst. Catalyst optimization was guided by a computational study of the DAC/D reaction mechanism over different alkali-exchanged faujasites using periodic density functional theory calculations complemented by microkinetic modeling. Two types of faujasite models were compared, i.e., a high-silica alkali-exchanged faujasite model representing isolated active cation sites and a low-silica alkali-exchanged faujasite in which the reaction involves several cations in the proximity. The mechanistic study points to a significant synergetic cooperative effect of the ensemble of cations in the faujasite supercage on the DAC/D reaction. Alignment of the reactants by their interactions with the cationic sites and stabilization of reaction intermediates contribute to the high catalytic performance. Experiments confirmed the prediction that KY is the most active catalyst among low-silica alkali-exchanged faujasites. This work is an example of how the catalytic reactivity of zeolites depends on multiple interactions between the zeolite and reagents.