In the relativistic mean field (RMF) calculations usually the basis expansi
on method is employed. For this one uses single harmonic oscillator (HO) ba
sis functions. A proper description of the ground state nuclear properties
of spherical nuclei requires a large (around 20) number of major oscillator
shells in the expansion. In halo nuclei where the nucleons have extended s
patial distributions, the use of single HO basis for the expansion is inade
quate for the correct description of the nuclear properties, especially tha
t of the surface region. In order to rectify these inadequacies, in the pre
sent work an orthonormal basis composed of two HO basis functions having di
fferent sizes is proposed. It has been shown that for a typical case of (A
= 11) the ground state constructed using two-HO wave functions extends much
beyond the second state or even third excited state of the single HO wave
function. To demonstrate its usefulness explicit numerical RMF calculations
have been carried out using this procedure for a set of representative sph
erical nuclei ranging from O-16 to Pb-208. The binding energies, charge rad
ii and density distributions have been correctly reproduced in the present
scheme using a much smaller number of major shells (around 10) in the expan
sion.