Group contribution methods (GCMs) provide a practical and computationally efficient approach for predicting thermodynamic properties of hydrocarbons, especially when experimental data are scarce. This review evaluates the evolution of GCMs from classical first-order schemes (e.g. Lydersen method, Joback method) to more advanced second-order frameworks (e.g. CG94, Sharma method), hybrid extensions, and emerging machine learning integrations. While first-order models are simple and widely used, these models struggle with branched and long-chain molecules. Second-order approaches significantly improve structural sensitivity and predictive accuracy, achieving deviations below 2–3% for critical properties and within 1 kcal/mol for formation enthalpies of branched alkanes. Nevertheless, challenges remain in extrapolating to highly complex molecules, underrepresented functional groups, and extreme conditions. Promising directions include reinforcement of second-order GCMs with molecular theory, systematic expansion of experimental and quantum-based datasets, and hybrid GCM–machine learning models that retain interpretability while improving generalisability. We recommend prioritising models that balance accuracy, robustness, simplicity, and transferability to accelerate sustainable process and product designs, particularly in applications such as fuel upgrading including hydroisomerisation, separation processes, and green chemical development.

Shape-selective adsorption in zeolites plays a pivotal role in catalytic hydroisomerization of long-chain alkanes, a key process in producing sustainable aviation fuels from Fischer–Tropsch products. Accurately predicting adsorption behavior for the large number of alkane isomers in different zeolite frameworks is computationally intensive. To address this, we have developed a machine learning framework that rapidly and accurately predicts Henry coefficients of linear (C1–C30) and branched (C4–C20) alkanes in one-dimensional zeolites. Using descriptors based on chain length, branching patterns, and molecular graphs, we evaluate multiple ML models, including Random Forest, XGBoost, CatBoost, TabPFN, and D-MPNN in MTT-, MTW-, MRE-, and AFI-type zeolites. TabPFN and D-MPNN offer the highest predictive accuracy. Active learning further boosts model performance by efficiently selecting diverse and structurally informative isomers. We also uncover activity cliffs, where small changes in molecular structure lead to sharp variations in adsorption, and demonstrate that targeted oversampling of these cases improves model robustness. Finally, we combine the ML-predicted Henry coefficients with gas-phase thermodynamics to compute reaction equilibrium distributions for C16 hydroisomerization. This integrated, data-driven approach enables efficient screening and design of shape-selective zeolite catalysts, thereby reducing the need for costly simulations

Accurate prediction of thermodynamic properties of hydrocarbons is essential for chemical process modelling. Conventional group contribution methods often are used to predict these properties. However, these methods often require extensive parameter sets to handle structural complexities. A refined group contribution method for predicting thermodynamic properties of hydrocarbon isomers with reduced complexity and improved accuracy is presented and discussed. By combining the structural framework of Constantinou and Gani (CG94) with a sensitivity-based selection of second-order groups, a reduced yet highly effective set of twelve second-order groups is identified. This reduced set retains the predictive power comparable to more complex models while significantly reducing the number of parameters. Linear regression is applied to model enthalpies and Gibbs free energies of formation for a wide temperature range. To test broader applicability, the model is further extended to properties that require nonlinear regression, including critical temperatures, critical pressures, acentric factors, and liquid densities. For all cases, the proposed model achieves high predictive accuracy, demonstrating its robustness and generalizability. This methodology balances interpretability, efficiency, and performance, making it suitable for both research and industrial thermodynamic modelling.

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.

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.

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.

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.