Multipôles

L'extraction des multipôles est codée dans quatre modules supplémentaires :

• un module génère un réseaux de points au dessus de la surface à étudier où est evalué le potentiel du système périodique. Ces points servent de référence pour l'adaptations des multipôles;
• un deuxième module calcule les intégrales monoélectroniques entre multipôles et électrons à l'aide de l'algorithme de McMurchie et Davidson;[2]L. E. McMurchie, E. R. Davidson, J.Comput.Phys., 26 (1978) 218
• un troisième calcule le potentiel électrostatique de l'agrégat et des multipôles aux points générés par la première étape;
• finalement, l'adaptation des multipôles est effectuée par un dernier module à l'aide des valeurs du potentiel électrostatique calculées pour la surface idéale.

Figure 33: Schéma d'itération de l'extraction des multipôles.

Le potentiel sur la surface est calculé par CRYSTAL (version 92 et 98), de même les position des atomes (remplacés par multipôles) sont trouvées par CRYSTAL. La fonction d'onde de l'agrégat est calculé par le programme ORS de la branche moléculaire des programmes de calcul en orbitales localisées, mais peut être confié également à une ancienne version de TURBOMOLE.[3]Turbomol version 1.0, R. Ahlrichs, M. Häser, M. Bär, H. Horn, Chr. Kölmel, Université de Karlsruhe, 1989 Comme le potentiel créé par les multipôles entre dans la matrice Fock par des intégrales monoélectroniques supplémentaires, les étapes du calcul de la corrélation électronique (transformtion des intégrales biélectroniques, perturbation ou IC) sont disponibles sans modifications.