Sharp estimates for the gradient of solutions to the heat equation

Various sharp pointwise estimates for the gradient of solutions to the heat equation are obtained. The Dirichlet and Neumann conditions are prescribed on the boundary of a half-space. All data belong to the Lebesgue space $L^p$. Derivation of the coefficients is based on solving certain optimization problems with respect to a vector parameter inside of an integral over the unit sphere.