I'm working on Problem 7-7 in Lee's "Introduction to Riemannian Manifolds", which asks us to prove Bochner's formula: for a Riemannian manifold $(M,g)$ and $u \in C^\infty(M)$, $$ \Delta \left(\frac 1 2 |\mathrm{grad}\: u|^2\right) = \left|\nabla^2 u\right|^2 + \left\langle \mathrm{grad}\:(\Delta u), \mathrm{grad}\:u\right\rangle + Rc(\mathrm{grad}\:u, \mathrm{grad}\:u) $$ where $\Delta u = \mathrm{div}\:\mathrm{grad}\:u$ is the Laplacian of $u \in C^\infty(M)$, $\nabla^2 u = u_{;ij} dx^i \otimes dx^j$ is the covariant Hessian (where $u_{;ij} = \partial_j\partial_i u - \Gamma_{ji}^k \partial_k u$), and $Rc = R_{ij} dx^i \otimes dx^j$ is the Ricci curvature, where $$ R_{ij} = R_{kij}^{\:\:\:\:k} $$ and $R_{ijk}^{\:\:\:\:l}$ are the coefficients of the curvature endomorphism $$ R(X,Y)Z = \nabla_X \nabla_Y Z - \nabla_Y \nabla_X Z - \nabla_{[X,Y]}Z. $$ Lee suggests using the following two facts:
- $\Delta u = g^{ij} u_{;ij} = u_{;i}^{\:\,i}$
- If $\beta$ is a smooth 1-form on $M$, then $$\nabla^2_{X,Y}\beta - \nabla^2_{Y,X} \beta = -R(X,Y)^*\beta,$$ or in coordinates, $$ \beta_{j;pq} - \beta_{j;qp} = R_{pqj}^{\:\:\:\,m}\beta_m $$ where $\beta_{j;pq}$ are the coefficients of $\nabla^2\beta$.
I've tried deriving Bochner's formula from a variety of calculations, mostly involving Riemannian normal coordinates $(x^i)$ at a point $x \in M$. I've used the first fact to expand both sides but the right side especially gets pretty hairy even with normal coordinates. I am really not sure where the second fact comes into play. Any suggestions?