A Fully ab Initio Kinetic Monte Carlo Approach for Modeling Adsorption and Diffusion in Interstellar Icy Grain Mantles: The Case of H <sub>2</sub> S

Understanding diffusion on interstellar ices is key to modeling the chemical evolution of cold molecular clouds, where low temperatures severely limit molecular mobility. In this study, we introduce a robust and fully automated multiscale computational framework to quantify diffusion processes of adsorbates at the surface of amorphous solid water (ASW). Using H<sub>2</sub>S as a test case, whose binding sites were previously studied at the ab initio level, we constructed a detailed network of 141 adsorption sites connected by over 270 transition states. All density functional energetics were benchmarked against DLPNO-CCSD-(T), achieving chemical accuracy in the description of diffusion barriers, which span from 0.1 to 27 kJ mol<sup>-1</sup> with a median value of 5.4 kJ mol<sup>-1</sup>. An off-lattice kinetic Monte Carlo (kMC) model adopting both the ab initio diffusion barriers and binding energies for the desorption processes was carried out to compute temperature-dependent diffusion coefficients and to reconstruct the temperature-programmed desorption (TPD) curve. Our simulations reveal that thermal diffusion of H<sub>2</sub>S is negligible below 20 K, with diffusion coefficients as low as 10<sup>-48</sup> cm<sup>2</sup> s<sup>-1</sup> at 10 K, thus excluding Langmuir-Hinshelwood surface encounters under typical dense cloud conditions. Moreover, under submonolayer conditions, diffusion was found to have negligible influence on the reconstructed TPD peak position. Furthermore, our results demonstrate that a universal scaling factor <i>f</i> to guess the diffusion barriers (Δ<i>E</i> <sub>diff</sub>) from the sole knowledge of BE: <i>f</i> = Δ<i>E</i> <sub>diff</sub>/BE does not apply as it exhibits wide variability across the sampled configurations. These findings highlight the need for incorporating statistically meaningful distributions of binding energies and diffusion barriers in astrochemical models to more accurately capture diffusion and surface reactivity on interstellar ices.

Authors

Institutions

Publication Details

Journal
ACS Earth and Space Chemistry
Published
2025-12-18
DOI
https://doi.org/10.1021/acsearthspacechem.5c00208
Primary Topic
Astrophysics and Star Formation Studies
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

A Fully ab Initio Kinetic Monte Carlo Approach for Modeling Adsorption and Diffusion in Interstellar Icy Grain Mantles: The Case of H 2 S

Piero Ugliengo, Stefano Pantaleone, Albert Rimola, C. Ceccarelli et al.
ACS Earth and Space Chemistry
Astrophysics and Star Formation Studies
article

A Fully ab Initio Kinetic Monte Carlo Approach for Modeling Adsorption and Diffusion in Interstellar Icy Grain Mantles: The Case of H 2 S

Piero Ugliengo, Stefano Pantaleone, Albert Rimola, C. Ceccarelli, Vittorio Bariosco
article en

Abstract

Understanding diffusion on interstellar ices is key to modeling the chemical evolution of cold molecular clouds, where low temperatures severely limit molecular mobility. In this study, we introduce a robust and fully automated multiscale computational framework to quantify diffusion processes of adsorbates at the surface of amorphous solid water (ASW). Using H<sub>2</sub>S as a test case, whose binding sites were previously studied at the ab initio level, we constructed a detailed network of 141 adsorption sites connected by over 270 transition states. All density functional energetics were benchmarked against DLPNO-CCSD-(T), achieving chemical accuracy in the description of diffusion barriers, which span from 0.1 to 27 kJ mol<sup>-1</sup> with a median value of 5.4 kJ mol<sup>-1</sup>. An off-lattice kinetic Monte Carlo (kMC) model adopting both the ab initio diffusion barriers and binding energies for the desorption processes was carried out to compute temperature-dependent diffusion coefficients and to reconstruct the temperature-programmed desorption (TPD) curve. Our simulations reveal that thermal diffusion of H<sub>2</sub>S is negligible below 20 K, with diffusion coefficients as low as 10<sup>-48</sup> cm<sup>2</sup> s<sup>-1</sup> at 10 K, thus excluding Langmuir-Hinshelwood surface encounters under typical dense cloud conditions. Moreover, under submonolayer conditions, diffusion was found to have negligible influence on the reconstructed TPD peak position. Furthermore, our results demonstrate that a universal scaling factor <i>f</i> to guess the diffusion barriers (Δ<i>E</i> <sub>diff</sub>) from the sole knowledge of BE: <i>f</i> = Δ<i>E</i> <sub>diff</sub>/BE does not apply as it exhibits wide variability across the sampled configurations. These findings highlight the need for incorporating statistically meaningful distributions of binding energies and diffusion barriers in astrochemical models to more accurately capture diffusion and surface reactivity on interstellar ices.

ACS Earth and Space ChemistryVol. 10(1)
Universitat Autònoma de Barcelona (ES), Institut de Planétologie et d'Astrophysique de Grenoble (FR), Ospedale Maria Vittoria (IT), Institut des langues et cultures d'Europe, Amérique, Afrique, Asie et Australie (FR), University of Turin (IT), Université Grenoble Alpes (FR)
Openalex Percentile: Top 72%
Astrophysics and Star Formation Studies
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.