In topological n-type thermoelectrics based on Bi2Te3 with atomic substitutions Bi → In, Te → Se, S, the morphology and the surface states of Dirac fermions on the interlayer (0001) surface of van der Waals were studied by scanning tunneling microscopy and spectroscopy (STM/STS) techniques. By the STM method, the dark and light spots on the surface were found, which intensities depend on the composition and thermoelectric properties of solid solutions such as the Seebeck coefficient and thermoelectric power factor. The observed surface morphology features in the solid solutions are explained by distortions of surface electronic states originated from atomic substitutions, the influence of doping impurity, and formation mainly of substitutional impurity defects in thermoelectrics. The dips associated with substitutional impurities and antisite defects were found from the analysis of the height profiles obtained on the (0001) surface. Fast Fourier transform of the morphology STM images of the (0001) surface were used to obtain the interference patterns of the quasiparticles excitation caused by surface electrons scattering by defects. The Dirac point energy and its fluctuations, peak energies of surface defects, the positions of the valence and conduction band edges, and the energy gap were determined from an analysis of tunneling spectra. A correlation between the parameters of surface states of Dirac fermions and thermoelectric properties was found. Thus, a contribution of the fermions surface states increases with rise of the surface concentration in solid solutions with high power factor, and the largest concentration value was observed in the Bi1.98In0.02Te2.85Se0.15 composition. The dependences of Fermi energy on the wave vector for different solid solutions are described by a set of Dirac cone sections located within the limits of the fluctuations of the Dirac point energy that explained by weak changes of the Fermi velocity for studied atomic substitutions in sublattices of bismuth telluride.