In this study we investigated the lowest energy structures of neutral and cationic Li6Hn (n=0–7) clusters using ab initio calculations based on DFT and CCSD(T) methods. We used the stochastic search method to obtain the global minimum on the potential energy surface of hydrogenated lithium clusters. For the lowest-energy isomers of the clusters, the relative stabilities, electronic and reactivity properties are presented depending on the binding energies, second energy difference, chemical hardness, chemical potential, softness, electrophilicity index, ionization potentials, electron affinities and HOMO–LUMO energy gaps. We also studied the fragmentation energies of neutral and cationic hydrogenated lithium clusters. In particular, the fragmentation energies of cationic clusters are compared with experimental data. The results are in good agreement with most available experimental findings.