I have a Fermionic Hubbard Hamiltonian in the second quantized form. How can I encode it in qiskit, so that I can use `VQE`

to find the ground state? The Hamiltonian is written in the following form:

$$

H = -t \sum_{i} (a_i^\dagger a_{i+1} + \frac{1}{2} U a_i^\dagger a_i) + h.c.

$$

where $a^\dagger$ is a fermionic creation operator, and $i$ is the location index on a 1-D chain from 0 to 3.