The Binding Energies of Atoms on Amorphous Silicate Dust: A Computational Study

TL;DR

This study computationally determines the binding energies of key interstellar atoms on amorphous silicate dust, providing crucial data for understanding dust stability, growth, and chemical reactions in space.

Contribution

First-principles calculations of atom binding energies on amorphous silicate dust, filling a gap in data for astrophysical dust models.

Findings

Silicon, aluminum, and calcium have high binding energies, enabling stability at high temperatures.

Most other atoms exhibit weaker binding energies but remain stable against sublimation in the ISM.

Binding energies follow a log-normal distribution for most elements, informing dust evolution models.

Abstract

Context. We investigate the binding energies of atoms to interstellar dust particles, which play a key role in their growth and evolution, as well as for the chemical reactions on their surfaces. Aims. We aim to compute the binding energies of abundant atoms in the interstellar medium (C, N, O, Mg, Al, Si, S, Ca, Fe, and Ni) to silicate dust. Methods. We used the Geometries, Frequencies, and Non-covalent Interactions Tight Binding (GFN1-xTB) method to compute the binding energies. An FeMgSiO$_4$ periodic surface model, containing 81 local minima on the surface, was used. Results. A range of binding energies was found for each element. The median of the binding energies follows the order Si (14.8 eV) > Al (12.8 eV) > Ca (12.7 eV) > C (9.5 eV) > O (8.1 eV) > N (6.2 eV) > Fe (6.0 eV) > S (5.2 eV) > Mg (2.4 eV). The probability distribution of binding energies for each element except Ca is statistically consistent with a log-normal distribution. Conclusions. In general, Si, Ca, and Al atoms have large binding energies. Thus, these atoms can stay on the silicate dust particles at high temperatures. The binding energies of the other atoms, C, N, O, Mg, S, Fe and Ni, are relatively weak. However, the computed binding energies for these elements are still far stronger than the energies associated with dust temperatures typical of the ambient interstellar medium (ISM), suggesting that silicate grains are generally stable against sublimation. We estimate sublimation temperatures for silicate grains to range from 1600 K to 3000K depending on assumed grain size and lifetime. These binding energies on silicate dust grains, estimated from first principles for the first time, provide invaluable input to models of dust evolution and dust-catalyzed chemical reactions in the interstellar medium.