Generate sparse matrix representations of the fermionic creation and annihilation operators for nmodes modes (or sites), based on Jordan-Wigner transformation.