For real $a,b,c$ and imaginary $d$ the matrix $A$ has chiral symmetry, meaning it anticommutes with thea matrix $$X=\left( \begin{array}{cccc} 0 & 0 & 0 & -i \\ 0 & 0 & i & 0 \\ 0 & -i & 0 & 0 \\ i & 0 & 0 & 0 \\ \end{array} \right),\;\;XA+AX=0,$$ that $X$ that squares to unitythe identity: $X^2=I$. Hence $$X=\left( \begin{array}{cccc} 0 & 0 & 0 & -i \\ 0 & 0 & i & 0 \\ 0 & -i & 0 & 0 \\ i & 0 & 0 & 0 \\ \end{array} \right),\;\;XA+AX=0,\;\;X^2=I.$$ Hence the spectrum of $A$ has $\pm$ symmetry: $$\det (\lambda-A)=\det(\lambda X^2-XAX)=\det(\lambda+X^2A)=\det(\lambda+A),$$ so if $\lambda$ is an eigenvalue then also $-\lambda$.