Let $L=-\mathsf{\Delta }+V$ be a Schrödinger operator, where the potential V belongs to the reverse Hölder class. By the subordinative formula, we introduce the fractional heat semigroup ${\left\{e^{-t{L}^{\alpha }}\right\}}_{t>0}, 0<\alpha <1$, associated with L. By the aid of the fundamental solution of the heat equation: ${\partial }_{t}u+Lu={\partial }_{t}u-\mathsf{\Delta }u+Vu=0,$ we estimate the gradient and the time-fractional derivatives of the fractional heat kernel $K_{\alpha ,t}^L(·,·)$, respectively. This method is independent of the Fourier transform, and can be applied to the second-order differential operators whose heat kernels satisfy the Gaussian upper bounds. As an application, we establish a Carleson measure characterization of the Campanato-type space $BM{O}_L^{\gamma }\left({\mathbb{R}}^n\right)$ via the fractional heat semigroup ${\left\{e^{-t{L}^{\alpha }}\right\}}_{t>0}$.