Next: Etat actuel
Up: Un logiciel de calcul
Previous: Dimères
Contents
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.
Next: Etat actuel
Up: Un logiciel de calcul
Previous: Dimères
Contents
Peter Reinhardt
2004-10-05