 Solution

For each complex fermion define two Majorana fermion operatorsThey obey $c_r^{†}=c_r$, $\{c_r,c_s\}=2{\delta }_{rs}$, and $a_j=(c_{2j-1}+i{c}_{2j})/2$. Substituting these inverse relations expands every hopping and pairing monomial as a bilinear in the $c_r$. Diagonal terms $c_r^2=1$ contribute only a constant.

For $r\ne s$, $c_r{c}_s$ is anti-Hermitian. Hermiticity of the original quadratic fermion Hamiltonian therefore makes its coefficient purely imaginary, so it can be written $i{A}_{rs}$ with $A_{rs}$ real. Sincethe symmetric part again changes only the constant. Absorbing conventional factors into $A$ gives