El operador lineal $B \to AB + BA$ corresponde a la $4 \times 4$ matriz
$$ M = \pmatrix{2a & b & c & 0\cr c & a+d & 0 & c\cr b & 0 & a+d & b\cr 0 & b & c & 2d\cr}$$
el uso de la base $\left(\pmatrix{1 & 0\cr 0 & 0\cr}, \pmatrix{0 & 1\cr 0 & 0\cr},\pmatrix{0 & 0\cr 1 & 0\cr}, \pmatrix{0 & 0\cr 0 & 1\cr}\right)$ de la $2 \times 2$ matrices.
Esto ha determinante $$-4\,bc{d}^{2}-8\,bcad-4\,ac{a}^{2}+4\, {d}^{3}+8\,{d}^{2}{a}^{2}+4\,d{
un}^{3}=4\, \a la izquierda( a+d \right) ^{2} \left( ad-bc \right) $$
así que hay distinto de cero soluciones para $B$ si y sólo si $a+d=0$ o $ad-bc=0$.
Si $a+d=0$, sustituyendo $d = -a$ $M$ tenemos
$$ \pmatrix{2a & b & c & 0\cr c & 0 & 0 & c\cr b & 0 & 0 & b\cr 0 & b & c & -2a\cr}$$
Mientras $a,b,c$ no son todos los $0$, esto tiene rango $2$, por lo que hay un $2$-dimensional del espacio lineal de $B$ que $AB + BA = 0$ en este caso. Por otro lado, si $ad - bc = 0$ pero $a+d \ne 0$ la matriz debe tener rango de $3$ y solo debemos tener un $1$-dimensional del espacio lineal de $B$: de hecho, la $(1,1)$, $(2,3)$, $(3,2)$ y $(4,4)$ entradas de la clásica adjunto de la matriz son $2(d^2 + ad-bc)(a+d)$, ($2 c^2 (a+d)$, $2 b^2 (a+d)$ y $2 (a^2 + ad-bc)(a+d)$
respectivamente: si $ad-bc=0$, $a+d \ne 0 $ y el rango es menor que $3$, todos estos serían $0$, lo que implica $a=b=c=d=0$.
EDIT: En el $3 \times 3$ de los casos, el factor determinante de la $9 \times 9$ matriz es
$$\eqalign{8\, y\left( a_{{3,3}}a_{{1,1}}a_{{2,2}}-a_{{2,3}}a_{{3,2}}a_{{1,1}}+a_{
{1,3}}a_{{2,1}}a_{{3,2}}-a_{{3,3}}a_{{1,2}}a_{{2,1}}+a_{{2,3}}a_{{1,2}
}a_{{3,1}}-a_{{1,3}}a_{{3,1}}a_{{2,2}} \right)\cr &( -a_{{2,3}}a_{{1
,2}}a_{{3,1}}-a_{{1,3}}a_{{3,1}}a_{{3,3}}-a_{{1,2}}a_{{2,1}}a_{{2,2}}+
{a_{{3,3}}}^{2}a_{{1,1}}+a_{{1,1}}{a_{{2,2}}}^{2}-a_{{1,3}}a_{{1,1}}a_
{{3,1}}-a_{{1,1}}a_{{1,2}}a_{{2,1}}+{a_{{3,3}}}^{2}a_{{2,2}}\cr&+2\,a_{{3,
3}}a_{{1,1}}a_{{2,2}}+a_{{3,3}}{a_{{2,2}}}^{2}-a_{{2,3}}a_{{3,2}}a_{{3
,3}}-a_{{1,3}}a_{{2,1}}a_{{3,2}}-a_{{2,3}}a_{{3,2}}a_{{2,2}}+{a_{{1,1}
}}^{2}a_{{2,2}}+a_{{3,3}}{a_{{1,1}}}^{2} ) ^{2}\cr}
$$
A mí me parece que hay un $1$-dimensional del espacio de solución, cuando uno de estos factores es $0$ $2$- dimensiones del espacio cuando el otro es $0$. A continuación, hay un buen número de casos en los que ambos factores son $0$, algunos de los cuales conducen a $3$-dimensiones de los espacios.