For scalar fields Φ and Ψ, the Laplacian is defined by
$$∇^2(ΦΨ)=(∇^2Φ)Ψ+2∇Φ\cdot∇Ψ+Φ∇^2Ψ$$
where $∇$ is the usual del operator and $∇^2$ is the Laplacian.
How can I prove this relation? I tried the brute-force method and it became rather messy by getting an expression for the RHS and LHS. I believe that using Einstien notation would be a better approach.
Doing this with index notation is like killing a fly with a hammer, but here we go.
First task: Show that $\Delta$ is a scalar, i.e, it does not transform under coordinate transformations. In other words, it is the same regardless what basis you use. I'll leave this as an exercise.
Now, choose a (normalized) basis such that all the Christoffel symbols vanish (i.e, $\mathbf g=\operatorname{diag(1,\dots,1)}$) and such that the covariant derivative is just $$\nabla_i=\partial_i$$ Now, $$\Delta=\nabla_i\nabla^i$$ So $$\Delta(\lambda\kappa)\\ =\partial_i(\partial^i(\lambda\kappa))\\ =\partial_i(\lambda\partial^i\kappa+\kappa\partial^i\lambda)\\ =\partial_i(\lambda\partial^i\kappa)+\partial_i(\kappa\partial^i\lambda)\\ =(\partial_i\lambda)(\partial^i\kappa)+\lambda\partial_i\partial^i\kappa+(\partial_i\kappa)(\partial^i\lambda)+\kappa\partial_i\partial^i\lambda\\=\lambda\Delta\kappa+2\langle\nabla\lambda,\nabla\kappa\rangle+\kappa\Delta\lambda$$ Done.