Les électrons dans les atomes lourds, en particulier ceux qui sont proches du noyau, sont soumis à des effets relativistes importants. Il est nécessaire de prendre en compte ces effets si l'on veut, par exemple, décrire précisément les niveaux d'énergies des atomes. L'étude des modèles atomiques quantiques relativistes remonte aux travaux fondateurs de P.A.M. Dirac, dès 1928. Ses travaux ont permis d'anticiper la découverte des antiparticules. En effet, le hamiltonien quantique qu'il obtient pour l'atome d'hydrogène n'a de sens physique que si l'on peut interpréter ses énergies négatives comme celles d'une mer infinie de particules virtuelles. Un « trou » dans le spectre des énergies négatives est alors interprété comme l'apparition d'une anti-particule : le positron. Peu après, en 1938, pour étudier les atomes à plusieurs électrons Swirles propose un modèle d'approximation qui donnera lieu aux fameuses équations de Dirac-Fock. Cette approche qui est auto-consistante, et pour laquelle les équations obtenues sont non linéaires, permet une étude numérique dont les résultats sont en très bon accord avec les mesures expérimentales. Pour autant, la motivation physique de cette approche reste incomplète. Elle s'appuie essentiellement sur l'analogue non relativiste des modèles atomiques quantiques, mais ne tient pas compte de l'interprétation de Dirac. De plus, le lien des équations de Dirac-Fock avec l'approche théorique donnée par l'électrodynamique quantique (QED) reste à établir clairement. En particulier, en QED, la question de la définition d'un espace qui décrit les états électroniques reste posée. Le travail présenté ici est une tentative d'apporter quelques réponses mathématiques rigoureuses sur ces problèmes. Nous commencerons par construire une famille de fonctionnelles à partir du hamiltonien formel de la QED qui dépendra du choix de l'espace à un électron. On se placera dans l'approximation de Hartree-Fock. On étudiera alors le problème de la stabilité, celui de l'existence de minima pour ces fonctionnelles (avec ou sans condition de charge totale fixée). On se consacrera ensuite à l'exposé des résultats obtenus qui permettent de comparer les deux approches : « Equations de Dirac-Fock » et « QED dans l'approximation de Hartree-Fock ». On distinguera en particulier le cas des couches pleines qui conduit aux mêmes résultats dans les deux cas, tout au moins pour des constantes de couplages faibles.