  In this paper, we establish a H\"older-type quantitative estimate of unique
continuation for solutions to the heat equation with Coulomb potentials in
either a bounded convex domain or a $C^2$-smooth bounded domain. The approach
is based on the frequency function method, as well as some parabolic-type Hardy
inequalities.
