ABSTRACT:The dependence of the harmonic oscillator (HO) energy level spacinghω on the particle number N is studied analytically for atomic (metal) clusters on the basis of their electronic densities, parametrizing Ekardt's results (for sodium clusters) by means of a Fermi distribution. An interesting feature of such an approach is that it leads, under the assumptions made, to "kinks," that is, to "marked discontinuities in the slope" ofhω at the closed shells. These discontinuities diminish as N increases. Key words: atomic (metal) clusters; sodium clusters; jellium model; harmonic oscillator; Fermi distribution T he dependence of the harmonic oscillator (HO) energy-level spacinghω on the particle number N (that is, on the number of the valence electrons of the atoms in the jellium model) for atomic (metal) clusters has attracted theoretical interest [1 -5]. Related experience from nuclear [6,7] or hypernuclear [8,9] physics has been very valuable, and the methods used in these fields have been adjusted to the cluster case. Two general approaches have been followed.In the first one [1 -3], which has been mainly discussed in the case of large values of N, the energy scale is related to the size of the cluster, as this is determined by the mean-square radius of the Correspondence to: M. E. Grypeos; e-mail: grypeos@auth.gr. electronic density distribution. It has been assumed that this density has a constant value equal to the bulk conduction electron density ρ B , and the size of the (spherical) cluster is normalized so that it contains N valence electrons, that is, R 0 = r s N 1/3 where r s = (3/4πρ B ) 1/3 is the Seitz-Wigner radius. The expression for large N which results in this way is where r s is in atomic units andhω in eV. An improved expression has also been derived by taking into account the so-called spill-out [3,10,11] of the electrons beyond the boundary R 0 . This spill-out