The structural and energetic proprieties for the Li + Xen (n = 1–18) clusters are investigated using both Basin- Hopping combined with Potential Model description (BH-PM) and DFT methods. A structural transition from tetrahedral (4 coordination) form to octahedral (6 coordination) one is observed for n = 6. Above this size, all structures have an octahedral core. The cubic-face-centered arrangement for xenon atoms is detected for Li + Xe14. To the best of our knowledge, the Li + Xen (n = 1–18) clusters are studied in the present work for the first time using the DFT theoretical approach. The M062X functional combined with aug-cc-pVDZ (for Li) and def2- TZVP (for Xe) basis sets reproduces accurately the CCSD(T) potential energy curve of Li + Xe system. Atom- Centered Density Matrix Propagation (ADMP) molecular dynamic calculations have been carried. Moreover, we investigate the larger sizes n = 31–35, 44, and 55 for the first time using the BH-PM theoretical approach. The closing of the first and second octahedron shells are proved for the n = 6 and 34 sizes, respectively. The relative stabilities of the Li + Xen molecules are also studied by computing the total energy, the binding energy per atoms for each size n. Then, the second energy difference between the size n and its two near neighbors allows iden- tifying the magic number series. Our present data are analyzed, discussed and compared with the available theoretical and experimental data.