Ce document est le fruit d'un long travail approuvé par le jury de ...

3.9. Cas d’une simulation de l’eau 97
moléculaire pour des molécules d’eau. Pour cela, j’ai réalise une simulation du modèle
polarisable de Dang et Chang [337, 339] au cours de laquelle j’ai prélevé 861 dimères à
des distances intermoléculaires inférieures à 10 Å. Ce dernier critère est pris car le rayon
de coupure (cut-off ) choisi pour la simulation est de 11 Å.
Pour chaque dimère de la simulation, un calcul quantique (corrigé de la BSSE, Basis
Set Superposition Error) permet d’obtenir l’énergie d’interaction. On retranche alors à
l’énergie totale la contribution de l’énergie d’interaction électrostatique et d’induction
de notre modèle empirique. Cela permet d’obtenir un nuage de points pour différentes
conformations. Un programme 9 nous permet, en utilisant un algorithme de Levenberg-
Marquardt, de trouver de nouveaux paramètres pour la molécule d’eau. Les nouveaux
paramètres sont :
– ε = 2845438.0
– A1ii = 4.7532
– A2ii = 2.6936
Une simulation de dynamique moléculaire de l’eau réalisée avec ce nouveau modèle
conduit à une densité de 0.98 g.cm 3 et une distribution radiale (voir figure 3.19) très proche
de l’expérience, excepté une intensité du premier pic trop forte. Ces premiers résultats,
très encourageants, du modèle pour les simulations de dynamique moléculaire montrent
qu’une approche toute ab initio est plausible pour obtenir des paramètres non biaisés par
une reproduction de quelques valeurs obtenues expérimentalement.
9. Ce programme, fitmultisite.py, a été développé en python spécifiquement au laboratoire pour
ajuster les paramètres de van der Waals. Celui-ci permettra aussi de déterminer des paramètres pour
divers types d’atomes pour la suite du projet.