Bimetallic clusters aroused tremendous interest because of the property changes that occur: structure, size, and composition. Herein, a structural search of the global minimum for anionic LiMg
n
– (n=2–11) clusters was performed using an efficient CALYPSO structural searching program with subsequent DFT calculations. A great variety of low energetic isomers converged, and the most stable ones were confirmed by comparing their total energy for each size. It is found that the iMg
n
– clusters are structurally consistent with corresponding Mg clusters anions except for LiMg5
– and LiMg7
–. In all the doped clusters, the Li atom prefers to occupy the convex position. Simulated PES, IR, and Raman spectra of iMg
n
– could be used as an essential evidence for identifying cluster structures experimentally in the future. Stability study revealed that a tower-like structure of LiMg9
– has prominent stability and can be identified as a magic number cluster. The reason might be that there are both closed-shell 1S21P61D102S2 electronic configurations and stronger Li-Mg bonds caused by sp hybridization in the LiMg9
– cluster.