Density functional calculations using B3LYP/6-311G method have been carried out for small to medium-sized lithium clusters (Li N , N ¼ 2-30). The optimized geometries of neutral and singly charged clusters, their binding energies, ionization potential, electron affinity, chemical potential, softness, hardness, highest occupied molecular orbital and lowest unoccupied molecular orbital (HOMO-LUMO) gap, and static dipole polarizability have been investigated systematically. In addition, we study the distribution of partial charges in detail using natural population analysis (NPA) in small-sized clusters (Li N , N ¼ 2-10), both neutral and cationic, and demonstrate the correlation between symmetry and charge. Uniform distribution of charges in cationic clusters confirms them to be energetically more favorable than the neutral counterparts. Whenever possible, results have been compared with available data. An excellent agreement in every case supports new results as reliable predictions. A careful study of optimized geometries shows that Li 9 is derivable from bulk Li structure, i.e., body centered cubic cell, and higher clusters have optimized shapes derived from this. Further, the turnover form two to three dimensional structure occurs at cluster size N ¼ 6. The quantity a 1/3 (a ¼ polarizability) per atom is found to be broadly proportional to softness (per atom) as well as inverse ionization potential (per atom). The present work forms a sound basis for further study of large-sized clusters as well as other atomic clusters.