Let $\Omega \subseteq \mathbb R^d$ be open, and let $u, v: \Omega \to \mathbb R$ be harmonic functions. I now want to prove: the product function $u v: \Omega \to \mathbb R$ is harmonic on $\Omega$, if and only if $\langle \nabla u, \nabla v \rangle \equiv 0$ on $\Omega$ (where $\langle \cdot,\cdot \rangle$ is the standard scalar product).
Now to show that, I tried to write out
$$\langle \nabla u, \nabla v \rangle = \sum_{i=1}^n \frac{\partial}{\partial x_i}u \cdot \frac{\partial}{\partial x_i} v \\ = \frac{\partial}{\partial x_1} (u v) + \dots + \frac{\partial}{\partial x_n} (u v)$$
but from there on, I'm not really sure how to continue. I would need to show that this sum over the first-order derivatives is equal to $0$ if and only if the sum over the second-order derivatives
$$\Delta u v = \sum_{i=1}^n \frac{\partial^2}{\partial x_n^2} (u v)$$
is also equal to $0$ everywhere. But I'm not sure how I can show that and how the harmonicity of $u, v$ themselves comes into play? Am I missing something obvious here?