The contributions of this thesis consist of three main results. The first result is concerned with analytic perturbation theory for Kohn-Sham type models. We prove, under some technical conditions, the existence, uniqueness and analyticity of the perturbed reduced Hartree-Fock ground state density matrix for regular perturbations arising from an external potential. Our analysis encompasses the case when the Fermi level of the unperturbed ground state is a degenerate eigenvalue of the mean-field operator and the frontier orbitals are partially occupied. The second result is concerned with the mathematical construction of pseudo potentials for Kohn-Sham models. We define a set of admissible semi local norm-conserving pseudo potentials of given local Sobolev regularity and prove that this set is non-empty and closed for an appropriate topology. This allows us to propose a new way to construct pseudo potentials, which consists in optimizing on the latter set some criterion taking into account both smoothness and transferability requirements. The third result is a numerical study of the reduced Hartree-Fock model of atoms. We propose a discretization method and an algorithm to solve numerically the Kohn-Sham equations for an atom subjected to a cylindrically-symmetric external potential. We report the computed occupied energy levels and the occupation numbers for all the atoms of the four first rows of the periodic table and consider the case of an atom subjected to a uniform electric-field