I'm trying to use index notation (i.e. Einstein summation notation) in order to show that $(\vec a \times \vec b) \cdot (\vec a \times \vec b)=|\vec a|^2|\vec b|^2-(\vec a \cdot \vec b)^2$.
Here's what I've done so far, but I'm stuck at a dead end: $$\begin{align}(\vec a \times \vec b) \cdot (\vec a \times \vec b)&=(\varepsilon_{ijk}a_jb_k) (\varepsilon_{ijk}a_jb_k) \\ \ \\&=\varepsilon_{ijk}\varepsilon_{ijk}a_ja_jb_kb_k \\ \ \\ &=[\delta_{jj}\delta_{kk}-\delta_{jk}\delta_{kj}] a_ja_jb_kb_k \\ \ \\ &=[3\cdot 3-\underbrace{\delta_{kk}}_{3} \ ]a_ja_jb_kb_k \\ \ \\&=6a_ja_jb_kb_k\end{align}$$ but where do I go from here? Thanks
To clarify my remark in comments: Using any summation index more than twice renders the summation convention is insensible. This becomes obvious if you write in terms of $\Sigma's$: $$(\vec{a}\times \vec{b})^2=\left(\sum_{j,k=1}^3\epsilon_{ijk}a_j b_k\right)^2\neq \sum_{j,k=1}^3(\epsilon_{ijk}a_j b_k)^2.$$ What is correct is $$(\vec{a}\times \vec{b})^2=\left(\sum_{j,k=1}^3\epsilon_{ijk}a_j b_k\right)\left(\sum_{l,m=1}^3\epsilon_{ilm}a_l b_m\right)= \sum_{jklm}^3 \epsilon_{ijk}a_j b_k\epsilon_{ilm}a_l b_m.$$
So a second pair of indices is absolutely crucial in carrying out the algebra correctly.