My professor showed me this in her office today but I didn't like her method and wanted to use another method.
So, I computed the characteristic polynomial of some arbitrary $2 \times 2$ matrix with nonnegative entries $a,b,c,d$, and solved for its two eigenvalues using the quadratic formula.
Then I equated the conjugate pair of complex numbers (with non-zero imaginary part) to the two eigenvalues, and am trying to derive some contradiction at this point. But, I'm stuck.
My professor had used some complexification process and showed that the realizing matrix must have certain trigonometric entries (she mentioned rotation matrices), and then derived an impossibility.
Thanks in advance,
I'll assume your matrix is
$$ M=\left( \begin{array}{ccc} a & b\\ c & d\\\end{array} \right). $$
The the characteristic polynomial is
$$ |M-\lambda I|=(a-x)(d-x)-bc $$
and setting the characteristic polynomial to zero yields the solution
$$ \lambda = \frac 12 \left(a + d \pm\sqrt{4bc + (a-d)^2}\right). $$
The eigenvalues are only imaginary if
$$ 4bc + (a-d)^2<0 $$
but by assumption $b>0$ and $c>0$. The second term is positive as long as $a$ and $d$ are real (or their complex parts cancel). Thus, the eigenvalues must be real.