We give two-sided, global (in all variables) estimates of the heat kernel and the Green function of the fractional Schrödinger operator with a non-negative and locally bounded potential $V$ such that $V(x) \to \infty$ as $|x| \to \infty$. We assume that $V$ is comparable to a radial profile with the doubling property. Our bounds are sharp with respect to spatial variables and qualitatively sharp with respect to time. The methods we use combine probabilistic and analytic arguments. They are based on the strong Markov property and the Feynman--Kac formula.