We investigate a control problem for the heat equation. The goal is to find an optimal heat transfer coefficient in the dynamic boundary condition such that a desired temperature distribution at the boundary is adhered. To this end we consider a function space setting in which the heat flux across the boundary is forced to be an L (p) function with respect to the surface measure, which in turn implies higher regularity for the time derivative of temperature. We show that the corresponding elliptic operator generates a strongly continuous semigroup of contractions and apply the concept of maximal parabolic regularity. This allows to show the existence of an optimal control and the derivation of necessary and sufficient optimality conditions.