Solvation of lithium ion in helium clusters: Structural properties and relative stabilities

Structural study and relative stabilities of Li⁺-doped helium clusters Li⁺He_n (n = 1-18) has been reported in this work using two theoretical protocols. The first one is based on the basin-hopping optimization technique, where the total energy of each cluster is described by an additive model describing Li⁺-He and He-He interactions. The second one is the DFT calculations, in which the initial structures are generated by ABCluster algorithm and CALYPSO software. The CSA shape was found where the first solvation shell is completed at n = 10. The relative stabilities of Li⁺He_n (n = 1-18) clusters have been discussed based on the variation of the binding energy, second-order difference in energy, fragmentation energy and HOMO-LUMO energy gap as a function of the cluster size. The results showed that Li⁺He₁₀ is the most stable cluster. The dipole moment is calculated and showed the polar character of the Li⁺He_n clusters. Finally, the interatomic interactions have been examined topologically by the means of AIM and non-covalent reduced density gradient (NC-RDG) analyses.