# Calculus/Vector calculus identities

In this chapter, numerous identities related to the gradient ($\nabla f$), directional derivative ($(\mathbf {V} \cdot \nabla )f$, $(\mathbf {V} \cdot \nabla )\mathbf {F}$), divergence ($\nabla \cdot \mathbf {F}$), Laplacian ($\nabla ^{2}f$, $\nabla ^{2}\mathbf {F}$), and curl ($\nabla \times \mathbf {F}$) will be derived.

## Notation

To simplify the derivation of various vector identities, the following notation will be utilized:
* The coordinates $x,y,z$ will instead be denoted with $x_{1},x_{2},x_{3}$ respectively.
* Given an arbitrary vector $\mathbf {F}$, then $F_{i}$ will denote the $i^{\text{th}}$ entry of $\mathbf {F}$ where $i=1,2,3$. All vectors will be assumed to be denoted by Cartesian basis vectors ($\mathbf {i} ,\mathbf {j} ,\mathbf {k}$) unless otherwise specified: $\mathbf {F} =F_{1}\mathbf {i} +F_{2}\mathbf {j} +F_{3}\mathbf {k}$.
* Given an arbitrary expression $f:\{1,2,3\}\to \mathbb {R}$ that assigns a real number to each index $i=1,2,3$, then $(i,f(i))$ will denote the vector whose entries are determined by $f$. For example, $\mathbf {F} =(i,F_{i})$.
* Given an arbitrary expression $f:\{1,2,3\}\to \mathbb {R}$ that assigns a real number to each index $i=1,2,3$, then $\sum _{i}f(i)$ will denote the sum $f(1)+f(2)+f(3)$. For example, $\nabla \cdot \mathbf {F} =\sum _{i}{\frac {\partial F_{i}}{\partial x_{i}}}$.
* Given an index variable $i\in \{1,2,3\}$, $i+1$ will rotate $i$ forwards by 1, and $i+2$ will rotate $i$ forwards by 2. In essence, $i+1=\left\{{\begin{array}{cc}i+1&(i=1,2)\\1&(i=3)\end{array}}\right.$ and $i+2=\left\{{\begin{array}{cc}3&(i=1)\\i-1&(i=2,3)\end{array}}\right.$. For example, $\mathbf {F} \times \mathbf {G} =(i,F_{i+1}G_{i+2}-F_{i+2}G_{i+1})$.
As an example of using the above notation, consider the problem of expanding the triple cross product $\mathbf {F} \times (\mathbf {G} \times \mathbf {H} )$.
$\mathbf {F} \times (\mathbf {G} \times \mathbf {H} )=\mathbf {F} \times (i,G_{i+1}H_{i+2}-G_{i+2}H_{i+1})$
$=(i,F_{i+1}(G_{i}H_{i+1}-G_{i+1}H_{i})-F_{i+2}(G_{i+2}H_{i}-G_{i}H_{i+2}))$
$=(i,G_{i}(F_{i+1}H_{i+1}+F_{i+2}H_{i+2})-(F_{i+1}G_{i+1}+F_{i+2}G_{i+2})H_{i})$
$=(i,G_{i}(F_{i}H_{i}+F_{i+1}H_{i+1}+F_{i+2}H_{i+2})-(F_{i}G_{i}+F_{i+1}G_{i+1}+F_{i+2}G_{i+2})H_{i})$
$=(i,G_{i}(\mathbf {F} \cdot \mathbf {H} )-(\mathbf {F} \cdot \mathbf {G} )H_{i})$
$=(\mathbf {F} \cdot \mathbf {H} )\mathbf {G} -(\mathbf {F} \cdot \mathbf {G} )\mathbf {H}$
Therefore: $\mathbf {F} \times (\mathbf {G} \times \mathbf {H} )=(\mathbf {F} \cdot \mathbf {H} )\mathbf {G} -(\mathbf {F} \cdot \mathbf {G} )\mathbf {H}$
As another example of using the above notation, consider the scalar triple product $\mathbf {F} \cdot (\mathbf {G} \times \mathbf {H} )$
$\mathbf {F} \cdot (\mathbf {G} \times \mathbf {H} )=\mathbf {F} \cdot (i,G_{i+1}H_{i+2}-G_{i+2}H_{i+1})$
$=\sum _{i}F_{i}(G_{i+1}H_{i+2}-G_{i+2}H_{i+1})$
$=(\sum _{i}F_{i}G_{i+1}H_{i+2})-(\sum _{i}F_{i}G_{i+2}H_{i+1})$
The index $i$ in the above summations can be shifted by fixed amounts without changing the sum. For example, $\sum _{i}F_{i}G_{i+1}H_{i+2}=\sum _{i}F_{i+1}G_{i+2}H_{i}=\sum _{i}F_{i+2}G_{i}H_{i+1}$. This allows:
$(\sum _{i}F_{i}G_{i+1}H_{i+2})-(\sum _{i}F_{i}G_{i+2}H_{i+1})=(\sum _{i}F_{i+2}G_{i}H_{i+1})-(\sum _{i}F_{i+1}G_{i}H_{i+2})=(\sum _{i}F_{i+1}G_{i+2}H_{i})-(\sum _{i}F_{i+2}G_{i+1}H_{i})$
$\implies \mathbf {F} \cdot (i,G_{i+1}H_{i+2}-G_{i+2}H_{i+1})=\mathbf {G} \cdot (i,H_{i+1}F_{i+2}-H_{i+2}F_{i+1})=\mathbf {H} \cdot (i,F_{i+1}G_{i+2}-F_{i+2}G_{i+1})$
$\implies \mathbf {F} \cdot (\mathbf {G} \times \mathbf {H} )=\mathbf {G} \cdot (\mathbf {H} \times \mathbf {F} )=\mathbf {H} \cdot (\mathbf {F} \times \mathbf {G} )$
which establishes the cyclical property of the scalar triple product.

## Gradient Identities

Given scalar fields, $f$ and $g$, then $\nabla (f+g)=(\nabla f)+(\nabla g)$.
$\nabla (f+g)=(i,{\frac {\partial }{\partial x_{i}}}(f+g))$ $=(i,{\frac {\partial f}{\partial x_{i}}}+{\frac {\partial g}{\partial x_{i}}})$ $=(i,{\frac {\partial f}{\partial x_{i}}})+(i,{\frac {\partial g}{\partial x_{i}}})$ $=(\nabla f)+(\nabla g)$
Given scalar fields $f$ and $g$, then $\nabla (fg)=(\nabla f)g+f(\nabla g)$. If $f$ is a constant $c$, then $\nabla (cg)=c(\nabla g)$.
$\nabla (fg)=(i,{\frac {\partial }{\partial x_{i}}}(fg))$ $=(i,{\frac {\partial f}{\partial x_{i}}}g+f{\frac {\partial g}{\partial x_{i}}})$ $=(i,{\frac {\partial f}{\partial x_{i}}})g+f(i,{\frac {\partial g}{\partial x_{i}}})$ $=(\nabla f)g+f(\nabla g)$
Given vector fields $\mathbf {F}$ and $\mathbf {G}$, then $\nabla (\mathbf {F} \cdot \mathbf {G} )=((\mathbf {F} \cdot \nabla )\mathbf {G} +\mathbf {F} \times (\nabla \times \mathbf {G} ))+((\mathbf {G} \cdot \nabla )\mathbf {F} +\mathbf {G} \times (\nabla \times \mathbf {F} ))$
$\nabla (\mathbf {F} \cdot \mathbf {G} )=(i,{\frac {\partial }{\partial x_{i}}}(\mathbf {F} \cdot \mathbf {G} ))$ $=(i,{\frac {\partial }{\partial x_{i}}}(\sum _{j}F_{j}G_{j}))$ $=(i,\sum _{j}({\frac {\partial F_{j}}{\partial x_{i}}}G_{j}+F_{j}{\frac {\partial G_{j}}{\partial x_{i}}}))$ $=(i,\sum _{j}F_{j}{\frac {\partial G_{j}}{\partial x_{i}}})+(i,\sum _{j}G_{j}{\frac {\partial F_{j}}{\partial x_{i}}})$
$=(i,F_{i}{\frac {\partial G_{i}}{\partial x_{i}}}+F_{i+1}{\frac {\partial G_{i+1}}{\partial x_{i}}}+F_{i+2}{\frac {\partial G_{i+2}}{\partial x_{i}}})$ $+(i,G_{i}{\frac {\partial F_{i}}{\partial x_{i}}}+G_{i+1}{\frac {\partial F_{i+1}}{\partial x_{i}}}+G_{i+2}{\frac {\partial F_{i+2}}{\partial x_{i}}})$
$=(i,(F_{i}{\frac {\partial G_{i}}{\partial x_{i}}}+F_{i+1}{\frac {\partial G_{i}}{\partial x_{i+1}}}+F_{i+2}{\frac {\partial G_{i}}{\partial x_{i+2}}})+((F_{i+1}{\frac {\partial G_{i+1}}{\partial x_{i}}}-F_{i+1}{\frac {\partial G_{i}}{\partial x_{i+1}}})+(F_{i+2}{\frac {\partial G_{i+2}}{\partial x_{i}}}-F_{i+2}{\frac {\partial G_{i}}{\partial x_{i+2}}})))$ $+(i,(G_{i}{\frac {\partial F_{i}}{\partial x_{i}}}+G_{i+1}{\frac {\partial F_{i}}{\partial x_{i+1}}}+G_{i+2}{\frac {\partial F_{i}}{\partial x_{i+2}}})+((G_{i+1}{\frac {\partial F_{i+1}}{\partial x_{i}}}-G_{i+1}{\frac {\partial F_{i}}{\partial x_{i+1}}})+(G_{i+2}{\frac {\partial F_{i+2}}{\partial x_{i}}}-G_{i+2}{\frac {\partial F_{i}}{\partial x_{i+2}}})))$
$=(i,\sum _{j}F_{j}{\frac {\partial G_{i}}{\partial x_{j}}})+(i,F_{i+1}({\frac {\partial G_{i+1}}{\partial x_{i}}}-{\frac {\partial G_{i}}{\partial x_{i+1}}})-F_{i+2}({\frac {\partial G_{i}}{\partial x_{i+2}}}-{\frac {\partial G_{i+2}}{\partial x_{i}}}))$ $+(i,\sum _{j}G_{j}{\frac {\partial F_{i}}{\partial x_{j}}})+(i,G_{i+1}({\frac {\partial F_{i+1}}{\partial x_{i}}}-{\frac {\partial F_{i}}{\partial x_{i+1}}})-G_{i+2}({\frac {\partial F_{i}}{\partial x_{i+2}}}-{\frac {\partial F_{i+2}}{\partial x_{i}}}))$
$=(i,(\mathbf {F} \cdot \nabla )G_{i})+\mathbf {F} \times (i,{\frac {\partial G_{i+2}}{\partial x_{i+1}}}-{\frac {\partial G_{i+1}}{\partial x_{i+2}}})+(i,(\mathbf {G} \cdot \nabla )F_{i})+\mathbf {G} \times (i,{\frac {\partial F_{i+2}}{\partial x_{i+1}}}-{\frac {\partial F_{i+1}}{\partial x_{i+2}}})$
$=((\mathbf {F} \cdot \nabla )\mathbf {G} +\mathbf {F} \times (\nabla \times \mathbf {G} ))+((\mathbf {G} \cdot \nabla )\mathbf {F} +\mathbf {G} \times (\nabla \times \mathbf {F} ))$
Given scalar fields $f_{1},f_{2},\dots ,f_{n}$ and an $n$ input function $g(y_{1},y_{2},\dots ,y_{n})$, then $\nabla (g(f_{1},f_{2},\dots ,f_{n}))={\frac {\partial g}{\partial y_{1}}}{\bigg |}_{y_{1}=f_{1}}(\nabla f_{1})+{\frac {\partial g}{\partial y_{2}}}{\bigg |}_{y_{2}=f_{2}}(\nabla f_{2})+\dots +{\frac {\partial g}{\partial y_{n}}}{\bigg |}_{y_{n}=f_{n}}(\nabla f_{n})$.
$\nabla (g(f_{1},f_{2},\dots ,f_{n}))=(i,{\frac {\partial }{\partial x_{i}}}(g(f_{1},f_{2},\dots ,f_{n})))$ $=(i,{\frac {\partial g}{\partial y_{1}}}{\bigg |}_{y_{1}=f_{1}}{\frac {\partial f_{1}}{\partial x_{i}}}+{\frac {\partial g}{\partial y_{2}}}{\bigg |}_{y_{2}=f_{2}}{\frac {\partial f_{2}}{\partial x_{i}}}+\dots +{\frac {\partial g}{\partial y_{n}}}{\bigg |}_{y_{n}=f_{n}}{\frac {\partial f_{n}}{\partial x_{i}}})$ $={\frac {\partial g}{\partial y_{1}}}{\bigg |}_{y_{1}=f_{1}}(i,{\frac {\partial f_{1}}{\partial x_{i}}})+{\frac {\partial g}{\partial y_{2}}}{\bigg |}_{y_{2}=f_{2}}(i,{\frac {\partial f_{2}}{\partial x_{i}}})+\dots +{\frac {\partial g}{\partial y_{n}}}{\bigg |}_{y_{n}=f_{n}}(i,{\frac {\partial f_{n}}{\partial x_{i}}})$ $={\frac {\partial g}{\partial y_{1}}}{\bigg |}_{y_{1}=f_{1}}(\nabla f_{1})+{\frac {\partial g}{\partial y_{2}}}{\bigg |}_{y_{2}=f_{2}}(\nabla f_{2})+\dots +{\frac {\partial g}{\partial y_{n}}}{\bigg |}_{y_{n}=f_{n}}(\nabla f_{n})$

## Directional Derivative Identities

Given vector fields $\mathbf {V}$ and $\mathbf {W}$, and scalar field $f$, then $((\mathbf {V} +\mathbf {W} )\cdot \nabla )f=(\mathbf {V} \cdot \nabla )f+(\mathbf {W} \cdot \nabla )f$.
When $\mathbf {F}$ is a vector field, it is also the case that: $((\mathbf {V} +\mathbf {W} )\cdot \nabla )\mathbf {F} =(\mathbf {V} \cdot \nabla )\mathbf {F} +(\mathbf {W} \cdot \nabla )\mathbf {F}$.
For scalar fields:
$((\mathbf {V} +\mathbf {W} )\cdot \nabla )f=\sum _{i}((V_{i}+W_{i}){\frac {\partial f}{\partial x_{i}}})$ $=\sum _{i}(V_{i}{\frac {\partial f}{\partial x_{i}}}+W_{i}{\frac {\partial f}{\partial x_{i}}})$ $=\sum _{i}(V_{i}{\frac {\partial f}{\partial x_{i}}})+\sum _{i}(W_{i}{\frac {\partial f}{\partial x_{i}}})$ $=(\mathbf {V} \cdot \nabla )f+(\mathbf {W} \cdot \nabla )f$
For vector fields:
$((\mathbf {V} +\mathbf {W} )\cdot \nabla )\mathbf {F} =(i,((\mathbf {V} +\mathbf {W} )\cdot \nabla )F_{i})$ $=(i,(\mathbf {V} \cdot \nabla )F_{i}+(\mathbf {W} \cdot \nabla )F_{i})$ $=(\mathbf {V} \cdot \nabla )\mathbf {F} +(\mathbf {W} \cdot \nabla )\mathbf {F}$
Given vector field $\mathbf {V}$, and scalar fields $v$ and $f$, then $((v\mathbf {V} )\cdot \nabla )f=v((\mathbf {V} \cdot \nabla )f)$.
When $\mathbf {F}$ is a vector field, it is also the case that: $((v\mathbf {V} )\cdot \nabla )\mathbf {F} =v((\mathbf {V} \cdot \nabla )\mathbf {F} )$.
For scalar fields:
$((v\mathbf {V} )\cdot \nabla )f=\sum _{i}(vV_{i}{\frac {\partial f}{\partial x_{i}}})$ $=v\sum _{i}(V_{i}{\frac {\partial f}{\partial x_{i}}})$ $=v((\mathbf {V} \cdot \nabla )f)$
For vector fields:
$((v\mathbf {V} )\cdot \nabla )\mathbf {F} =(i,((v\mathbf {V} )\cdot \nabla )F_{i})$ $=(i,v((\mathbf {V} \cdot \nabla )F_{i}))$ $=v((\mathbf {V} \cdot \nabla )\mathbf {F} )$
Given vector field $\mathbf {V}$, and scalar fields $f$ and $g$, then $(\mathbf {V} \cdot \nabla )(f+g)=(\mathbf {V} \cdot \nabla )f+(\mathbf {V} \cdot \nabla )g$.
When $\mathbf {F}$ and $\mathbf {G}$ are vector fields, it is also the case that: $(\mathbf {V} \cdot \nabla )(\mathbf {F} +\mathbf {G} )=(\mathbf {V} \cdot \nabla )\mathbf {F} +(\mathbf {V} \cdot \nabla )\mathbf {G}$.
For scalar fields:
$(\mathbf {V} \cdot \nabla )(f+g)=\sum _{i}(V_{i}{\frac {\partial }{\partial x_{i}}}(f+g))$ $=\sum _{i}(V_{i}{\frac {\partial f}{\partial x_{i}}}+V_{i}{\frac {\partial g}{\partial x_{i}}})$ $=\sum _{i}(V_{i}{\frac {\partial f}{\partial x_{i}}})+\sum _{i}(V_{i}{\frac {\partial g}{\partial x_{i}}})$ $=(\mathbf {V} \cdot \nabla )f+(\mathbf {V} \cdot \nabla )g$
For vector fields:
$(\mathbf {V} \cdot \nabla )(\mathbf {F} +\mathbf {G} )=(i,(\mathbf {V} \cdot \nabla )(F_{i}+G_{i}))$ $=(i,(\mathbf {V} \cdot \nabla )F_{i}+(\mathbf {V} \cdot \nabla )G_{i})$ $=(\mathbf {V} \cdot \nabla )\mathbf {F} +(\mathbf {V} \cdot \nabla )\mathbf {G}$
Given vector field $\mathbf {V}$, and scalar fields $f$ and $g$, then $(\mathbf {V} \cdot \nabla )(fg)=((\mathbf {V} \cdot \nabla )f)g+f((\mathbf {V} \cdot \nabla )g)$
If $\mathbf {G}$ is a vector field, it is also the case that: $(\mathbf {V} \cdot \nabla )(f\mathbf {G} )=((\mathbf {V} \cdot \nabla )f)\mathbf {G} +f((\mathbf {V} \cdot \nabla )\mathbf {G} )$
For scalar fields:
$(\mathbf {V} \cdot \nabla )(fg)=\sum _{i}V_{i}{\frac {\partial }{\partial x_{i}}}(fg)$ $=\sum _{i}V_{i}({\frac {\partial f}{\partial x_{i}}}g+f{\frac {\partial g}{\partial x_{i}}})$ $=(\sum _{i}V_{i}{\frac {\partial f}{\partial x_{i}}})g+f(\sum _{i}V_{i}{\frac {\partial g}{\partial x_{i}}})$ $=((\mathbf {V} \cdot \nabla )f)g+f((\mathbf {V} \cdot \nabla )g)$
For vector fields:
$(\mathbf {V} \cdot \nabla )(f\mathbf {G} )=(i,(\mathbf {V} \cdot \nabla )(fG_{i}))$ $=(i,((\mathbf {V} \cdot \nabla )f)G_{i}+f((\mathbf {V} \cdot \nabla )G_{i}))$ $=((\mathbf {V} \cdot \nabla )f)\mathbf {G} +f((\mathbf {V} \cdot \nabla )\mathbf {G} )$
Given vector fields $\mathbf {V}$, $\mathbf {F}$, and $\mathbf {G}$, then $(\mathbf {V} \cdot \nabla )(\mathbf {F} \cdot \mathbf {G} )=((\mathbf {V} \cdot \nabla )\mathbf {F} )\cdot \mathbf {G} +\mathbf {F} \cdot ((\mathbf {V} \cdot \nabla )\mathbf {G} )$
$(\mathbf {V} \cdot \nabla )(\mathbf {F} \cdot \mathbf {G} )=\sum _{i}V_{i}{\frac {\partial }{\partial x_{i}}}(\mathbf {F} \cdot \mathbf {G} )$ $=\sum _{i}V_{i}{\frac {\partial }{\partial x_{i}}}\sum _{j}(F_{j}G_{j})$ $=\sum _{i}\sum _{j}V_{i}{\frac {\partial }{\partial x_{i}}}(F_{j}G_{j})$ $=\sum _{i}\sum _{j}V_{i}({\frac {\partial F_{j}}{\partial x_{i}}}G_{j}+F_{j}{\frac {\partial G_{j}}{\partial x_{i}}})$ $=\sum _{j}((\sum _{i}V_{i}{\frac {\partial F_{j}}{\partial x_{i}}})G_{j})+\sum _{j}(F_{j}(\sum _{i}V_{i}{\frac {\partial G_{j}}{\partial x_{i}}}))$ $=\sum _{j}(((\mathbf {V} \cdot \nabla )F_{j})G_{j})+\sum _{j}(F_{j}((\mathbf {V} \cdot \nabla )G_{j}))$ $=((\mathbf {V} \cdot \nabla )\mathbf {F} )\cdot \mathbf {G} +\mathbf {F} \cdot ((\mathbf {V} \cdot \nabla )\mathbf {G} )$
Given vector fields $\mathbf {V}$, $\mathbf {F}$, and $\mathbf {G}$, then $(\mathbf {V} \cdot \nabla )(\mathbf {F} \times \mathbf {G} )=((\mathbf {V} \cdot \nabla )\mathbf {F} )\times \mathbf {G} +\mathbf {F} \times ((\mathbf {V} \cdot \nabla )\mathbf {G} )$
$(\mathbf {V} \cdot \nabla )(\mathbf {F} \times \mathbf {G} )=(i,(\mathbf {V} \cdot \nabla )(F_{i+1}G_{i+2}-F_{i+2}G_{i+1}))$ $=(i,\sum _{j}V_{j}{\frac {\partial }{\partial x_{j}}}(F_{i+1}G_{i+2}-F_{i+2}G_{i+1}))$ $=(i,\sum _{j}V_{j}(({\frac {\partial F_{i+1}}{\partial x_{j}}}G_{i+2}+F_{i+1}{\frac {\partial G_{i+2}}{\partial x_{j}}})-({\frac {\partial F_{i+2}}{\partial x_{j}}}G_{i+1}+F_{i+2}{\frac {\partial G_{i+1}}{\partial x_{j}}})))$ $=(i,(\sum _{j}V_{j}{\frac {\partial F_{i+1}}{\partial x_{j}}})G_{i+2}-(\sum _{j}V_{j}{\frac {\partial F_{i+2}}{\partial x_{j}}})G_{i+1})+(i,F_{i+1}(\sum _{j}V_{j}{\frac {\partial G_{i+2}}{\partial x_{j}}})-F_{i+2}(\sum _{j}V_{j}{\frac {\partial G_{i+1}}{\partial x_{j}}}))$ $=(i,((\mathbf {V} \cdot \nabla )F_{i+1})G_{i+2}-((\mathbf {V} \cdot \nabla )F_{i+2})G_{i+1})+(i,F_{i+1}((\mathbf {V} \cdot \nabla )G_{i+2})-F_{i+2}((\mathbf {V} \cdot \nabla )G_{i+1}))$ $=((\mathbf {V} \cdot \nabla )\mathbf {F} )\times \mathbf {G} +\mathbf {F} \times ((\mathbf {V} \cdot \nabla )\mathbf {G} )$

## Divergence Identities

Given vector fields $\mathbf {F}$ and $\mathbf {G}$, then $\nabla \cdot (\mathbf {F} +\mathbf {G} )=(\nabla \cdot \mathbf {F} )+(\nabla \cdot \mathbf {G} )$.
$\nabla \cdot (\mathbf {F} +\mathbf {G} )=\sum _{i}({\frac {\partial }{\partial x_{i}}}(F_{i}+G_{i}))$ $=(\sum _{i}{\frac {\partial F_{i}}{\partial x_{i}}})+(\sum _{i}{\frac {\partial G_{i}}{\partial x_{i}}})$ $=(\nabla \cdot \mathbf {F} )+(\nabla \cdot \mathbf {G} )$
Given a scalar field $f$ and a vector field $\mathbf {G}$, then $\nabla \cdot (f\mathbf {G} )=(\nabla f)\cdot \mathbf {G} +f(\nabla \cdot \mathbf {G} )$. If $f$ is a constant $c$, then $\nabla \cdot (c\mathbf {G} )=c(\nabla \cdot \mathbf {G} )$. If $\mathbf {G}$ is a constant $\mathbf {C}$, then $\nabla \cdot (f\mathbf {C} )=(\nabla f)\cdot \mathbf {C}$.
$\nabla \cdot (f\mathbf {G} )=\sum _{i}{\frac {\partial }{\partial x_{i}}}(fG_{i})$ $=\sum _{i}({\frac {\partial f}{\partial x_{i}}}G_{i}+f{\frac {\partial G_{i}}{\partial x_{i}}})$ $=\sum _{i}({\frac {\partial f}{\partial x_{i}}}G_{i})+f\sum _{i}{\frac {\partial G_{i}}{\partial x_{i}}}$ $=(\nabla f)\cdot \mathbf {G} +f(\nabla \cdot \mathbf {G} )$
Given vector fields $\mathbf {F}$ and $\mathbf {G}$, then $\nabla \cdot (\mathbf {F} \times \mathbf {G} )=(\nabla \times \mathbf {F} )\cdot \mathbf {G} -\mathbf {F} \cdot (\nabla \times \mathbf {G} )$.
$\nabla \cdot (\mathbf {F} \times \mathbf {G} )=\sum _{i}{\frac {\partial }{\partial x_{i}}}(F_{i+1}G_{i+2}-F_{i+2}G_{i+1})$ $=\sum _{i}(({\frac {\partial F_{i+1}}{\partial x_{i}}}G_{i+2}+F_{i+1}{\frac {\partial G_{i+2}}{\partial x_{i}}})-({\frac {\partial F_{i+2}}{\partial x_{i}}}G_{i+1}+F_{i+2}{\frac {\partial G_{i+1}}{\partial x_{i}}}))$ $=\sum _{i}(({\frac {\partial F_{i+2}}{\partial x_{i+1}}}G_{i}+F_{i}{\frac {\partial G_{i+1}}{\partial x_{i+2}}})-({\frac {\partial F_{i+1}}{\partial x_{i+2}}}G_{i}+F_{i}{\frac {\partial G_{i+2}}{\partial x_{i+1}}}))$ $=\sum _{i}(({\frac {\partial F_{i+2}}{\partial x_{i+1}}}-{\frac {\partial F_{i+1}}{\partial x_{i+2}}})G_{i}-F_{i}({\frac {\partial G_{i+2}}{\partial x_{i+1}}}-{\frac {\partial G_{i+1}}{\partial x_{i+2}}}))$ $=\sum _{i}(\nabla \times \mathbf {F} )_{i}G_{i}-\sum _{i}F_{i}(\nabla \times \mathbf {G} )_{i}$ $=(\nabla \times \mathbf {F} )\cdot \mathbf {G} -\mathbf {F} \cdot (\nabla \times \mathbf {G} )$
In the above derivation, the third equality is established by cycling the terms inside a sum. For example: $\sum _{i}{\frac {\partial F_{i+1}}{\partial x_{i}}}G_{i+2}=\sum _{i}{\frac {\partial F_{i+2}}{\partial x_{i+1}}}G_{i}$ by replacing $i$ with $i+1$. Different terms can be cycled independently: $\sum _{i}({\frac {\partial F_{i+1}}{\partial x_{i}}}G_{i+2}+F_{i+1}{\frac {\partial G_{i+2}}{\partial x_{i}}})=\sum _{i}({\frac {\partial F_{i+2}}{\partial x_{i+1}}}G_{i}+F_{i}{\frac {\partial G_{i+1}}{\partial x_{i+2}}})$
The following identity is a very important property regarding vector fields which are the curl of another vector field. A vector field which is the curl of another vector field is divergence free. Given vector field $\mathbf {F}$, then $\nabla \cdot (\nabla \times \mathbf {F} )=0$
$\nabla \cdot (\nabla \times \mathbf {F} )=\nabla \cdot (i,{\frac {\partial F_{i+2}}{\partial x_{i+1}}}-{\frac {\partial F_{i+1}}{\partial x_{i+2}}})$ $=\sum _{i}{\frac {\partial }{\partial x_{i}}}({\frac {\partial F_{i+2}}{\partial x_{i+1}}}-{\frac {\partial F_{i+1}}{\partial x_{i+2}}})$ $=\sum _{i}({\frac {\partial ^{2}F_{i+2}}{\partial x_{i}\partial x_{i+1}}}-{\frac {\partial ^{2}F_{i+1}}{\partial x_{i}\partial x_{i+2}}})$ $=\sum _{i}{\frac {\partial ^{2}F_{i+2}}{\partial x_{i}\partial x_{i+1}}}-\sum _{i}{\frac {\partial ^{2}F_{i+1}}{\partial x_{i+2}\partial x_{i}}}$ $=\sum _{i}{\frac {\partial ^{2}F_{i+2}}{\partial x_{i}\partial x_{i+1}}}-\sum _{i}{\frac {\partial ^{2}F_{i+2}}{\partial x_{i}\partial x_{i+1}}}$ $=0$

## Laplacian Identities

Given scalar fields $f$ and $g$, then $\nabla ^{2}(f+g)=(\nabla ^{2}f)+(\nabla ^{2}g)$
When $\mathbf {F}$ and $\mathbf {G}$ are vector fields, it is also the case that: $\nabla ^{2}(\mathbf {F} +\mathbf {G} )=(\nabla ^{2}\mathbf {F} )+(\nabla ^{2}\mathbf {G} )$
For scalar fields:
$\nabla ^{2}(f+g)=\sum _{i}{\frac {\partial ^{2}}{\partial x_{i}^{2}}}(f+g)$ $=\sum _{i}({\frac {\partial ^{2}f}{\partial x_{i}^{2}}}+{\frac {\partial ^{2}g}{\partial x_{i}^{2}}})$ $=(\sum _{i}{\frac {\partial ^{2}f}{\partial x_{i}^{2}}})+(\sum _{i}{\frac {\partial ^{2}g}{\partial x_{i}^{2}}})$ $=(\nabla ^{2}f)+(\nabla ^{2}g)$
For vector fields:
$\nabla ^{2}(\mathbf {F} +\mathbf {G} )=(i,\nabla ^{2}(F_{i}+G_{i}))$ $=(i,(\nabla ^{2}F_{i})+(\nabla ^{2}G_{i}))$ $=(\nabla ^{2}\mathbf {F} )+(\nabla ^{2}\mathbf {G} )$
Given scalar fields $f$ and $g$, then $\nabla ^{2}(fg)=(\nabla ^{2}f)g+2(\nabla f)\cdot (\nabla g)+f(\nabla ^{2}g)$
When $\mathbf {G}$ is a vector field, it is also the case that $\nabla ^{2}(f\mathbf {G} )=(\nabla ^{2}f)\mathbf {G} +2((\nabla f)\cdot \nabla )\mathbf {G} +f(\nabla ^{2}\mathbf {G} )$
For scalar fields:
$\nabla ^{2}(fg)=\sum _{i}{\frac {\partial ^{2}}{\partial x_{i}^{2}}}(fg)$ $=\sum _{i}{\frac {\partial }{\partial x_{i}}}({\frac {\partial f}{\partial x_{i}}}g+f{\frac {\partial g}{\partial x_{i}}})$ $=\sum _{i}({\frac {\partial ^{2}f}{\partial x_{i}^{2}}}g+2{\frac {\partial f}{\partial x_{i}}}{\frac {\partial g}{\partial x_{i}}}+f{\frac {\partial ^{2}g}{\partial x_{i}^{2}}})$ $=(\sum _{i}{\frac {\partial ^{2}f}{\partial x_{i}^{2}}})g+2\sum _{i}({\frac {\partial f}{\partial x_{i}}}{\frac {\partial g}{\partial x_{i}}})+f(\sum _{i}{\frac {\partial ^{2}g}{\partial x_{i}^{2}}})$ $=(\nabla ^{2}f)g+2(\nabla f)\cdot (\nabla g)+f(\nabla ^{2}g)$
For vector fields:
$\nabla ^{2}(f\mathbf {G} )=(i,\nabla ^{2}(fG_{i}))$ $=(i,(\nabla ^{2}f)G_{i}+2(\nabla f)\cdot (\nabla G_{i})+f(\nabla ^{2}G_{i}))$ $=(i,(\nabla ^{2}f)G_{i})+2(i,((\nabla f)\cdot \nabla )G_{i})+(i,f(\nabla ^{2}G_{i}))$ $=(\nabla ^{2}f)\mathbf {G} +2((\nabla f)\cdot \nabla )\mathbf {G} +f(\nabla ^{2}\mathbf {G} )$

## Curl Identities

Given vector fields $\mathbf {F}$ and $\mathbf {G}$, then $\nabla \times (\mathbf {F} +\mathbf {G} )=(\nabla \times \mathbf {F} )+(\nabla \times \mathbf {G} )$
$\nabla \times (\mathbf {F} +\mathbf {G} )=(i,{\frac {\partial }{\partial x_{i+1}}}(F_{i+2}+G_{i+2})-{\frac {\partial }{\partial x_{i+2}}}(F_{i+1}+G_{i+1}))$ $=(i,({\frac {\partial F_{i+2}}{\partial x_{i+1}}}+{\frac {\partial G_{i+2}}{\partial x_{i+1}}})-({\frac {\partial F_{i+1}}{\partial x_{i+2}}}+{\frac {\partial G_{i+1}}{\partial x_{i+2}}}))$ $=(i,{\frac {\partial F_{i+2}}{\partial x_{i+1}}}-{\frac {\partial F_{i+1}}{\partial x_{i+2}}})+(i,{\frac {\partial G_{i+2}}{\partial x_{i+1}}}-{\frac {\partial G_{i+1}}{\partial x_{i+2}}})$ $=(\nabla \times \mathbf {F} )+(\nabla \times \mathbf {G} )$
Given scalar field $f$ and vector field $\mathbf {G}$, then $\nabla \times (f\mathbf {G} )=(\nabla f)\times \mathbf {G} +f(\nabla \times \mathbf {G} )$. If $f$ is a constant $c$, then $\nabla \times (c\mathbf {G} )=c(\nabla \times \mathbf {G} )$. If $\mathbf {G}$ is a constant $\mathbf {C}$, then $\nabla \times (f\mathbf {C} )=(\nabla f)\times \mathbf {C}$.
$\nabla \times (f\mathbf {G} )=(i,{\frac {\partial }{\partial x_{i+1}}}(fG_{i+2})-{\frac {\partial }{\partial x_{i+2}}}(fG_{i+1}))$ $=(i,({\frac {\partial f}{\partial x_{i+1}}}G_{i+2}+f{\frac {\partial G_{i+2}}{\partial x_{i+1}}})-({\frac {\partial f}{\partial x_{i+2}}}G_{i+1}+f{\frac {\partial G_{i+1}}{\partial x_{i+2}}}))$ $=(i,{\frac {\partial f}{\partial x_{i+1}}}G_{i+2}-{\frac {\partial f}{\partial x_{i+2}}}G_{i+1})+f(i,{\frac {\partial G_{i+2}}{\partial x_{i+1}}}-{\frac {\partial G_{i+1}}{\partial x_{i+2}}})$ $=(\nabla f)\times \mathbf {G} +f(\nabla \times \mathbf {G} )$
Given vector fields $\mathbf {F}$ and $\mathbf {G}$, then $\nabla \times (\mathbf {F} \times \mathbf {G} )=((\nabla \cdot \mathbf {G} )\mathbf {F} +(\mathbf {G} \cdot \nabla )\mathbf {F} )-((\nabla \cdot \mathbf {F} )\mathbf {G} +(\mathbf {F} \cdot \nabla )\mathbf {G} )$
$\nabla \times (\mathbf {F} \times \mathbf {G} )=\nabla \times (i,F_{i+1}G_{i+2}-F_{i+2}G_{i+1})$ $=(i,{\frac {\partial }{\partial x_{i+1}}}(F_{i}G_{i+1}-F_{i+1}G_{i})-{\frac {\partial }{\partial x_{i+2}}}(F_{i+2}G_{i}-F_{i}G_{i+2}))$ $=(i,(({\frac {\partial F_{i}}{\partial x_{i+1}}}G_{i+1}+F_{i}{\frac {\partial G_{i+1}}{\partial x_{i+1}}})-({\frac {\partial F_{i+1}}{\partial x_{i+1}}}G_{i}+F_{i+1}{\frac {\partial G_{i}}{\partial x_{i+1}}}))-(({\frac {\partial F_{i+2}}{\partial x_{i+2}}}G_{i}+F_{i+2}{\frac {\partial G_{i}}{\partial x_{i+2}}})-({\frac {\partial F_{i}}{\partial x_{i+2}}}G_{i+2}+F_{i}{\frac {\partial G_{i+2}}{\partial x_{i+2}}})))$ $=(i,F_{i}({\frac {\partial G_{i+1}}{\partial x_{i+1}}}+{\frac {\partial G_{i+2}}{\partial x_{i+2}}})-({\frac {\partial F_{i+1}}{\partial x_{i+1}}}+{\frac {\partial F_{i+2}}{\partial x_{i+2}}})G_{i}-(F_{i+1}{\frac {\partial G_{i}}{\partial x_{i+1}}}+F_{i+2}{\frac {\partial G_{i}}{\partial x_{i+2}}})+({\frac {\partial F_{i}}{\partial x_{i+1}}}G_{i+1}+{\frac {\partial F_{i}}{\partial x_{i+2}}}G_{i+2}))$
$=(i,F_{i}({\frac {\partial G_{i}}{\partial x_{i}}}+{\frac {\partial G_{i+1}}{\partial x_{i+1}}}+{\frac {\partial G_{i+2}}{\partial x_{i+2}}})-({\frac {\partial F_{i}}{\partial x_{i}}}+{\frac {\partial F_{i+1}}{\partial x_{i+1}}}+{\frac {\partial F_{i+2}}{\partial x_{i+2}}})G_{i}$ $-(F_{i}{\frac {\partial G_{i}}{\partial x_{i}}}+F_{i+1}{\frac {\partial G_{i}}{\partial x_{i+1}}}+F_{i+2}{\frac {\partial G_{i}}{\partial x_{i+2}}})+({\frac {\partial F_{i}}{\partial x_{i}}}G_{i}+{\frac {\partial F_{i}}{\partial x_{i+1}}}G_{i+1}+{\frac {\partial F_{i}}{\partial x_{i+2}}}G_{i+2}))$
$=(i,F_{i}(\nabla \cdot \mathbf {G} )-(\nabla \cdot \mathbf {F} )G_{i}-(\mathbf {F} \cdot \nabla )G_{i}+(\mathbf {G} \cdot \nabla )F_{i})$ $=(\nabla \cdot \mathbf {G} )\mathbf {F} -(\nabla \cdot \mathbf {F} )\mathbf {G} -(\mathbf {F} \cdot \nabla )\mathbf {G} +(\mathbf {G} \cdot \nabla )\mathbf {F}$ $=((\nabla \cdot \mathbf {G} )\mathbf {F} +(\mathbf {G} \cdot \nabla )\mathbf {F} )-((\nabla \cdot \mathbf {F} )\mathbf {G} +(\mathbf {F} \cdot \nabla )\mathbf {G} )$
The following identity is a very important property of vector fields which are the gradient of a scalar field. A vector field which is the gradient of a scalar field is always irrotational. Given scalar field $f$, then $\nabla \times (\nabla f)=\mathbf {0}$
$\nabla \times (\nabla f)=\nabla \times (i,{\frac {\partial f}{\partial x_{i}}})$ $=(i,{\frac {\partial }{\partial x_{i+1}}}({\frac {\partial f}{\partial x_{i+2}}})-{\frac {\partial }{\partial x_{i+2}}}({\frac {\partial f}{\partial x_{i+1}}}))$ $=(i,{\frac {\partial ^{2}f}{\partial x_{i+1}\partial x_{i+2}}}-{\frac {\partial ^{2}f}{\partial x_{i+2}\partial x_{i+1}}})$ $=(i,0)$ $=\mathbf {0}$
The following identity is a complex, yet popular identity used for deriving the Helmholtz decomposition theorem. Given vector field $\mathbf {F}$, then $\nabla \times (\nabla \times \mathbf {F} )=\nabla (\nabla \cdot \mathbf {F} )-\nabla ^{2}\mathbf {F}$
$\nabla \times (\nabla \times \mathbf {F} )=\nabla \times (i,{\frac {\partial F_{i+2}}{\partial x_{i+1}}}-{\frac {\partial F_{i+1}}{\partial x_{i+2}}})$ $=(i,{\frac {\partial }{\partial x_{i+1}}}({\frac {\partial F_{i+1}}{\partial x_{i}}}-{\frac {\partial F_{i}}{\partial x_{i+1}}})-{\frac {\partial }{\partial x_{i+2}}}({\frac {\partial F_{i}}{\partial x_{i+2}}}-{\frac {\partial F_{i+2}}{\partial x_{i}}}))$ $=(i,({\frac {\partial ^{2}F_{i+1}}{\partial x_{i}\partial x_{i+1}}}+{\frac {\partial ^{2}F_{i+2}}{\partial x_{i}\partial x_{i+2}}})-({\frac {\partial ^{2}F_{i}}{\partial x_{i+1}^{2}}}+{\frac {\partial ^{2}F_{i}}{\partial x_{i+2}^{2}}}))$ $=(i,({\frac {\partial ^{2}F_{i}}{\partial x_{i}\partial x_{i}}}+{\frac {\partial ^{2}F_{i+1}}{\partial x_{i}\partial x_{i+1}}}+{\frac {\partial ^{2}F_{i+2}}{\partial x_{i}\partial x_{i+2}}})-({\frac {\partial ^{2}F_{i}}{\partial x_{i}^{2}}}+{\frac {\partial ^{2}F_{i}}{\partial x_{i+1}^{2}}}+{\frac {\partial ^{2}F_{i}}{\partial x_{i+2}^{2}}}))$ $=(i,{\frac {\partial }{\partial x_{i}}}({\frac {\partial F_{i}}{\partial x_{i}}}+{\frac {\partial F_{i+1}}{\partial x_{i+1}}}+{\frac {\partial F_{i+2}}{\partial x_{i+2}}})-\nabla ^{2}F_{i})$ $=(i,{\frac {\partial }{\partial x_{i}}}(\nabla \cdot \mathbf {F} )-\nabla ^{2}F_{i})$ $=\nabla (\nabla \cdot \mathbf {F} )-\nabla ^{2}\mathbf {F}$

## Basis Vector Identities

The Cartesian basis vectors $\mathbf {i}$, $\mathbf {j}$, and $\mathbf {k}$ are the same at all points in space. However, in other coordinate systems like cylindrical coordinates or spherical coordinates, the basis vectors can change with respect to position.
In cylindrical coordinates, the unit-length mutually perpendicular basis vectors are ${\hat {\mathbf {\rho } }}=(\cos \phi )\mathbf {i} +(\sin \phi )\mathbf {j}$, ${\hat {\mathbf {\phi } }}=(-\sin \phi )\mathbf {i} +(\cos \phi )\mathbf {j}$, and ${\hat {\mathbf {z} }}=\mathbf {k}$ at position $(\rho ,\phi ,z)$ which corresponds to Cartesian coordinates $(\rho \cos \phi ,\rho \sin \phi ,z)$.
In spherical coordinates, the unit-length mutually perpendicular basis vectors are ${\hat {\mathbf {r} }}=(\sin \theta \cos \phi )\mathbf {i} +(\sin \theta \sin \phi )\mathbf {j} +(\cos \theta )\mathbf {k}$, ${\hat {\mathbf {\theta } }}=(\cos \theta \cos \phi )\mathbf {i} +(\cos \theta \sin \phi )\mathbf {j} +(-\sin \theta )\mathbf {k}$, and ${\hat {\mathbf {\phi } }}=(-\sin \phi )\mathbf {i} +(\cos \phi )\mathbf {j}$ at position $(r,\theta ,\phi )$ which corresponds to Cartesian coordinates $(r\sin \theta \cos \phi ,r\sin \theta \sin \phi ,r\cos \theta )$.
It should be noted that ${\hat {\mathbf {\phi } }}$ is the same in both cylindrical and spherical coordinates.
This section will compute the directional derivative and Laplacian for the following vectors since these quantities do not immediately follow from the formulas established for the directional derivative and Laplacian for scalar fields in various coordinate systems.
    ${\hat {\mathbf {\rho } }}$ which is the unit length vector that points away from the z-axis and is perpendicular to the z-axis.
    ${\hat {\mathbf {\phi } }}$ which is the unit length vector that points around the z-axis in a counterclockwise direction and is both parallel to the xy-plane and perpendicular to the position vector projected onto the xy-plane.
    ${\hat {\mathbf {r} }}$ which is the unit length vector that points away from the origin.
    ${\hat {\mathbf {\theta } }}$ which is the unit length vector that is perpendicular to the position vector and points "south" on the surface of a sphere that is centered on the origin.
The following quantities are also important:
    $\rho$ which is the perpendicular distance from the z-axis.
    $\phi$ which is the azimuth: the counterclockwise angle of the position vector relative to the x-axis after being projected onto the xy-plane.
    $r$ which is the distance from the origin.
    $\theta$ which is the angle of the position vector to the z-axis.

### Vector Rho

${\hat {\mathbf {\rho } }}$ only changes with respect to $\phi$: ${\frac {\partial {\hat {\mathbf {\rho } }}}{\partial \phi }}={\hat {\mathbf {\phi } }}$.
Given vector field $\mathbf {V} =\mathbf {V} _{\perp }+v_{\phi }{\hat {\mathbf {\phi } }}$ where $\mathbf {V} _{\perp }$ is always orthogonal to ${\hat {\mathbf {\phi } }}$, then $(\mathbf {V} \cdot \nabla ){\hat {\mathbf {\rho } }}={\frac {v_{\phi }}{\rho }}{\hat {\mathbf {\phi } }}$
Using cylindrical coordinates, let $\mathbf {V} _{\perp }=v_{\rho }{\hat {\mathbf {\rho } }}+v_{z}{\hat {\mathbf {z} }}$
The cylindrical coordinate version of the directional derivative gives:
$(\mathbf {V} \cdot \nabla ){\hat {\mathbf {\rho } }}=((v_{\rho }{\hat {\mathbf {\rho } }}+v_{\phi }{\hat {\mathbf {\phi } }}+v_{z}{\hat {\mathbf {z} }})\cdot \nabla ){\hat {\mathbf {\rho } }}$ $=v_{\rho }{\frac {\partial {\hat {\mathbf {\rho } }}}{\partial \rho }}+{\frac {v_{\phi }}{\rho }}{\frac {\partial {\hat {\mathbf {\rho } }}}{\partial \phi }}+v_{z}{\frac {\partial {\hat {\mathbf {\rho } }}}{\partial z}}$ $=v_{\rho }\mathbf {0} +{\frac {v_{\phi }}{\rho }}{\hat {\mathbf {\phi } }}+v_{z}\mathbf {0}$ $={\frac {v_{\phi }}{\rho }}{\hat {\mathbf {\phi } }}$
$\nabla ^{2}{\hat {\mathbf {\rho } }}=-{\frac {1}{\rho ^{2}}}{\hat {\mathbf {\rho } }}$
Using the cylindrical coordinate version of the Laplacian,
$\nabla ^{2}{\hat {\mathbf {\rho } }}={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}(\rho {\frac {\partial {\hat {\mathbf {\rho } }}}{\partial \rho }})+{\frac {1}{\rho ^{2}}}{\frac {\partial ^{2}{\hat {\mathbf {\rho } }}}{\partial \phi ^{2}}}+{\frac {\partial ^{2}{\hat {\mathbf {\rho } }}}{\partial z^{2}}}$ $={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}(\rho \mathbf {0} )+{\frac {1}{\rho ^{2}}}{\frac {\partial {\hat {\mathbf {\phi } }}}{\partial \phi }}+{\frac {\partial \mathbf {0} }{\partial z}}$ $=-{\frac {1}{\rho ^{2}}}{\hat {\mathbf {\rho } }}$

### Vector Phi

${\hat {\mathbf {\phi } }}$ only changes with respect to $\phi$: ${\frac {\partial {\hat {\mathbf {\phi } }}}{\partial \phi }}=-{\hat {\mathbf {\rho } }}$.
Given vector field $\mathbf {V} =\mathbf {V} _{\perp }+v_{\phi }{\hat {\mathbf {\phi } }}$ where $\mathbf {V} _{\perp }$ is always orthogonal to ${\hat {\mathbf {\phi } }}$, then $(\mathbf {V} \cdot \nabla ){\hat {\mathbf {\phi } }}=-{\frac {v_{\phi }}{\rho }}{\hat {\mathbf {\rho } }}$
Using cylindrical coordinates, let $\mathbf {V} _{\perp }=v_{\rho }{\hat {\mathbf {\rho } }}+v_{z}{\hat {\mathbf {z} }}$
The cylindrical coordinate version of the directional derivative gives:
$(\mathbf {V} \cdot \nabla ){\hat {\mathbf {\phi } }}=((v_{\rho }{\hat {\mathbf {\rho } }}+v_{\phi }{\hat {\mathbf {\phi } }}+v_{z}{\hat {\mathbf {z} }})\cdot \nabla ){\hat {\mathbf {\phi } }}$ $=v_{\rho }{\frac {\partial {\hat {\mathbf {\rho } }}}{\partial \phi }}+{\frac {v_{\phi }}{\rho }}{\frac {\partial {\hat {\mathbf {\phi } }}}{\partial \phi }}+v_{z}{\frac {\partial {\hat {\mathbf {\phi } }}}{\partial z}}$ $=v_{\rho }\mathbf {0} +{\frac {v_{\phi }}{\rho }}(-{\hat {\mathbf {\rho } }})+v_{z}\mathbf {0}$ $=-{\frac {v_{\phi }}{\rho }}{\hat {\mathbf {\rho } }}$
$\nabla ^{2}{\hat {\mathbf {\phi } }}=-{\frac {1}{\rho ^{2}}}{\hat {\mathbf {\phi } }}$
Using the cylindrical coordinate version of the Laplacian,
$\nabla ^{2}{\hat {\mathbf {\phi } }}={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}(\rho {\frac {\partial {\hat {\mathbf {\phi } }}}{\partial \rho }})+{\frac {1}{\rho ^{2}}}{\frac {\partial ^{2}{\hat {\mathbf {\phi } }}}{\partial \phi ^{2}}}+{\frac {\partial ^{2}{\hat {\mathbf {\phi } }}}{\partial z^{2}}}$ $={\frac {1}{\rho }}{\frac {\partial }{\partial \rho }}(\rho \mathbf {0} )-{\frac {1}{\rho ^{2}}}{\frac {\partial {\hat {\mathbf {\rho } }}}{\partial \phi }}+{\frac {\partial \mathbf {0} }{\partial z}}$ $=-{\frac {1}{\rho ^{2}}}{\hat {\mathbf {\phi } }}$

### Vector r

${\hat {\mathbf {r} }}$ changes with respect to $\theta$ and $\phi$: ${\frac {\partial {\hat {\mathbf {r} }}}{\partial \theta }}={\hat {\mathbf {\theta } }}$ and ${\frac {\partial {\hat {\mathbf {r} }}}{\partial \phi }}=(\sin \theta ){\hat {\mathbf {\phi } }}$
Given vector field $\mathbf {V} =v_{r}{\hat {\mathbf {r} }}+v_{\theta }{\hat {\mathbf {\theta } }}+v_{\phi }{\hat {\mathbf {\phi } }}$, then $(\mathbf {V} \cdot \nabla ){\hat {r}}={\frac {1}{r}}(v_{\theta }{\hat {\mathbf {\theta } }}+v_{\phi }{\hat {\mathbf {\phi } }})$
The spherical coordinate version of the directional derivative gives:
$(\mathbf {V} \cdot \nabla ){\hat {\mathbf {r} }}=((v_{r}{\hat {\mathbf {r} }}+v_{\theta }{\hat {\mathbf {\theta } }}+v_{\phi }{\hat {\mathbf {\phi } }})\cdot \nabla ){\hat {\mathbf {r} }}$ $=v_{r}{\frac {\partial {\hat {\mathbf {r} }}}{\partial r}}+{\frac {v_{\theta }}{r}}{\frac {\partial {\hat {\mathbf {r} }}}{\partial \theta }}+{\frac {v_{\phi }}{r\sin \theta }}{\frac {\partial {\hat {\mathbf {r} }}}{\partial \phi }}$ $=v_{r}\mathbf {0} +{\frac {v_{\theta }}{r}}{\hat {\mathbf {\theta } }}+{\frac {v_{\phi }}{r\sin \theta }}(\sin \theta {\hat {\mathbf {\phi } }})$ $={\frac {1}{r}}(v_{\theta }{\hat {\mathbf {\theta } }}+v_{\phi }{\hat {\mathbf {\phi } }})$
$\nabla ^{2}{\hat {\mathbf {r} }}=-{\frac {2}{r^{2}}}{\hat {\mathbf {r} }}$
The spherical coordinate version of the Laplacian gives:
$\nabla ^{2}{\hat {\mathbf {r} }}={\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}(r^{2}{\frac {\partial {\hat {\mathbf {r} }}}{\partial r}})+{\frac {1}{r^{2}\sin \theta }}{\frac {\partial }{\partial \theta }}(\sin \theta {\frac {\partial {\hat {\mathbf {r} }}}{\partial \theta }})+{\frac {1}{r^{2}\sin ^{2}\theta }}{\frac {\partial ^{2}{\hat {\mathbf {r} }}}{\partial \phi ^{2}}}$ $={\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}(r^{2}\mathbf {0} )+{\frac {1}{r^{2}\sin \theta }}{\frac {\partial }{\partial \theta }}(\sin \theta {\hat {\mathbf {\theta } }})+{\frac {1}{r^{2}\sin ^{2}\theta }}{\frac {\partial }{\partial \phi }}(\sin \theta {\hat {\mathbf {\phi } }})$ $={\frac {1}{r^{2}\sin \theta }}(\cos \theta {\hat {\mathbf {\theta } }}+\sin \theta (-{\hat {\mathbf {r} }}))+{\frac {1}{r^{2}\sin \theta }}(-\sin \theta {\hat {\mathbf {r} }}-\cos \theta {\hat {\mathbf {\theta } }})$ $=-{\frac {2}{r^{2}}}{\hat {\mathbf {r} }}$

### Vector Theta

${\hat {\mathbf {\theta } }}$ changes with respect to $\theta$ and $\phi$: ${\frac {\partial {\hat {\mathbf {\theta } }}}{\partial \theta }}=-{\hat {\mathbf {r} }}$ and ${\frac {\partial {\hat {\mathbf {\theta } }}}{\partial \phi }}=(\cos \theta ){\hat {\mathbf {\phi } }}$
Given vector field $\mathbf {V} =v_{r}{\hat {\mathbf {r} }}+v_{\theta }{\hat {\mathbf {\theta } }}+v_{\phi }{\hat {\mathbf {\phi } }}$, then $(\mathbf {V} \cdot \nabla ){\hat {\theta }}={\frac {1}{r}}(-v_{\theta }{\hat {\mathbf {r} }}+\cot \theta v_{\phi }{\hat {\mathbf {\phi } }})$
The spherical coordinate version of the directional derivative gives:
$(\mathbf {V} \cdot \nabla ){\hat {\mathbf {\theta } }}=((v_{r}{\hat {\mathbf {r} }}+v_{\theta }{\hat {\mathbf {\theta } }}+v_{\phi }{\hat {\mathbf {\phi } }})\cdot \nabla ){\hat {\mathbf {\theta } }}$ $=v_{r}{\frac {\partial {\hat {\mathbf {\theta } }}}{\partial r}}+{\frac {v_{\theta }}{r}}{\frac {\partial {\hat {\mathbf {\theta } }}}{\partial \theta }}+{\frac {v_{\phi }}{r\sin \theta }}{\frac {\partial {\hat {\mathbf {\theta } }}}{\partial \phi }}$ $=v_{r}\mathbf {0} +{\frac {v_{\theta }}{r}}(-{\hat {\mathbf {r} }})+{\frac {v_{\phi }}{r\sin \theta }}(\cos \theta {\hat {\mathbf {\phi } }})$ $={\frac {1}{r}}(-v_{\theta }{\hat {\mathbf {r} }}+\cot \theta v_{\phi }{\hat {\mathbf {\phi } }})$
$\nabla ^{2}{\hat {\mathbf {\theta } }}=-{\frac {1}{r^{2}\sin \theta }}(2\cos \theta {\hat {\mathbf {r} }}+\csc \theta {\hat {\mathbf {\theta } }})=-{\frac {1}{r^{2}\sin ^{2}\theta }}(\sin(2\theta ){\hat {\mathbf {r} }}+{\hat {\mathbf {\theta } }})$
The spherical coordinate version of the Laplacian gives:
$\nabla ^{2}{\hat {\mathbf {\theta } }}={\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}(r^{2}{\frac {\partial {\hat {\mathbf {\theta } }}}{\partial r}})+{\frac {1}{r^{2}\sin \theta }}{\frac {\partial }{\partial \theta }}(\sin \theta {\frac {\partial {\hat {\mathbf {\theta } }}}{\partial \theta }})+{\frac {1}{r^{2}\sin ^{2}\theta }}{\frac {\partial ^{2}{\hat {\mathbf {\theta } }}}{\partial \phi ^{2}}}$ $={\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}(r^{2}\mathbf {0} )+{\frac {1}{r^{2}\sin \theta }}{\frac {\partial }{\partial \theta }}(\sin \theta (-{\hat {\mathbf {r} }}))+{\frac {1}{r^{2}\sin ^{2}\theta }}{\frac {\partial }{\partial \phi }}(\cos \theta {\hat {\mathbf {\phi } }})$ $=-{\frac {1}{r^{2}\sin \theta }}(\cos \theta {\hat {\mathbf {r} }}+\sin \theta {\hat {\mathbf {\theta } }})+{\frac {\cos \theta }{r^{2}\sin ^{2}\theta }}(-\sin \theta {\hat {\mathbf {r} }}-\cos \theta {\hat {\mathbf {\theta } }})$ $=-{\frac {1}{r^{2}\sin \theta }}(2\cos \theta {\hat {\mathbf {r} }}+(\sin \theta +{\frac {\cos ^{2}\theta }{\sin \theta }}){\hat {\mathbf {\theta } }})$ $=-{\frac {1}{r^{2}\sin \theta }}(2\cos \theta {\hat {\mathbf {r} }}+\csc \theta {\hat {\mathbf {\theta } }})$