Nonrelativistic energies for the low-lying states of lithium are calculated using the variational method in Hylleraas coordinates. Variational eigenvalues for the infinite nuclear mass case with up to 34 020 terms are −7.478 060 323 910 147(1) a.u. for 1s 2 2s 2 S, −7.354 098 421 444 37(1) a.u. for 1s 2 3s 2 S, −7.318 530 845 998 91(1) a.u. for 1s 2 4s 2 S, −7.410 156 532 652 41(4) a.u. for 1s 2 2p 2 P , and −7.335 523 543 524 688(3) a.u. for 1s 2 3d 2 D. The selection of the minimum set of angular momentum configurations is discussed, with the 2P and 3D states as examples to demonstrate the impact of various configurations on the variational energies. It is shown by numerical example that the second spin function (i.e., coupled to form a triplet intermediate state) has no significant effect on either the variational energies or the spin-dependent Fermi contact term. Results of greatly improved accuracy for the Fermi contact term are presented for all the states considered.