Effect of Ice Charging on the Astrochemistry of Interstellar Sulfur Bearing Species on Amorphous Solid Water

We aim to derive statistically robust and physically interpretable BE distributions for atomic S and the sulfur-bearing molecules H$_2$S, SO$_2$, and OCS on neutral and negatively charged ASW, and assess how excess negative charge alters their retention and potential role in the sulfur reservoir of cold dense molecular clouds. Basis set superposition error (BSSE) and zero-point energy (ZPE) corrected BEs of S, SO$_2$, OCS and H$_2$S, are calculated by density functional theory (DFT), using the ORCA software. Molecule-specific DFT levels of theory were first selected from benchmark calculations on neutral and charged small water complexes, against coupled-cluster reference energies. The selected protocols were then applied to study adsorption on neutral and charged ASW clusters. A range of adsorption sites were sampled on five independent amorphous ice clusters, yielding BE distributions that account for the site heterogeneity of ASW. Neutral ASW yields broad, site-dependent BE distributions consistent with previous water-ice estimates. On charged ASW, three general cases are identified: the BE will either always increase due to electron transfer (S-atom, SO$_2$), the BE increases slightly without any electron transfer (H$_2$S) or the BE remains the same unless an electron transfer occurs under specific conditions (OCS). These findings are inherently linked to the molecular properties of these molecules.