AbstractGradient enhanced theories of crystal plasticity enjoy great research interest. The focus of this work is on thermodynamically consistent modeling of grain size dependent hardening effects. In this contribution, we develop a model framework for damage coupled to gradient enhanced crystal thermoplasticity. The damage initiation is directly linked to the accumulated plastic slip. The theoretical setting is that of finite strains. Numerical results on single-crystalline metal showing the development of damage conclude the paper.