Hardy perturbations of subordinated Bessel heat kernels

Motivated by the spectral theory of relativistic atoms, we prove matching upper and lower bounds for the transition density of Hardy perturbations of subordinated Bessel heat kernels. The analysis is based on suitable supermedian functions, in particular invariant functions.