I am simply trying to understand where the factor of $2$ comes from in the RHS of $$\epsilon_{ijk} \left[\hat a_j, \hat b_k \right]=\epsilon_{ijk}\left(\hat a_j \hat b_k-\hat b_k \hat a_j\right)=\epsilon_{ijk}\left(2\hat a_j\hat b_k\right)\tag{a}$$
Here, $\left[\hat a_j, \hat b_k \right]$ is the commutator of operators $\hat a_j,$ and $\hat b_k$ and the Einstein summation convention is being used. Here the operators must not be assumed to commute or anticommute.
I once asked my lecturer about why eqn. $(\mathrm{a})$ is true and the usual story is that the "Levi-Civita is anti-symmetric under exchange of $j$ and $k$, so the two terms are the same up to a minus sign which gives the extra factor of $2$". But I'm finding it hard to understand this statement.
So in trying to convince myself that $(\mathrm{a})$ is true, starting with the LHS of $(\mathrm{a})$, I write it out for $i$th component by summing over $j$ and $k$ which run from $1$ to $2$:
$$\epsilon_{ijk} \left[\hat a_j, \hat b_k \right]=\epsilon_{ijk}\left(\hat a_j \hat b_k-\hat b_k \hat a_j\right)\tag{1}$$ $$=\epsilon_{i12}\left(\hat a_1 \hat b_2-\hat b_2 \hat a_1\right)+\epsilon_{i21}\left(\hat a_2 \hat b_1-\hat b_1 \hat a_2\right)$$ $$=\epsilon_{i12}\left(\hat a_1 \hat b_2-\hat b_2 \hat a_1\right)-\epsilon_{i12}\left(\hat a_2 \hat b_1-\hat b_1 \hat a_2\right)$$ $$=\epsilon_{i12}\left((\hat a_1 \hat b_2-\hat a_2 \hat b_1)-(\hat b_2\hat a_1-\hat b_1 \hat a_2)\right)$$ $$\stackrel{\color{red}{?}}{=}\epsilon_{ijk}\left(\hat a_j \hat b_k-\hat b_k \hat a_j\right)\tag{2}$$ In going from the second to the third equality I used the fact that $\epsilon_{i21}=-\epsilon_{i12}$ due to the fact that the Levi-Civita is totally antisymmetric under interchange of $2$ indices. But the final expression, $(2)$, doesn't show where the factor of $2$ comes from in $(a)$. All the final expression in $(2)$ shows is that the starting equation, $(1)$ is true (assuming my reasoning is correct). However, I'm unsure about the final equality and indicated this with a question mark.
So instead, starting from the RHS of $(\mathrm{a})$ $$\epsilon_{ijk}\left(2\hat a_j\hat b_k\right)=2\left(\epsilon_{i12}\hat a_1\hat b_2+\epsilon_{i21}\hat a_2 \hat b_1 \right)\tag{3}$$ $$=2\left(\epsilon_{i12}\hat a_1\hat b_2-\epsilon_{i12}\hat a_2 \hat b_1 \right)=2\epsilon_{i12}\left(\hat a_1\hat b_2-\hat a_2 \hat b_1 \right)$$ $$=2\epsilon_{i12}\left(\hat a_1\hat b_2-\hat a_2 \hat b_1 \right)=2\epsilon_{ijk}\hat a_j\hat b_k \tag{4}$$
But once again, in eqn. $(4)$ I've proven absolutely nothing, assuming my manipulations are correct I have simply proven that the starting equation, $(3)$ is true.
So to summarize, can anyone please show me where this factor of $2$ comes from? Or put another way, show that $$\fbox{$\epsilon_{ijk}\left(\hat a_j \hat b_k-\hat b_k \hat a_j\right)=2\epsilon_{ijk}\hat a_j\hat b_k$}\tag{a}$$
The proof does not have to be rigorous, it need only apply to a particular case, (like I tried to calculate above for the case of $j$ and $k$ running from $1$ to $2$).
Starting with: $\varepsilon_{ijk}a_jb_ k-\varepsilon_{ijk}a_kb_j $
We take the second term: $-\varepsilon_{ijk}a_kb_j=+\varepsilon_{ikj}a_kb_j$. This works because $ijk\to ikj$ is an odd permutation, so it generates a negative sign. We’re summing over $j$ and $k$, so provided we stay consistent we can rename them. Thus, the second term becomes $+\varepsilon_{ijk}a_jb_k$, which is identical to the first term, hence the factor of 2.
EDIT
Starting with: $\varepsilon_{ijk}[a_j,b_k]=\varepsilon_{ijk}a_jb_k-\varepsilon_{ijk}b_ka_j$ (Compare this to the expression above; they are different)
Then, manipulating this:
$$\varepsilon_{ijk}a_jb_k-\varepsilon_{ijk}b_ka_j=\varepsilon_{ijk}a_jb_k+\varepsilon_{ikj}b_ka_j=\varepsilon_{ijk}a_jb_k+\varepsilon_{ijk}b_ja_k$$
(Note how this is different from the math before the edit. The idea is the same, but the operators are never interchanged)
Continuing forward by adding "$0$":
$$\varepsilon_{ijk}(a_jb_k+b_ja_k+a_jb_k-a_jb_k)=2\varepsilon_{ijk}a_jb_k+\varepsilon_{ijk}(b_ja_k-a_jb_k)$$
$$=2\varepsilon_{ijk}a_jb_k-\varepsilon_{ikj}b_ja_k-\varepsilon_{ijk}a_jb_k=2\varepsilon_{ijk}a_jb_k-\varepsilon_{ijk}b_ka_j-\varepsilon_{ijk}a_jb_k=2\epsilon_{ijk}a_jb_k-\varepsilon_{ijk}\{a_j,b_k\}$$
Where $\{a_k,b_j\}$ denotes an anticommutator. But this means:
$$\varepsilon_{ijk}[a_j,b_k]=2\varepsilon_{ijk}a_jb_k-\varepsilon_{ijk}\{a_j,b_k\}$$
Thus, the property you give in your question is only possible if $\{a_j,b_k\}=0$.