$\underline{\textbf{Model:}}$
Let we have the $t-V$ model for spinless fermions on a 1D lattice, which is defined in second quantization operators as follows:
$$H_1 = -t\sum_i \big(c_i^\dagger c_{i+1} + H.C.\big) + V\sum_i n_i n_{i+1} \qquad (eq.1)$$ here $c_i^\dagger (c_i)$ are creation (annahilation) operators and $n_i$ is number operator.
$\underline{\textbf{JW transformation:}}$
$$c_i^\dagger \rightarrow S_i^+$$ $$c_i \rightarrow S_i^-$$ $$n_i \rightarrow S_i^z + \frac{1}{2}$$ where $S_i$ are spin-1/2 operators. (eq.1) become
$$H_2 = -t\sum_i \big(S_i^+ S_{i+1}^- + H.C.\big) + V\sum_i \big(S_i^z S_{i+1}^z + S_i^z\big) +\frac{NV}{4}\qquad (eq.2)$$ where $N$ is total number of spins.
Ground state energy of $(eq.1)$ can be found at any $N/L$-filling ($N=$ total number of particles, $L=$ total number of sites) using exact diagonalization (as explained here).
$\underline{\textbf{Question:}}$
$(eq.2)$ can also be numerically diagonalized (for small system size). But how can I find its ground state energy at some specific filling, let's say half-filling?