The efficient relativistic Dirac-Hartree-Fock (DHF) and Dirac-Kohn-Sham (DKS) methods are proposed by an application of the pseudospectral (PS) approach. The present PS-DHF/DKS method is a relativistic extension of the PS-HF/KS method of Friesner, though we aim at higher numerical accuracy by elimination of superfluous arbitrariness. The relativistic PS-DHF/DKS method is implemented into our REL4D programs. Several PS applications to molecular systems show that the relativistic PS-DHF/DKS approach is more efficient than the traditional approach without a loss of accuracy. The present PS-DKS method successfully assigns and predicts the photoelectron spectra of hexacarbonyl complexes of tungsten and seaborgium theoretically. (C) 2004 American Institute of Physics.