By using density functional theory with Becke-Lee-Yang-Parr (BLYP) gradient
correction, we have studied the structural and electronic properties of be
ryllium clusters up to 21 atoms. The theoretical calculations successfully
reproduce the experimental bond length and vibrational frequency of the Bee
dimer and the overestimation of the binding energy is significantly reduce
d compared to that in previous theoretical work. The Be clusters with n = 4
, 10, 17 show particularly high stability, consistently with the electron s
hell model. The size evolution of the electronic properties from van der Wa
als to covalent and bulk metallic behaviour in Be clusters is discussed.