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

96
Chapitre 3. Développement d’un potentiel intermoléculaire polarisable
simple et précis
trimère tétramère pentamère hexamère
KM Modèle KM Modèle KM Modèle KM Modèle
Électrostatique -27.51 -19.24 -51.74 -33.54 -67.41 -43.03 -85.27 -53.56
Induction -4.83 -4.36 -11.41 -11.62 -16.20 -18.01 -21.62 -25.39
van der Waals 21.03 15.23 42.49 26.73 55.81 34.68 72.36 45.62
Totale -11.31 -8.37 -20.66 -18.42 -27.80 -26.36 -34.53 -33.33
Table 3.13 – Multimères d’eau obtenus à la référence [338]. L’énergie est donnée pour le
modèle et la référence quantique (KM avec le programme Gamess) au niveau de théorie
HF/Sadlej en kcal/mol.
en compte explicitement dans notre modèle. Par contre, nous avons défini un potentiel
de van der Waals à même de compenser cette approximation. Or, plus le nombre de
multimères augmente, moins la compensation est importante. Cela veut donc dire que
le modèle de van der Waals ne reproduit pas fidèlement les interactions selon plusieurs
directions comme c’est le cas avec les multimères.
3.9.1 Nouveaux paramètres de van der Waals
Après avoir fait ces premières analyses sur des multimères, nous avons voulu les utiliser
pour obtenir des paramètres de van der Waals. J’ai donc extrait les 35 dimères possibles
à partir des multimères de la référence [338]. Pour chaque dimère, un calcul SAPT a été
effectué et à partir de cela, un nouveau jeu de paramètres a été obtenu :
– ε = 244565.0
– A1ii = 3.7604
– A2ii = 3.1182
La densité du modèle est de 0.84 g.cm 3 ce qui ne convient toujours pas. Cela provient
du manque d’interactions à courte distance étant donné que seuls les états d’équilibre
sont considérés pour ces dimères mais aussi de la faible statistique avec les 35 dimères
considérés.
Afin d’obtenir des paramètres de van der Waals plus représentatifs la méthode
employée ici est l’extraction des géométries de dimère issues de simulations de dynamique