Atomistic Simulations of O2, N2, and CO2 Diffusion in Hydrochar─Quantifying the Effect of the Pore Size Regime
- 23 July 2026
- journal article
- Published by American Chemical Society (ACS) in Energy & Fuels
- Vol. 40 (31) , 16989-16998
Abstract
A rigorous yet flexible methodology was developed to generate realistic molecular char models that take experimental data on char geometry and elemental composition into account parametrically. Following this methodology, three upper-nanoscale models of an activated hydrochar derived from cellulose, with the composition C17708H5036O5962 were constructed, with densities from 0.90 to 1.15 gcm–3, porosities ranging from 20% to 36%, and maximum pore sizes between 20 and 30 Å. These models were employed to investigate the mobility of the light gases O2, N2, and CO2, which are typically present in oxy-fuel atmospheres confined within char nanopores and for which literature data are sparse and rather inconclusive. To achieve this, nonreactive molecular dynamics simulations were performed for several nanoseconds, using a modified formulation of the generic Universal Force Field. The studied combinations of char systems and gas molecules were investigated systematically with respect to gas molecule load and char pore size regime, encoded by the char density, porosity, and maximum pore diameter, to resolve the dependence of the diffusion dynamics of the confined gases on these parameters. Their interplay was quantified by computing self-diffusion coefficients for the confined gas species at 298.15 K, revealing absolute values in the 10–8 m2s–1 range for O2 and N2 and in the 10–9 m2s–1 range for CO2. The obtained values are comparable to the effective diffusion constants measured in previous experimental studies on related porous media, indicating that the models are sufficiently realistic within the chosen size regime. Furthermore, a linear decay of the gas diffusivities with the char density, linear growth with the char porosity, and superlinear growth with the maximum pore diameter were observed. Regression models fitted to the diffusion data may be incorporated into general diffusion equations, taking the pore size regime into account and supplementing future experimental studies on gas diffusion in char. The employed computational workflow can be used to characterize char and therein confined gases under experimentally inaccessible conditions.Keywords
Funding Information
- Deutsche Forschungsgemeinschaft (215035359)
This publication has 68 references indexed in Scilit:
- Atomistic Simulations2026
- Predicting Atomic Charges in MOFs by Topological Charge EquilibrationAdvanced Functional Materials, 2026
- An Introduction to the Combustion of Carbon MaterialsChemistry – A European Journal, 2022
- Extended tight‐binding quantum chemistry methodsWIREs Computational Molecular Science, 2020
- A complete review based on various aspects of pulverized coal combustionInternational Journal of Energy Research, 2019
- New approaches to measuring biochar density and porosityBiomass and Bioenergy, 2014
- Molecular simulation of the CH4/CO2/H2O adsorption onto the molecular structure of coalScience China Earth Sciences, 2014
- PACKMOL: A package for building initial configurations for molecular dynamics simulationsJournal of Computational Chemistry, 2009
- VMD: Visual molecular dynamicsJournal of Molecular Graphics, 1996
- Heterogeneous kinetics of coal char gasification and combustionProgress in Energy and Combustion Science, 1978