Let m ∈ N, P(D) := |α|=2m (−1) m a α D α be a 2m-order homogeneous elliptic operator with real constant coefficients on R n , and V a measurable function on R n . In this article, the authors introduce a new generalized Schechter class concerning V and show that the higher order Schrödinger operator L := P(D) + V possesses a heat kernel that satisfies the Gaussian upper bound and the Hölder regularity when V belongs to this new class. The Davies-Gaffney estimates for the associated semigroup and their local versions are also given. These results pave the way for many further studies on the analysis of L.