Abstract: | Two electron systems confined by harmonic potential is known as harmonium. Such a system has been studied for many reasons in the literature. In this work we study harmonium under Debye potential. We use higher order finite difference method for the solution of Schrodinger equation. Complete energy spectrum of harmonium and harmonium under Debye potential is studied. Debye screening length shows considerable effect on the energy levels and the radial matrix elements. The results are analysed in the light of existing results and the comparison with available results shows remarkable agreement. |