Euclidean scalar-gradient stress identity (${\partial }_k{\Sigma }_{kj}=(\Delta \psi ){\partial }_j\psi$)

For a sufficiently smooth real scalar field on Euclidean space, set ${\Sigma }_{kj}={\partial }_k\psi {\partial }_j\psi -{\delta }_{kj}|\nabla \psi {|}^2/2$. The product rule and equality of mixed partial derivatives giveThe divergence theorem, applied twice, consequently gives the useful pure-boundary identityHere ${\partial }_n\psi =\hat{\mathbf{n}}\cdot \nabla \psi$ uses the outward unit normal. This is the Euclidean zero-potential counterpart of the scalar stress-energy divergence identity.