Here: $$ \
abla \cdot (D\
abla\psi)= (-{\rm i})\
abla \cdot (D\psi \vec{k})=(-{\rm i})\sum_i \frac{\partial}{\partial x_i}(D\psi \vec{k})_i=(-{\rm i})\sum_{i,j} \frac{\partial}{\partial x_i}D_{ij}\psi k_j$$
$$=-\sum_{i,j} D_{ij}\psi k_ik_j=-(k^TDk)\psi$$