Density functional theory was used to examine the solvation number and aggregation state of several alkyllithium compounds in clusters with tetrahydrofuran molecules coordinated to each lithium atom. We then made the microsolvation approximation and approximated the bulk free energy of solvation by the free energy of clustering with solvent molecules in the gas phase. The trends in the computed results are in reasonable agreement with the available experimental data.