|
Même si on est débutant , il vaut mieux apprendre à classer et noter ses exercices dans un recueil.
En voici un . Il a été prévu pour interroger à l'examen (cf [[Mécanique, enseignée via l'Histoire des Sciences|Mécanique , enseignée via l'Histoire des Sciences]] , remarques-en-vrac)
== qq Exercices d'électrostatique ==
=== Moment octupolaire d'un cube ===
Soit un cube de Na+ etCl- de côté 2a : montrer que loin du cube le poentiel (à 1/ 4Pi epsilon près) s'écrit :
V(x,y,z) = 120 xyz a^3/ r^7 .
Solution 1 : choisir dans les orbitales f celles qui ont la symétrie du cube : cela donne tout de suite le résultat xyz/ r^7 : il reste à trouver le facteur 120 a^3 . On peut alors développer les 8 termes +(ou-) q/AiM selon la formule usuelle de Legendre : il ne restera que les huit termes :
<math>{5 \over 2} (\vec{OA_i}\cdot \vec{OM})^3</math>
soit f(z) = +(x+y+z)^3 - (-x +y+z)^3 + (-x-y +z)^3 - (x-y+z)^3 et -f(-z)
donc un polynôme du 3eme degré impair en z , mais aussi en x et en y , donc 48 xyz (fois 5/2) =120xyz
Solution 2 : considérer la formule générale de la distribution dipolaire V(D') = E(D).translation.
On part du dipole 2a <math>\vec{}I</math> , on translate selon J de 2a , puis selon K de 2a .
On n'aura donc besoin que de la composante J du champ du dipole soit :E_y = 2a 3xy/r^5 d'où le potentiel du quadrupole : V =2a.E_y = 2a .2a 3xy/r^5 ; on n'a besoin que de la composante K du champ correspondant, soit 2a.2a 3.5 xyz/r^7 d'où le potentiel final : (2a)^3 . 3.5 . xyz/r^7.(solution très élégante de Melba Philips).
=== Constante de Madelung ===
Il existe un superbe article de Tissier dans la RMS.Ici, on se contente des méthodes classiques.
Calculer la cste de Madelung de la chaîne linéaire d'une mole de NaCl :
Solution : L'énergie est donc du type (en négligeant les effets de bords) 1/2 sigma (qiVi) =
1/2 2N (qo Vo) = - Ne^2/a. 2 .(1-1/2 +1/3 -1/4 +) = -Ne^2/(2a) . 4 Ln2 (Attention : certains prennent la période du cristal , soit 2a , et d'autres la distance Na---Cl = a).
Calculer la cste de Madelung du cristal cubique de NaCl : -Ne^2/a .(1.7476)
Solution : cette fois le décompte est moins simple : le potentiel Vo au centre sera :
on compte les carrés des distances de plus en plus grandes et les ions
*1 (1,0,0) : +6 (n'existe pas pour la boîte cubique)
*2 (1,1,0) : -12 (idem)
*3 (1,1,1) : +8
*4 (2,0,0) : -6 (idem)
*5 (2,1,0) : +24 (idem)
*6 (2,1,1) : -24
*7 : 0 (regarder pourquoi)
*8 (2,2,0) : -12
*9 (2,2,1) : +24 et (3,0,0) : +6 (dégénérescence , qui n'existe pas pour la particule dans la boîte cubique)
*10 (3,1,0): -24(idem)
*11 (3,1,1): +24
*12 (2,2,2): -8
Vo = -e^2/a [ 6/1 -12/sqrt(2) +8/sqrt(3) -6/2 +24/sqrt(5) -24/sqrt(6)-12/2sqrt(2) + ...]: désespéré.
d'où
==== la méthode d'Evjen ====
(comptage par cellule "neutre"): -1.212 ; puis -1.751 ; c'est plus encourageant.
les trois résultats importants sont :
#CsCl = -1.7627
#NaCl = -1.7476
#ZnS = -1.6381
==== la méthode d'Ewald ====
'''Beaucoup plus sophistiquée''' :l'idée fondamentale est de remplacer chaque charge ponctuelle par une gaussienne, et de considérer que le potentiel d'une gaussienne est q/r .erf(r/sqrt(2<math>\sigma</math>)qui a le bon comportement à longue distance et donc prendra bien en compte la sommation (si on sait la faire) aux grandes distances (donc la lenteur de la sommation aura disparu). Et bien sûr comme erf(x) + erfc(x) = 1 , on prendra soin du comportement à courte distance par le calcul avec les erfc(x).
Tout ceci exige une bonne connaissance des transformées de fourier et de la relation de fermeture des ondes planes. Moyennant quoi (cf exo d'électrostatique correspondant), le calcul (reférence par exemple Grosso & Parravicini)donne une sommation sur les cations et sur les anions avec erfc :
V1(0) = <math>\Sigma_{ca} {e^2 \over {|\vec a + \vec t_i|}} erfc({{|\vec a+\vec t_i|} \over \sqrt(2)\sigma}) -\Sigma_{\vec t_i non nul}{e^2 \over |t_i| }erfc ({|t_i| \over \sqrt(2)\sigma})</math>
et une autre sommation sur les longues distances , mais qu'on transforme astucieusement via la formule de Poisson en une sommation sur le réseau réciproque (donc sur de courts nombres d'ondes) :
V2(0) = <math>{-2e^2 \over \sqrt{2\pi \sigma^2}}</math> + <math>{4\pi e^2\over \Omega} </math><math>\Sigma_{\vec{k_i}non nul} (1- e^{-i \vec k_i \vec a}){1 \over k_i^2}e^{-k_i^2\cdot \sigma^2/2}</math>.
Un bon exercice est de vérifier le premier exercice (celui de la ligne : 2 Ln2): bien du plaisir !
=== Electrostatique de la charge d'Ewald ===
Pour mieux comprendre l'exercice précédent, une petite mise en jambe est nécessaire :
Se rappeler au moins la formule fondamentale du cours sur les gaussiennes :
<math><e^{ikx}> = <cos kx > = e^{-k^2 \over 2}</math>
*Exercice : montrer qu'une charge à l'origine , gaussienne ( de largeur <math>\sigma</math>
) a pour potentiel <math>V(r) = q/r \cdot erf(r/\sqrt{2 \sigma^2})</math>, puis vérifier le théorème de Gauss et les théorèmes sur l'énergie.
*Exercice : prendre un atome d'Ewald : charge +q à l'origine et électron en orbitale gaussienne :
comparer attentivement avec l'atome 1s . Pourquoi la méthode des orbitales gaussiennes a-t-elle eu un succès aussi immense?
== Mécanique quantique, rudiments ==
| 