Exploring the stability and aromaticity of rare earth doped tin cluster MSn16− (M = Sc, Y, La)

molecular dynamics simulations, the spherical jellium model, adaptive natural density partitioning, localized orbital locator, and electron density difference are employed to analyze the stability of these clusters. The aromaticity of the clusters is examined using iso-chemical shielding surfaces and the gauge-including magnetically induced currents. This study demonstrates that the stability and aromaticity of a tin cage can be systematically adjusted through doping.