There are many questions on the proof of the mentioned formula already on the site, for example: Question 1, Question 2. Somehow I fail to prove it while trying to perform the computation in normal coordinates. To be precise, we want to show
$ \Delta(f_i) = (\Delta f)_i + r_{ij}f_j$
where $r_{ij}$ is Ricci curvature components in normal coordinates and derivatives are coordinate derivatives. One easily sees in normal coordinates:
$\Delta(f_i)=f_{ikk}$
while
$ ( \Delta f )_i = f_{kki} - \Gamma_{kk,i}^m f_m$
so it would suffice to show:
$r_{im} - \Gamma_{kk,i}^m=0.$
Where the problem in my calculations occur as those two seem not to cancel out as:
$r_{im} = \Gamma^a_{im,a} - \Gamma^a_{ai,m} = \frac 1 2 (g_{ia,ma}+g_{ma,ia}-g_{im,aa}+g_{ia,am}+g_{gg,im}-g_{ma,ai})$
while
$\Gamma_{kk,i}^m = \frac 1 2 (g_{km,ki}+g_{mk,ki}-g_{kk,mi})$
so it is not clear to me why the terms cancel out since formally they do not. I would really appreciate pointing out the reason.