Gradient identity involving cross product I found this particular identity in an extract of my physics course : $$\nabla \times\nabla \times\vec{F}=-\nabla^2\vec{F}+\nabla(\nabla \vec{F}) $$
However, nothing is said about $\vec{F}$, nor any physical interpretation is associated with this particular equation.
How can one prove this ? And what does this identity actually mean ?
 A: It's pretty straightforward if you use index notation.

A Brief Explanation of Index Notation
The rules for index notation are:


*

*summation is always implied when the same index occurs exactly twice in a single term.  Such indices are called summing indices or dummy indices.  Example: $\vec v \cdot \vec w=\sum_{i=1}^3v_i w_i = v_iw_i$.  Here $i$ is the only summing index.

*if there is an index which is unpaired, then it is not summed over and must appear in every term on BOTH sides of the equation.  Such indices are called free indices.  Example: $x_a = R_{ab} v_b = \sum_{b=1}^3 R_{ab} v_b$.  Here $a$ is the free index and $b$ is the summing index.  Notice here that $x_a$ is the $a$th component of the vector $\vec x$ and $R_{ab}$ is a scalar indexed by $2$ indices -- where, in general, $R_{1,2} \ne R_{2,1}$ (this is how you denote the entries of a matrix in index notation).

*if your sum includes more than $2$ of the same index then you CANNOT use this implied summation convention.  Example: $A_{ab}B_{cb}x_bv_c$ is not a valid expression using the implied summation convention because there are $3$ $b$ indices in this single term and thus the summation is ambiguous.  If this summand occurs, you'd have to explicitly use the $\sum$ symbol.


A couple of symbols also need to be explained:


*

*The permutation symbol: $\epsilon_{pqr} = \begin{cases} 1, & pqr \in \{123, 231, 312\} \\ -1, & pqr \in \{321, 213, 132\} \\ 0, & \text{otherwise}\end{cases}$

*The Kronecker delta: $\delta_{pq} = \begin{cases} 1, & p=q \\ 0, & \text{otherwise}\end{cases}$

*The partial derivative wrt the $p$th argument of the function $f(x_1, x_2, \dots)$ will be denoted $\partial_p f$.  For example, we'll denote $\dfrac{\partial f(x,y,z)}{\partial y}$ as $\partial_2 f$.



So now, let's start:
$$\begin{array}\ [\nabla \times (\nabla \times \vec F)]_p  & \text{Just consider the $p$th component}
\\ = \epsilon_{pqr}\partial_q(\nabla \times \vec F)_r & \text{$[a\times b]_p := \epsilon_{pqr}a_qb_r$ --  this can be used as}
\\ & \text{the definition of the cross product}
\\ = \epsilon_{pqr}\partial_q\epsilon_{rst}\partial_sF_t & \text{Applying the same definition}
\\ = \epsilon_{pqr}\epsilon_{rst}\partial_q\partial_sF_t & \text{$\epsilon_{rst}$ is just a scalar and so it commutes}
\\ & \text{with derivative operator}
\\ = \epsilon_{rpq}\epsilon_{rst}\partial_q\partial_sF_t &  \text{$\epsilon_{pqr}=\epsilon_{qrp}=\epsilon_{rpq}$}
\\ = (\delta_{ps}\delta_{qt} - \delta_{pt}\delta_{qs})(\partial_q\partial_sF_t) & \text{Apply the identity (which you can prove}
\\ & \text{if you want) $\epsilon_{rpq}\epsilon_{rst} = \delta_{ps}\delta_{qt} - \delta_{pt}\delta_{qs}$}
\\ = \partial_t\partial_pF_t -\partial_s\partial_sF_p & \text{Apply the identity $\delta_{qp}a_p = \delta_{pq}a_p=a_q$}
\\ = \partial_p\partial_tF_t - \partial_s\partial_sF_p & \text{Clairaut's theorem}
\\ = [\nabla (\nabla \cdot \vec F)]_p -[(\nabla \cdot \nabla)\vec F]_p & \text{We're almost done now}
\\ = [\nabla(\nabla\cdot \vec F)-\nabla^2\vec F]_p\ \ \ \ \ \ \square\end{array}$$
Once you get used to index notation, you'll see that there wasn't really any thinking required for this -- it was all pretty mechanical.
BTW, this process also shows you that $\vec F$ must have continuous second partials or else this identity doesn't hold (because we had to make use of Clairaut's theorem).
A: First of all I woul like to point out two typo: the right formula is
\begin{equation}
\nabla \times\nabla \times\vec{F}=-\nabla^2\vec{F}+\nabla(\nabla \cdot \vec{F})
\end{equation}
Now, the curl of $F=(F_1,F_2,F_3)$ can be computed by using the standard definition, i.e.
\begin{equation}
\triangledown \times F = \left(\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z}  \right) \times  \left(F_1,F_2,F_3 \right) = 
\begin{vmatrix}
i & j & k \\
\frac{\partial}{\partial x} &\frac{\partial}{\partial y} &\frac{\partial}{\partial z} \\
F_1 & F_2 & F_3
\end{vmatrix}
\end{equation}
while
\begin{equation}
\triangledown \cdot F = \left(\frac{\partial}{\partial x},\frac{\partial}{\partial y},\frac{\partial}{\partial z}  \right) \cdot  \left(F_1,F_2,F_3 \right) = \frac{\partial F_1}{\partial x}+\frac{\partial F_2}{\partial y}+\frac{\partial F_3}{\partial z}
\end{equation}
When you apply the curl to F you get a new, vector. Then apply once more the definition of curl to compute $\triangledown \times \triangledown \times F$ and use the expression of the divergence of F to simplify the computation.
A: I want to give a hint:
One may easy proof:
$$a \times (b \times c) = \left( {a \cdot c} \right){\rm{b}} - \left( {a \cdot b} \right)c$$
for any $a,b,c \in {R^3}$, and it does depent on this order to get a linear combination.
In the sense of using differential operators, one need this order:
$$\begin{array}{c}
a \times (b \times c) = {\rm{b}}\left( {a \cdot c} \right) - \left( {a \cdot b} \right)c\\
\nabla  \times (\nabla  \times \overrightarrow F ) = \nabla \left( {\nabla  \cdot \overrightarrow F } \right) - \left( {\nabla  \cdot \nabla } \right)\overrightarrow F \\
 = \nabla \left( {\nabla  \cdot \overrightarrow F } \right) - {\nabla ^2}\overrightarrow F 
\end{array}$$
