  This paper is devoted to the study of gradient estimates for the Dirichlet
problem of the heat equation in the exterior domain of a compact set. Our
results describe the time decay rates of the derivatives of solutions to the
Dirichlet problem. Applications of these estimates to bilinear type commutator
estimates for Laplace operator with Dirichlet boundary condition in exterior
domain are discussed too.
