We present RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. In this article, we introduce the new additions and changes compared to RASPA2. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub. Using this code, we compare four Monte Carlo (MC) insertion/deletion techniques: unbiased Metropolis MC, Configurational-Bias Monte Carlo (CBMC), Continuous Fractional Component MC (CFCMC), and CB/CFCMC. We compare particle distribution shapes, acceptance ratios, accuracy and speed of isotherm computation, enthalpies of adsorption, and chemical potentials, over a wide range of loadings and systems, for the grand canonical ensemble and for the Gibbs ensemble.

@en

Entropies for alkane isomers longer than C10 are computed using our recently developed linear regression model for thermochemical properties which is based on second-order group contributions. The computed entropies show excellent agreement with experimental data and data from Scott’s tables which are obtained from a statistical mechanics-based correlation. Entropy production and heat input are calculated for the hydroisomerization of C7 isomers in various zeolites (FAU-, ITQ-29-, BEA-, MEL-, MFI-, MTW-, and MRE-types) at 500 K at chemical equilibrium. Small variations in these properties are observed because of the differences in reaction equilibrium distributions for these zeolites. The effect of chain length on heat input and entropy production is also studied for the hydroisomerization of C7, C8, C10, and C14 isomers in MTW-type zeolite at 500 K. For longer chains, both heat input and entropy production increase. Enthalpies and absolute entropies of C7 hydroisomerization reaction products in MTW-type zeolite increase with higher temperatures. These findings highlight the accuracy of our linear regression model in computing entropies for alkanes and provide insight for designing and optimizing zeolite-catalyzed hydroisomerization processes.@en

Linear regression (LR) is used to predict thermochemical properties of alkanes at temperatures (0–1000) K to study chemical reaction equilibria inside zeolites. The thermochemical properties of C1 until C10 isomers reported by Scott are used as training data sets in the LR model which is used to predict these properties for alkanes longer than C10 isomers. Second-order groups are used as independent variables which account for the interactions between the neighboring groups of atoms. This model accurately predicts Gibbs free energies, enthalpies, Gibbs free energies of formation, and enthalpies of formation for alkanes which exceeds the chemical accuracy of 1 kcal/mol and outperforms the group contribution methods developed by Benson et al., Joback and Reid, and Constantinou and Gani. Predictions from our model are used to compute the reaction equilibrium distribution of hydroisomerization of C10 and C14 isomers in MTW-type zeolite. Calculation of reaction equilibrium distribution inside zeolites also requires Henry coefficients of the isomers which can be computed using classical force field-based molecular simulations using the RASPA2 software for which we created an automated workflow. The reaction equilibrium distribution for C10 isomers obtained using the LR model and the training data set for this model are in very good agreement. The tools developed in this study will enable the computational study of hydroisomerization of long-chain alkanes (>C10).@en

We study important aspects of shape selectivity effects of zeolites for hydroisomerization of linear alkanes, which produces a myriad of isomers, particularly for long chain hydrocarbons. To investigate the conditions for achieving an optimal yield of branched hydrocarbons, it is important to understand the role of chemical equilibrium in these reversible reactions. We conduct an extensive analysis of shape selectivity effects of different zeolites for the hydroisomerization of C₇ and C₈ isomers at chemical reaction equilibrium conditions. The reaction ensemble Monte Carlo method, coupled with grand-canonical Monte Carlo simulations, is commonly used for computing reaction equilibrium of heterogeneous reactions. The computational demands become prohibitive for a large number of reactions. We used a faster alternative in which reaction equilibrium is obtained by imposing chemical equilibrium in the gas phase and phase equilibrium between the gas phase components and the adsorbed phase counterparts. This effectively mimics the chemical equilibrium distribution in the adsorbed phase. Using Henry’s law at infinite dilution and mixture adsorption isotherm models at elevated pressures, we calculate the adsorbed loadings in the zeolites. This study shows that zeolites with cage or channel-like structures exhibit significant differences in selectivity for alkane isomers. We also observe a minimal impact of pressure on the gas-phase equilibrium of these reactions at typical experimental reaction temperatures 400 − 700 K . This study marks initial strides in understanding the reaction product distribution for long-chain alkanes.

@en

Ideal Adsorbed Solution Theory (IAST) is a common method for modelling mixture adsorption isotherms based on pure component isotherms. When the adsorbent has distinct adsorption sites, the segregated version of IAST (SIAST) provides improved adsorbed loadings compared to IAST. We have adopted the concept of SIAST and applied it to an explicit isotherm model which takes into account the different sizes of the adsorbates: the so called Segregated Explicit Isotherm (SEI). The purpose of SEI is to have an explicit adsorption model that can consider both size-effects of the co-adsorbed molecules and surface heterogeneities. In sharp contrast to IAST and SIAST, no iterative scheme is required in case of SEI, which leads to much faster simulations. A comparative study has been performed to analyse the adsorption isotherms calculated using these three methods. The adsorbed loadings predicted by SEI and SIAST are in excellent agreement with the Grand-Canonical Monte Carlo (GCMC) simulation data. The loadings estimated by IAST show considerable deviations from the GCMC data at high pressures. Breakthrough curve modelling is used to compare the effects of these three models at dynamic conditions. The explicit model (SEI) leads to the fastest simulation run time, followed by SIAST.

@en

We present the RUPTURA code (https://github.com/iraspa/ruptura) as a free and open-source software package (MIT license) for (1) the simulation of gas adsorption breakthrough curves, (2) mixture prediction using methods like the Ideal Adsorption Solution Theory (IAST), segregated-IAST and explicit isotherm models, and (3) fitting of isotherm models on computed or measured adsorption isotherm data. The combination with the RASPA software enables computation of breakthrough curves directly from adsorption simulations in the grand-canonical ensemble. RUPTURA and RASPA have similar input styles. IAST is implemented near machine precision but we also provide several explicit mixture prediction methods that are non-iterative and potentially faster than IAST. The code supports a wide variety of isotherm models like Langmuir, Anti-Langmuir, BET, Henry, Freundlich, Sips, Langmuir-Freundlich, Redlich-Peterson, Toth, Unilan, O'Brian & Myers, Asymptotic Temkin, and Bingel & Walton. The isotherm model parameters can easily be obtained by the fitting module. Breakthrough plots and animations of the column properties are automatically generated. In addition to highlighting the code, we also review all the developed techniques from literature for mixture prediction, breakthrough simulations, and isotherm model fitting, and provide a tutorial discussing the workflows.

@en

A reactor model that deconvolutes thermodynamics of adsorption of hydrocarbon in the pores of zeolite Beta, obtained by Configurational-bias Monte Carlo simulations, from intrinsic, intraporous kinetics of hydroisomerization and hydrocracking reactions, provides a good quantitative description of all significant reactions in the kinetic network for interconversion and cracking of different heptane isomers. Activation enthalpies obtained for intraporous reactions follow the expected order according to the carbenium ion formalism: methyl shift< ethyl shift < isom(B) ∼ crack(B2) < crack(B1) < crack(C) ∼ crack(D) < crack(E) and apparently within each isomerization class, in terms of carbenium ions formally involved: sec → tert < sec → sec ∼ tert → tert < tert → sec. except for the ethyl shift reaction forming 3-ethylpentane. Cracking happens primarily through 2,4-dimethylpentane (type B2), regardless of the initial reactant. The model can be subsequently used to separate the effect of pore structure on selective adsorption and on intraporous reaction kinetics. Zeolite Beta will serve as a base case for a comparison of different zeolite structures.

@en