Hamiltonien moléculaire - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Hamiltonien du mouvement nucléaire de Watson

Afin d'obtenir un hamiltonien pour les mouvements externes (translation et rotation) couplés aux mouvements internes (vibration), il est usuel de faire un retour à ce niveau à la mécanique classique et de formuler lénergie cinétique classique correspondant à ces mouvements des noyaux. Il est facile de séparer de manière classique le mouvement de translation (du centre de masse) des autres mouvements. Cependant, séparer le mouvement de rotation du mouvement vibrationnel est plus difficile, et en réalité, pas entièrement possible. Cette séparation rovibrationelle fut pour la première fois réalisée par Eckart en 1935 en imposant ce qui est maintenant connu sous le nom de conditions d'Eckart. Lorsque le problème est décrit dans un référentiel (dit d'Eckart) qui tourne avec la molécule, et donc non galiléen, les énergies associées aux forces d'inertie (force centrifuge et force de Coriolis) apparaissent dans l'énergie cinétique.
En général, l'énergie cinétique classique T définit le tenseur métrique g = (gij) associé avec les coordonnées curvilignes s = (si) par la relation :

 2T = \sum_{ij} g_{ij} \dot{s}_i \dot{s}_j .

L'étape de quantification est la transformation de cette énergie cintéique classique en un opérateur de chimie quantique. Il est usuel de suivre la méthode proposée par Podolsky en réécrivant l'opérateur de Laplace-Beltrami dans les mêmes coordonnées s (généralisées, curvilignes) utilisées pour la forme classique. L'équation pour cet opérateur nécessite l'inversion du tenseur métrique g et son déterminant. La multiplication de l'opérateur de Laplace-Beltrami par -\hbar^2 procure l'opérateur d'énergie cinétique quantique requis. Lorsque l'on applique cette méthode aux coordonnées cartésiennes, qui ont une métrique unitaire, la même énergie cinétique est obtenue par application des règles de quantification.
Le hamiltonien de mouvement nucléaire a été décrit par Wilson et Howar en 1936, qui suivirent cette méthode, puis il fut redéfini plus tard par Darling et Dennison en 1940. Cette méthode resta le standard jursqu'en 1968, lorsque Watson fut capable de le simpflier de manière drastique en commutant par les dérivées le déterminant du tenseur métrique. On donnera par la suite le hamiltonien rovibrationnel obtenu par Watson, parfois appelé hamiltonien de Watson. On se doit de préciser que la dérivation de ce Hamiltonien est aussi possible en commençant depuis l'opérateur laplacien dans sa forme cartésienne, en appliquant les transformations de coordonnées, puis en utilisant les règles des dérivation des fonctions composées. Le hamiltonien de Watson, décrivant tous les mouvements des N noyaux, est :

  \hat{H} = -\frac{\hbar^2}{2M_\mathrm{tot}} \sum_{\alpha=1}^3 \frac{\partial^2}{\partial X_\alpha^2} +\frac{1}{2} \sum_{\alpha,\beta=1}^3 \mu_{\alpha\beta} (\mathcal{P}_\alpha - \Pi_\alpha)(\mathcal{P}_\beta - \Pi_\beta) +U  -\frac{\hbar^2}{2} \sum_{s=1}^{3N-6} \frac{\partial^2}{\partial q_s^2} + V .

Le premier terme est celui du centre de masse

 \mathbf{X} \equiv \frac{1}{M_\mathrm{tot}} \sum_{i=1}^N M_i \mathbf{R}_i \quad\mathrm{avec}\quad M_\mathrm{tot} \equiv \sum_{i=1}^N M_i.

Le second terme est le terme rotationnel semblable à l'énergie cinétique du rotateur rigide. Ici, \mathcal{P}_\alpha est la composante α de l' opérateur de moment angulaire du rotateur rigide fixe. L'opérateur \Pi_\alpha\, est une composante de l'opérateur dconnu sous le nom d'opérateur de moment angulaire vibrationnel (bien qu'il ne satisfasse pas aux relations de commutation du moment angulaire),

 \Pi_\alpha = -i\hbar \sum_{s,t=1}^{3N-6}  \zeta^{\alpha}_{st} \; q_s \frac{\partial}{\partial q_t}

avec la constante de couplage de Coriolis :

 \zeta^{\alpha}_{st} = \sum_{i=1}^N  \sum_{\beta,\gamma=1}^3 \epsilon_{\alpha\beta\gamma}   Q_{s, i\beta}\,Q_{t,i\gamma} \;\; \mathrm{et}\quad\alpha=1,2,3.

Ici, εαβγ est le symbole de Levi-Civita. Les termes quadratiques dans \mathcal{P}_\alpha sont des termes centrifuges, les termes bilinéaires dans \mathcal{P}_\alpha et \Pi_\beta\, sont des termes de Coriolis.
Les quantités Q s, iγ sont des composantes des coordonnées normales introduites ci-dessus. Les coordonnées normales peuvent être introduites de manière alternative par l'application de la méthode GF de Wilson. La matrice symétrique 3 x 3 \boldsymbol{\mu} est appelée tenseur effectif réciproque d'inertie. Si tous les qs étaient nuls (molécule rigide), le référentiel d'Eckart coinciderait avec le repère des axes principaux (voir rotateur rigide) et \boldsymbol{\mu} serait diagonale, avec les moments d'inertie réciproques d'équilibre sur la diagonale. Si tous les qs étaient nuls, seule les énergies cinétiques de translation et de rotation rigide seraient exprimées.
Le terme U de type potentiel est le terme de Watson :

 U = -\frac{1}{8} \sum_{\alpha=1}^3 \mu_{\alpha\alpha}

proportionnel à la trace du tenseur effectif réciproque d'inertie. Le quatrième terme dans le hamiltonien de Watson est l'énergie cinétique associée aux vibrations atomiques (nucléaires) exprimée en coordonnées normales qs, qui comme indiqué ci-dessus, sont données en termes de déplacements nucléaires ρ par :

 q_s = \sum_{i=1}^N \sum_{\alpha=1}^3 Q_{s, i\alpha} \rho_{i\alpha}\quad\mathrm{pour}\quad s=1,\ldots, 3N-6.

Enfin, V est l'énergie potentiel non développée par définition dépendante des seules coordonnées internes. Dans l'approximation harmonique, elle prend la forme :

 V \approx \frac{1}{2} \sum_{s=1}^{3N-6} f_s q_s^2.
Page générée en 0.096 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales
Version anglaise | Version allemande | Version espagnole | Version portugaise