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