How to encode a second-quantized Hamiltonian in qiskit?

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.