We study a mean-field relativistic model which is able to describe both the behavior of finitely many spin-1/2 particles like electrons and of the Dirac sea which is self-consistently polarized in the presence of the real particles. The model is derived from the QED Hamiltonian in Coulomb gauge neglecting the photon field. All our results are non-perturbative and mathematically rigorous.
