This is a branched question of my previous question What's the condition for a matrix $A$ ($2N\times 2N$ dimension) to have eigenvalues in pairs $\pm\lambda$?.
Let's look at this $2N\times2N$ matrix $A$:
$A=\left(\begin{array}{cc} B & C\\ -C^{*} & -B^{*} \end{array}\right)$,
where $B$ is a hermitian matrix, $C$ is a complex symmetric matrix, and $M^\ast$ denotes the complex conjugate (but not the conjugate transpose) of a matrix $M$.
If $A$ is a $2\times2$ matrix, the eigenvalues are obviously in pairs, and I also proved the $4\times4$ case by showing the determinant of $A-\lambda I$ only depends on $\lambda^{2}$.
Then I use Mathematica to calculate the eigenvalues of this matrix with randomly generated matrices $B$ and $C$ for higher dimension cases, the eigenvalues also come in pairs. So I guess this matrix always has eigenvalues in pairs, but I don't know how to prove it analytically.