In this paper, an h∕p spectral element method with least-square formulation for parabolic interface problem will be presented. The regularity result of the parabolic interface problem is proven for non-homogeneous interface data. The differentiability estimates and the main stability estimate theorem, using non-conforming spectral element functions, are proven. Error estimates are derived for h and p versions of the proposed method. Specific numerical examples are given to validate the theory. © 2018 Elsevier B.V.