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^{†}{a}_{i+1}+\frac{1}{2}U{a}_i^{†}{a}_i)+h.c.$$
where $a^{†}$ is a fermionic creation operator, and $i$ is the location index on a 1-D chain from 0 to 3.