Catalytic formation of H2 on carbonaceous dust grains - implications for interstellar observations
4 affiliations
- 1Park Tudor School , 7200 N College Ave, Indianapolis, IN 46240 ,
- 2Fairview High School , 1515 Greenbriar Blvd, Boulder, CO 80305 ,
- 3Texas Southern University Department of Physics, , 3100 Cleburne Street, Houston, TX 77004 ,
- 4ITAMP, Center for Astrophysics | Harvard & Smithsonian , 60 Garden Street, Cambridge, MA 02138 ,
Abstract
We use kinetic Monte Carlo (KMC) simulations to study molecular hydrogen formation on carbonaceous dust grain surfaces, validated against recent laboratory measurements of H2 formation on coronene films at temperatures from 10 to 250 K. The model uses a three-dimensional amorphous carbon lattice with heterogeneous physisorption (45 ± 5 meV) and chemisorption (1.75 ± 0.25 eV) sites, and tracks both Langmuir–Hinshelwood (LH) and Eley–Rideal (ER) formation channels within a stochastic Gillespie event-driven framework. The model reproduces the measured efficiency curve within the experimental uncertainties, including the isothermal (constant surface temperature) measurements at 100–250 K. The simulations correctly describe the phase boundary between the LH and ER driven processes as functions of grain temperature and the observed crossover. Under interstellar medium conditions, 10–250 K and nH = 10–104 cm−3, the model predicts three distinct regimes for the formation efficiency ε, the fraction of impinging H atoms released as H2. At 10 K diffusion is slow and ε ≈ 0.06. Between 20 K and 80 K, LH dominates and ε ≈ 0.28. Above 150 K, an ER plateau at ε = 0.19 is sustained by chemisorption-trapped H atoms. The LH-to-ER crossover occurs between 100 and 120 K. At 100 K we observe a 16 per cent density-dependent stochastic enhancement, which rate-equation models cannot capture. At Tdust = 60 K, nH = 103 cm−3 we find the ratio of H2 formation to free-fall time tH2/tff ≈ 0.93, so dust-catalysed H2 chemistry can keep pace with gravitational collapse in high-redshift star-forming environments.
Related
- Carbonaceous Chondrites Provide Evidence for Late-stage Planetesimal Formation in a Pressure Bump
- Visible, near‐, and thermal infrared spectra of asteroid Bennu samples: Relationship to and implications for remote sensing of carbonaceous asteroids
- The bar-driven central complex of NGC 6946: A panchromatic and spatially resolved view of its star formation history and interstellar medium
- JWST/MIRI Observations of Amorphous Silicates in Interstellar Comet 3I/ATLAS
- A Search for Persistent Rotational Modulation in the Dedicated TESS Observations of Interstellar Comet 3I/ATLAS
BibTeX
@article{das2026catalytic,
title = {{Catalytic formation of H2 on carbonaceous dust grains - implications for interstellar observations}},
author = {Aryav Das and Julia Gao and Daniel Vrinceanu and H R Sadeghpour},
journal = {Monthly Notices of the Royal Astronomical Society},
year = {2026},
month = oct,
doi = {10.1093/mnras/stag2052},
url = {https://academic.oup.com/mnras/advance-article/doi/10.1093/mnras/stag2052/9004096},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Catalytic formation of H2 on carbonaceous dust grains - implications for interstellar observations
AU - Aryav Das
AU - Julia Gao
AU - Daniel Vrinceanu
AU - H R Sadeghpour
PY - 2026
DA - 2026/10/28
JO - Monthly Notices of the Royal Astronomical Society
DO - 10.1093/mnras/stag2052
UR - https://academic.oup.com/mnras/advance-article/doi/10.1093/mnras/stag2052/9004096
AB - We use kinetic Monte Carlo (KMC) simulations to study molecular hydrogen formation on carbonaceous dust grain surfaces, validated against recent laboratory measurements of H2 formation on coronene films at temperatures from 10 to 250 K. The model uses a three-dimensional amorphous carbon lattice with heterogeneous physisorption (45 ± 5 meV) and chemisorption (1.75 ± 0.25 eV) sites, and tracks both Langmuir–Hinshelwood (LH) and Eley–Rideal (ER) formation channels within a stochastic Gillespie event-driven framework. The model reproduces the measured efficiency curve within the experimental uncertainties, including the isothermal (constant surface temperature) measurements at 100–250 K. The simulations correctly describe the phase boundary between the LH and ER driven processes as functions of grain temperature and the observed crossover. Under interstellar medium conditions, 10–250 K and nH = 10–104 cm−3, the model predicts three distinct regimes for the formation efficiency ε, the fraction of impinging H atoms released as H2. At 10 K diffusion is slow and ε ≈ 0.06. Between 20 K and 80 K, LH dominates and ε ≈ 0.28. Above 150 K, an ER plateau at ε = 0.19 is sustained by chemisorption-trapped H atoms. The LH-to-ER crossover occurs between 100 and 120 K. At 100 K we observe a 16 per cent density-dependent stochastic enhancement, which rate-equation models cannot capture. At Tdust = 60 K, nH = 103 cm−3 we find the ratio of H2 formation to free-fall time t_H_2/t_ff ≈ 0.93, so dust-catalysed H2 chemistry can keep pace with gravitational collapse in high-redshift star-forming environments.
SN - 0035-8711
SN - 1365-2966
ER -