In OpenFermion you can create a Hamiltonian in terms of creation and annihilation pretty easily:
ham = ((FermionOperator('1^ 1', .5)) + (FermionOperator('2^ 2', .25)))
And getting the eigenvalues of the Hamiltonian is pretty straightforward as well:
vals = eigenspectrum(ham)
But I don't understand how to get the eigenvectors of the Hamiltonian. FermionOperator doesn't return a numpy object so I can't use usual linear algebra libraries. What am I missing?