Rôle des glaces interstellaires dans la complexité moléculaire de l’espace : modélisation par les méthodes de la chimie théorique

It is postulated nowadays that complex organic molecules in space form on the surface and in the volume of interstellar ices. These ices can catalyse chemical reactions what could explain the formation of prebiotic molecules. However, because of the low temperatures, the diffusion of the reactants one towards another is slow, limiting their reactivity. The objectif of this thesis is to treat the reactivity and the diffusion of simple molecules in interstellar ices. I present in this thesis the results of the study of the chemical reaction between carbon dioxide and ammonia in interstellar ices. This study was conducted in the Laboratoire Univers et Particules de Montpellier : it combines various theoretical chemistry methods and confront the results to experimental ones I participated in at the Laboratoire de Physique des Interactions Ioniques etMoléculaires of the Aix-Marseille university. In a first part, classical molecular dynamic simulations are used to simulate a low-density amorphous ice model and to calculate the diffusion coefficients at various temperatures of a series of small molecules (NH3, CO, CO2, H2CO). These results are compared to the experimental diffusion coeficients of the carbon dioxide validating the theoretical approach used and allowing to porpose a mechanism for the diffusion process. In a second part, the reaction between carbon dioxide and ammonia is studied in the frame of Density-Functional Theory using a « super-molecular » approach. In this approach, the reaction energy profile for the molecular complexes xNH3:CO2:yH2O, is studied. Two reaction products are localised : the ammonium carbamate and the carbamic acid. The reaction energy barrier obtained by these calculations is similar to the one obtained experimentally, and the ammonium carbamate is confirmed as the major reaction product. The obtained energy profile is therfore investigated by constrained ab initio molecular dynamics and the free energy profile is computed with the Thermodynamics Integration method. These calculations confirme the general form of the previous energy profile and enlight the strong entropic effect of the water network. Finally, conclusions are drawn, and perspectives on methodological as well as on astrochemical aspects, as the inclusion of the reaction and diffusion energy barriers I calculated in astrochemical models, are given.Les glaces du milieu interstellaire sont invoquées comme l'une des origines de la formation de molécules organiques complexes dans l'espace. En effet, elles constituent un support catalytique pour des réactions chimiques et pourraient ainsi expliquer la formation de molécules d'intérêt prébiotique. Toutefois, en raison de la faible température des milieux considérés, la vitesse de la réaction est contrainte par le déplacement des différentes réactifs l'un vers l'autre. L'objectif de cette thèse est donc de traiter la réactivité et la diffusion de molécules simples dans les glaces interstellaires. Je présente dans cette thèse les résultats de l'étude de la réaction entre le dioxyde de carbone et l'ammoniac dans les glaces interstellaires. Cette étude a été effectuée au Laboratoire Univers et Particules de Montpellier ; elle combine différentes méthodes de la chimie théorique et confronte les résultats avec ceux issus d'expériences que j'ai réalisées au Laboratoire de Physique des interactions Ioniques et Moléculaires de l'Université d'Aix-Marseille. Dans une première partie, des simulations de dynamique moléculaire classiques sont employées pour simuler un modèle de glace amorphe basse densité. Ce modèle est utilisé pour simuler la trajectoire de petites molécules (NH3, CO, CO2, H2CO) dans cette glace et en déduire des coefficients de diffusion à différentes températures. Ces résultats sont comparés à des résultats expérimentaux de diffusion du dioxyde de carbone ce qui valide la méthode théorique utilisée et permet de suggérer un mécanisme pour ce processus de diffusion. Dans une deuxième partie, la réaction entre le dioxyde de carbone et l'ammoniac est traitée dans le cadre de la théorie de la fonctionnelle densité par une approche « super-molécule ». Dans cette approche, le profil d'énergie et le mécanisme de la réaction dans des complexes moléculaires xNH3:CO2:yH2O sont étudiés. Deux produits de la réaction sont localisés : le carbamate d'ammonium et l'acide carbamique. La barrière d'énergie de la réaction obtenue est similaire à celle obtenue expérimentalement, et le carbamate d'ammonium est confirmé comme produit majoritaire de la réaction. Le profil d'énergie obtenu par cette approche « super-molécule » est ensuite étudié par dynamique moléculaire ab initio contrainte et le profil d'énergie libre est calculé par la méthode d'Intégration Thermodynamique. Cette approche confirme la forme générale du profil d'énergie et met en évidence un fort effet entropique du réseau d'eau. Je dresse finalement des conclusions sur les plans méthodologiques et astrochimiques permettant d'inscrire ma thèse dans des perspectives futures, notamment en incluant les barrières d'énergie de diffusion et de réaction, dans des modèles astrochimiques prenant en compte directement les réactions chimiques dans les glaces interstellaires