Reaction Ensemble Monte Carlo Simulation of Xylene Isomerization in Bulk Phases and under Confinement

Competitive Adsorption of Xylenes at Chemical Equilibrium in Zeolites

The separation of xylenes is one of the most important processes in the petrochemical industry. In this article, the competitive adsorption from a fluid-phase mixture of xylenes in zeolites is studied. Adsorption from both vapor and liquid phases is considered. Computations of adsorption of pure xylenes and a mixture of xylenes at chemical equilibrium in several zeolite types at 250 °C are performed by Monte Carlo simulations. It is observed that shape and size selectivity entropic effects are predominant for small one-dimensional systems. Entropic effects due to the efficient arrangement of xylenes become relevant for large one-dimensional systems. For zeolites with two intersecting channels, the selectivity is determined by a competition between enthalpic and entropic effects. Such effects are related to the orientation of the methyl groups of the xylenes. m-Xylene is preferentially adsorbed if xylenes fit tightly in the intersection of the channels. If the intersection is much larger than the adsorbed molecules, p-xylene is preferentially adsorbed. This study provides insight into how the zeolite topology can influence the competitive adsorption and selectivity of xylenes at reaction conditions. Different selectivities are observed when a vapor phase is adsorbed compared to the adsorption from a liquid phase. These insight have a direct impact on the design criteria for future applications of zeolites in the industry. MRE-type and AFI-type zeolites exclusively adsorb p-xylene and o-xylene from the mixture of xylenes in the liquid phase, respectively. These zeolite types show potential to be used as high-performing molecular sieves for xylene separation and catalysis.

The effects of confinement in pores built of folded graphene sheets on the equilibrium of nitrogen monoxide dimerisation reaction

In the current work we have used reactive Monte Carlo simulations to systematically study the effects of graphene folding on equilibria of NO dimerisation occurring at isolated surfaces and in porous networks built of corrugated graphene sheets. It has been demonstrated that the folding of isolated graphene sheets significantly improves the yield of reactions occurring on their surface. Then, it has also been shown that in slit-like pores formed by the folded graphene sheets the reaction yield depends on the corrugation and arrangement of the pore walls. It has been found that the reaction yield increases when the walls' corrugation is high because of the appearance of narrow regions and/or wedge-like regions in the pores. The condensation of reacting fluid in such places, where the bulges at both walls are close one to another, leads to much higher reaction yield than on the surface of isolated sheets. Thus, we recommended the highly corrugated graphene to control the chemical reactions.

Carbon Nanohorns as Reaction Nanochambers - a Systematic Monte Carlo Study

Carbon nanohorns (CNHs, one of the newest carbon allotropes) have been subjected to intensive experimental and theoretical studies due to their potential applications. One of such applications can be their use as reaction nanochambers. However, experimental studies on the reaction equilibria under confinement are extremely challenging since accurate measurements of the concentrations of reacting species in pores are a very hard task. So, the main ways to examine such phenomena are theoretical methods (e.g. the reactive Monte Carlo, RxMC). We have presented the first systematic RxMC study on the influence of the CNH’s geometric parameters (the apex angle, the diameter, and the length) on reaction equilibria, taking the nitrogen monoxide dimerisation as an example. All the investigated parameters significantly affect the reaction yield at low and moderate coverages. Short and narrow CNHs have been found to be preferred. However, the key factor influencing the reaction equilibria is the presence of a conical part. Energetics of interactions between the reacting molecules in this fragment of a nanohorn maximises the effects of confinement. In consequence, CNHs have the advantage over their nanotube counterparts of the same diameter. The obtained results have confirmed that CNHs can be considered as potential reaction nanochambers.

Molecular Simulation of Chemical Reaction Equilibrium by Computationally Efficient Free Energy Minimization

The molecular simulation of chemical reaction equilibrium (CRE) is a challenging and important problem of broad applicability in chemistry and chemical engineering. The primary molecular-based approach for solving this problem has been the reaction ensemble Monte Carlo (REMC) algorithm [Turner et al. Molec. Simulation2008, 34, (2), 119-146], based on classical force-field methodology. In spite of the vast improvements in computer hardware and software since its original development almost 25 years ago, its more widespread application is impeded by its computational inefficiency. A fundamental problem is that its MC basis inhibits the implementation of significant parallelization, and its successful implementation often requires system-specific tailoring and the incorporation of special MC approaches such as replica exchange, expanded ensemble, umbrella sampling, configurational bias, and continuous fractional component methodologies. We describe herein a novel CRE algorithm (reaction ensemble molecular dynamics, ReMD) that exploits modern computer hardware and software capabilities, and which can be straightforwardly implemented for systems of arbitrary size and complexity by exploiting the parallel computing methodology incorporated within many MD software packages (herein, we use GROMACS for illustrative purposes). The ReMD algorithm utilizes these features in the context of a macroscopically inspired and generally applicable free energy minimization approach based on the iterative approximation of the system Gibbs free energy function by a mathematically simple convex ideal solution model using the composition at each iteration as a reference state. Finally, we additionally describe a simple and computationally efficient a posteriori method to estimate the equilibrium concentrations of species present in very small amounts relative to others in the primary calculation. To demonstrate the algorithm, we show its application to two classic example systems considered previously in the literature: the N₂-O₂-NO system and the ammonia synthesis system.

Reaction Ensemble Monte Carlo Simulations of CO2 Absorption in the Reactive Ionic Liquid Triethyl(octyl)phosphonium 2-Cyanopyrrolide

The absorption of CO₂ into an aprotic heterocyclic anion ionic liquid (IL) is modeled using reaction ensemble Monte Carlo (RxMC) with the semigrand reaction move. RxMC has previously been unable to sample chemical equilibrium involving molecular ions in nanostructured liquids due to the high free-energy requirements to open and close cavities and restructure the surrounding environment. Our results are validated by experiments in the modeled IL, triethyl(octyl)phosphonium 2-cyanopyrrolide ([P₂₂₂₈][cnp]), and in a close analog with longer alkyl chains on the cation. Heats of absorption and reaction from both experiment and simulation are exothermic and of comparable magnitude. Replacing experimental Henry's constants with their simulated counterparts improves the accuracy of a Langmuir-type model at moderate pressures. Nonidealities that affect chemical equilibrium are identified and calculated with high precision.

Alchemical free energy calculations and umbrella sampling with local molecular field theory

Understanding the thermodynamic driving forces underlying any chemical process requires a description of the underlying free energy surface. However, computation of free energies is difficult, often requiring advanced sampling techniques. Moreover, these computations can be further complicated by the evaluation of any long-ranged interactions in the system of interest, such as Coulomb interactions in charged and polar media. Local molecular field theory is a promising approach to avoid many of the conceptual and computational difficulties associated with long-ranged interactions. We present frameworks for performing alchemical free energy calculations and non-Boltzmann sampling with local molecular field theory. We demonstrate that local molecular field theory can be used to perform these free energy calculations with accuracy comparable to traditional methodologies while eliminating the need for explicit treatment of long-ranged interactions in simulations.

Understanding the Reactive Adsorption of H2 S and CO2 in Sodium-Exchanged Zeolites

Purifying sour natural gas streams containing hydrogen sulfide and carbon dioxide has been a long-standing environmental and economic challenge. In the presence of cation-exchanged zeolites, these two acid gases can react to form carbonyl sulfide and water (H₂ S+CO₂ ⇌H₂ O+COS), but this reaction is rarely accounted for. In this work, we carry out reactive first-principles Monte Carlo (RxFPMC) simulations for mixtures of H₂ S and CO₂ in all-silica and Na-exchanged forms of zeolite beta to understand the governing principles driving the enhanced conversion. The RxFPMC simulations show that the presence of Na⁺ cations can change the equilibrium constant by several orders of magnitude compared to the gas phase or in all-silica beta. The shift in the reaction equilibrium is caused by very strong interactions of H₂ O with Na⁺ that reduce the reaction enthalpy by about 20 kJ mol^-1 . The simulations also demonstrate that the siting of Al atoms in the framework plays an important role. The RxFPMC method presented here is applicable to any chemical conversion in any confined environment, where strong interactions of guest molecules with the host framework and high activation energies limit the use of other computational approaches to study reaction equilibria.

Adsorption equilibrium of nitrogen dioxide in porous materials

The effect of confinement on the equilibrium reactive system containing nitrogen dioxide and dinitrogen tetroxide is studied by molecular simulation and the reactive Monte Carlo (RxMC) approach. The bulk-phase reaction was successfully reproduced and five all-silica zeolites (i.e. FAU, FER, MFI, MOR, and TON) with different topologies were selected to study their adoption behavior. Dinitrogen tetroxide showed a stronger affinity than nitrogen dioxide in all the zeolites due to size effects, but exclusive adsorption sites in MOR allowed the adsorption of nitrogen dioxide with no competition at these sites. From the study of the adsorption isotherms and isobars of the reacting mixture, confinement enhanced the formation of dimers over the full range of pressure and temperature, finding the largest deviations from bulk fractions at low temperature and high pressure. The channel size and shape of the zeolite have a noticeable influence on the dinitrogen tetroxide formation, being more important in MFI, closely followed by TON and MOR, and finally FER and FAU. Preferential adsorption sites in MOR lead to an unusually strong selective adsorption towards nitrogen dioxide, demonstrating that the topological structure has a crucial influence on the composition of the mixture and must be carefully considered in systems containing nitrogen dioxide.

Monte Carlo study of chemical reaction equilibria in pores of activated carbons

Systematic Monte Carlo studies concerning relationships between the porous structure of activated carbons and the equilibria of reactions under confinement are presented.