In this paper, we will give a horizontal gradient estimate of positive solutions of $Delta_b u = - lambda u$ on complete noncompact pseudo-Hermitian manifolds. As a consequence, we recapture the Liouville theorem of positive pseudo-harmonic functions on Sasakian manifolds with nonnegative pseudo-Hermitian Ricci curvature.