  We prove that the integral operators associated with the layer heat
potentials depend smoothly upon a parametrization of the support of
integration. The analysis is carried out in the optimal H\"older setting.
