The free fermion Hamiltonian for the 1D quantum Ising model is $$H = -J\sum_i (c_{i}^{\dagger }c_{i+1} +c_{i+1}^{\dagger }c_{i}+c_{i}^{\dagger }c_{i+1}^{\dagger }+c_{i+1}c_{i}-2gc_{i}^{\dagger }c_{i} +g)$$ where the sum is over $i$, the site index.
My question is - what are the basis states of this system? I read somewhere that the matrix representation of the Hamiltonian of this form has $L \times L$ elements where $L$ are the total sites. I don't see how that is. (As opposed to the spin states where the matrix would be $2^{L}\times 2^{L}$)
I eventually want to find the ground state but do not want to resort to a Fourier transform to the momentum states, so I want the matrix in the fermionic basis states form only.