Quelques exercices de Chimie Théorique

P.Reinhardt
Laboratoire de Chimie Théorique, Université Paris VI,
email: Peter.Reinhardt@upmc.fr

0.1  Li2 en base minimale, orienté sur l'axe z

Questions:
  1. Combien y'a-t-il d'OA? Lesquelles?
  2. Combien y'a-t-il d'OM? Donner leur dénomination (symétrie) et leur occupation dans l'état fondamental.
  3. Donner les contributions à l'énergie total Hartree-Fock. Ecrire l'énergie totale Hartree-Fock plus simple en utilisant l' opérateur de Fock.
  4. Ecrire l'énergie totale Hartree-Fock en n'utilisant que des sommes sur les AO. Quelle quantité est généralement introduite pour cela?
  5. Comment calcule-t-on l'énergie MP2? L'énergie ICSD?
  6. Rappeler le théorème de Brillouin. Comment des déterminants monoexcités peuvent-ils interagir avec l'état fondamental?
  7. STO-2G est une base minimale un peu ancienne. Que veut dire le sigle? Sachant qu'une fonction de Slater d'exposant 1 peut être approchée par la combinaison linéaire de deux Gaussiennes avec "exposants" 0.8518 et 0.1516 et coefficients 0.43 et 0.67, donner l'expression de deux Gaussiennes approchant la fonction de Slater d'exposant 1.625. Dites quel avantage il y a à utiliser la fonction approchée.
  8. Parmi les intégrales de recouvrement en OA suivantes, quelles sont celles qui sont strictement zéro: á1sA | 2sAñ,á 1sA | 2sBñ, á2pxA | 2sBñ, á2pzA | 2sBñ, á2pzA | 2sAñ, á2pyA | 2pxBñ. Dessiner schématiquement la matrice de recouvrement S.
  9. Comment pourrait-on organiser les OA pour concentrer les éléments de valeur zéro dans la matrice de recouvrement?

0.2  Molécules en C2v, utilisation de la symétrie de la molécule

La symétrie aide à économiser les ressources informatiques. Ce n'est pas seulement utile dans le sens pragmatique ou calculatoire. Prenons une molécule d'ozone:
  1. Ecrire les deux diagrammes extrêmes de Lewis pour cette molécule ayant une résonance dans les électrons p dans l'état fondamental.
  2. La symétrie est bien C2v, correspondant au tableau de caractères donné. Quelles orbitales moléculaires seront occupées par les électrons de cur?
  3. Parmi les orbitales a1, a2, b1 et b2 quelles sont celles qui appartiennent au système p?
  4. La configuration
    (c\oeur)(3a1)2(2b2)2(4a1)2(3b2)2(1b1)2(5a1)2(1a2)2 (4b2)2(2b1)2(6a1)0
    est-elle l'état fondamental de l'ozone? Sinon donner la configuration de l'état fondamental correct. Que faire pour y arriver?
L'eau étant une molécule très fréquemment rencontrée dans des applications numériques, il est utile de connaître quelques caractéristiques (molécule dans le plan yz, axe C2 sur z):
  1. Il y a 10 électrons dans cette molécule. Dessiner la structure Lewis, et trouver l'occupation des orbitales à symétrie adaptée.
  2. Quelle est (à peu près, en a.u.) l'énergie totale d'une molécule d'eau? Quel sera à peu près l'erreur de l'énergie totale, si on se trompait dans les unités de longueur (Ångstrøm ou unités atomiques)?
  3. L'ordre énergétique des orbitales est (1b2)2 (3a1)2 (1b1)2 (4a1)0 (2b2)0. L'orbitale 1b1 est-elle liante, anti-liante ou non-liante?
  4. On voudrait faire une interaction de configurations limitée aux excitations simples et doubles pour calculer l'énergie de l'état fondamental. Enumérer les configuration utiles dans l'espace des cinq orbitales données.
  5. Quels éléments de la matrice H faut-il calculer pour obtenir une énergie MP2?
Le tableau de caractères et les produits directs des représentations irréductibles du groupe C2v:

C2v
E
C2
sv(xy)
s¢v(yz)

A1
A2
B1
B2

a1
1
1
1
1
A1
A1
A2
B1
B2
a2
1
1
-1
-1
A2

A1
B2
B1
b1
1
-1
1
-1
B1


A1
A2
b2
1
-1
-1
1
B2



A1
Voilà un vecteur de l'IC simples et doubles, tel qu'il nous est donné par un programme de calcul.


0
0
0
0
1.000000000000

-2
0
-4
0
-0.004244665774
2
0
4
0
-0.004244665774
1
0
5
0
-0.003596647246
-1
0
-5
0
-0.003596647246

-1
-2
-4
-5
0.007107179853
-1
2
-4
5
-0.003223039016
-1
2
4
-5
-0.010330222102
1
-2
-4
5
-0.010330222102
1
-2
4
-5
-0.003223039016
1
2
4
5
0.007107179853

1
-1
4
-4
-0.008520757366
1
-1
5
-5
-0.020204959919
2
-2
4
-4
-0.010996418956
2
-2
5
-5
-0.004103218411
3
-3
4
-4
-0.008239966503
3
-3
5
-5
-0.001902303646

0.3  Orbitales et énergie

  1. Trouvez les expressions d'énergie des occupations suivantes, en utilisant h11, h22, J11, J12, J22 et K12: Figure

  2. Le système d'équations linéaires d'un IC de 2 électrons sur 2 orbitales (H2 en base minimale par exemple) s'écrit

    æ
    ç
    ç
    è
    0
    H01
    H01
    H11

    ö
    ÷
    ÷
    ø

    æ
    ç
    ç
    è
    1
    c1

    ö
    ÷
    ÷
    ø
      =  E  æ
    ç
    ç
    è
    1
    c1

    ö
    ÷
    ÷
    ø
           .
    En déduire comment le système d'équations linéaires s'écrit pour N molécules d'hydrogène sans interaction. Et en MP2, que faut-il écrire?
  3. Comment calcule-t-on des propriétes telles que le moment dipolaire d'une molécule, des fréquences de vibration, une analyse de population?
  4. Que dit le théorème de Koopmans?

0.4  Couplage de plusieurs spins

L'opérateur S2 du spin total d'un ensemble de spins a et b s'écrit ainsi:
S2 y  =  (h/2p)2 ì
í
î

å
P 
Pab + 1

4
[(na - nb)2+2 (na+nb) ] ü
ý
þ
y       .
Pour qu'une fonction de spin soit fonction propre de S2 il faut que
S2y  =  (h/2p)2S(S+1)y       .
Pour obtenir la contribution d'une fonction de spin quelconque correspondant à la valeur propre Sm, on peut faire une projection:
OSm  =  
Õ
l ¹ m 
[ S2 - (h/2p)2 Sl(Sl+1)]       .
  1. Construire les 4 fonctions propres (normées) de spin pour S=3/2 avec 3 spins.
  2. Combien de quartets (S=3/2) indépendants peut-on construire à partir de 5 spins?
  3. Quelles configurations possibles seraient intéressantes pour faire un ICSD d'une molécule (de plus que 2 électrons) à couche fermée?

0.5  Généralités

  1. Rappeler le principe de la DFT et les différences avec Hartree-Fock. Que signifie l'expression "LDA"?
  2. Que veut-on exprimer quand on parle d'un BSSE?
  3. Que signifie le mot "corrélation dynamique"? Et "corrélation non-dynamique"?
  4. Qu'apelle-t-on MCSCF?
  5. Voyons un exemple d'un calcul MCSCF, même si nous n'optimisons qu'un seul déterminant : nous démarrons avec quelques orbitales données (structure Lewis, orbitales construites "à la main", orbitales Hückel), et nous (1) orthogonalisons ces orbitales, puis construisons une matrice Fock. Ensuite nous diagonalisons une matrice d'un IC des monoexcitations (montrer que pour des déterminants adaptés au spin

    áF0|H|F0ñ 
    =
    E0


    áF0|H|Fiañ 
    =
    Ö2 Fia


    áFia|H|Fiañ 
    =
    E0 + Faa - F ii + ( a a | i i ) + ( a i | i a )


    áFia|H|Fjañ 
    =
    -Fij - ( a a | i j ) + ( a i | j a )


    áFia|H|Fibñ 
    =
    Fab - ( i i | a b ) + ( i a | b i )     ).

    Ceci baisse l'énergie totale. Ensuite nous prenons les coefficients cia de l'IC et nous modifions les orbitales occupées [(f)\tilde]i=fi+åa ciafa et les orbitales virtuelles [(f)\tilde]a=fa-åi ciafi. Montrer que Sia = á [(f)\tilde]i|[(f)\tilde]añ = 0. On revient sur (1) et on itère jusqu'à la convergence. Cette procédure nous amène à une solution Hartree-Fock. Pourquoi?



File translated from TEX by TTH, version 3.66.
On 12 Feb 2005, 12:40.