Fiaure 1. Schematic représentation of the molecular model for water used in this work.

Considérons à présent les principales étapes de la méthode de Monte Carlo. La pression de vapeur et les densités de vapeur calculées ont été comparées avec les données expérimentales. Une nouvelle expression a été tirée des calculs de la pression dans les systèmes où a été utilisée la méthode d'Ewald. Afin d'améliorer l'accord de la simulation avec les données expérimentales, une méthode a été présentée pour la réévaluation des paramètres du modèle moléculaire. Pour les modèles examinés, il n'y a pas d'ensemble de paramètres unique qui donne simultanément une bonne description de la pression de vapeur et des densités expérimentales. Ensuite, l'enthalpie de vaporisation, les fonctions de distribution radiales et un coefficient caractéristique appelé second coefficient virial ont été calculés pour chaque modèle et comparés aux données expérimentales.

Il en a été déduit que les résultats du modèle SPC/E constituaient une bonne prédiction pour la température critique et la densité critique. Quant aux valeurs de pression critique, les deux modèles sont assez loin des données expérimentales.

(*)Méthode de sommation d'Ewald = méthode calculatoire appliquée aux énergies.

b) Par une autre méthode

Un autre groupe de chercheurs a réalisé des simulations moléculaires dynamiques en utilisant des modèles simples non polarisables de l'eau et un nouveau modèle polarisable. Il en est ressorti que l'utilisation de potentiels intermoléculaires issus de la littérature apportait de la précision dans les simulations. Afin d'accroître encore cette précision, il est intéressant d'inclure les effets de la polarisabilité dans le potentiel de l'eau.

25