Search

Sputtering by solar wind ions is a key process driving the ejection of high-energy particles into the exospheres of airless bodies like asteroids, Mercury and the Moon. In view of upcoming missions which will deliver new in-situ data on these exospheres like the Artemis program at the Moon and BepiColombo at Mercury, a deeper understanding of sputtering effects is crucial. In this work, we combine sensitive quartz crystal microbalance measurements and numerical simulations to quantify sputter yields of Apollo soil 68501 under solar wind relevant conditions. We find that none of the commonly used simulation codes can reliably predict laboratory sputter yields without experimental benchmarks. All of the employed packages significantly overestimate the sputter yields of flat samples by at least a factor of 2 for the case of hydrogen. When accounting for surface roughness and regolith-like porosity, sputter yields are decreased even further to 7.3E-3 atoms\ion and 7.6E-2 atoms\ion for H and He at solar wind energies of 1 keV\amu, respectively. The reduced yields of porous regolith structures are largely independent of the ion incidence angle, making them applicable across a wide range of lunar latitudes. This study highlights the need for experimental validation of sputtering models to ensure accurate predictions for space weathering and lunar exosphere composition.

We perform Low Energy Ion Scattering with 1\,keV He ions on CaSiO$_3$ using a commercial electrostatic detector system and determine the charge fraction of scattered ions from comparison with Binary Collision Approximation simulations. The simulations take dynamical surface changes due to surface cleaning Ar sputtering into account and scattered He particles are separated into single, dual, and multiple scattering trajectories. We find that the charge fraction of single and dual scattered He is about 10 times higher than the one for multiple collisions. Our results show that quantitative concentration profiles can be inferred from this method, if the charge fraction components are determined first.

Rocky planets and moons experiencing solar wind sputtering are continuously supplying their enveloping exosphere with ejected neutral atoms. To understand the quantity and properties of the ejecta, well established Binary Collision Approximation Monte Carlo codes like TRIM with default settings are used predominantly. Improved models such as SDTrimSP have come forward and together with new experimental data the underlying assumptions have been challenged. We introduce a hybrid model, combining the previous surface binding approach with a new bulk binding model akin to Hofs\"ass & Stegmaier (2023). In addition, we expand the model implementation by distinguishing between free and bound components sourced from mineral compounds such as oxides or sulfides. The use of oxides and sulfides also enables the correct setting of the mass densities of minerals, which was previously limited to the manual setting of individual atomic densities of elements. All of the energies and densities used are thereby based on tabulated data, so that only minimal user input and no fitting of parameters are required. We found unprecedented agreement between the newly implemented hybrid model and previously published sputter yields for incidence angles up to 45{\deg} from surface normal. Good agreement is found for the angular distribution of mass sputtered from enstatite MgSiO$_3$ compared to latest experimental data. Energy distributions recreate trends of experimental data of oxidized metals. Similar trends are to be expected from future mineral experimental data. The model thus serves its purpose of widespread applicability and ease of use for modelers of rocky body exospheres., Comment: 23 pages, 6 figures, 3 tables