The electronic work function and surface energy are calculated for the
principal surfaces of body-centered-cubic Li. Results for three-layer
Li slabs with the bulk lattice parameter show that the work function
increases in the order (100), (111), and (110) in the surface orientat
ion, which differs from the prediction of the jellium model with pseud
opotential corrections. The results of the structural optimization sho
w that the dependence of the work function and surface energy on the s
urface orientation almost vanishes for very thin slabs. Energetically
the most stable one of these three investigated slabs is the (111) sla
b. In all three cases, the central layer seems to have a tendency towa
rd bulklike behavior, which Significantly affects the relaxation of th
e slab.