Méthode GF de Wilson

Méthode GF de Wilson

Méthode GF

La méthode GF de Wilson, parfois appelée méthode FG, est une méthode de mécanique classique d'obtention de certaines coordonnées internes pour une molécule semi-rigide vibrante, les coordonnées normales Qk. Les coordonnées normales décomposent les mouvements classiques de vibrations de la molécule et permettent ainsi d'obtenir des amplitudes vibrationnelles des atomes en fonction du temps. On postule, dans la méthode GF de Wilson, que l'énergie cinétique moléculaire est due seulement aux vibrations harmoniques des atomes, c'est-à-dire que l'énergie globale de rotation et translation est ignorée. Les coordonnées normales apparaissent aussi dans la description quantique des mouvements de vibration de la molécule et dans le couplage de Coriolis entre rotations et vibrations.

Il en découle qu'après application des conditions d'Eckart que la matrice G-1 donne l'énergie cinétique en termes de coordonnées internes linéaires arbitraires, alors que F représente l'énergie potentielle (harmonique) en fonction de ces coordonnées. La méthode GF donne la transformation linéaire des coordonnées internes générales à l'ensemble spécial de coordonnées normales.

Sommaire

La méthode GF

Une molécule non-linéaire de N atomes possède 3N-6 degrés de liberté internes, car localiser une molécule dans un espace tridimensionnel nécessite trois degrés de liberté et son orientation spatiale nécessite trois autres degrés de liberté. Ces degrés de liberté doivent être soustraits des 3N degrés de liberté d'un système de N particules.

Les atomes d'une molécule sont « liés » par une surface d'énergie potentielle (SEP) (ou un champ de force dans une approche classique)[1] qui est une fonction de 3N-6 coordonnées.

Les degrés internes de liberté q1, ..., q3N-6 décrivant la SEP de façon optimale sont parfois non-linéaires, comme les coordonnées de valence telles que les angles de liaison, de torsion ou les étirements de liaisons. Il est possible d'écrire l'opérateur d'énergie cinétique quantique pour de telles coordonnées curvilignes, mais il est difficile de formuler une théorie générale applicable à toute molécule. C'est pourquoi E.B. Wilson[2] linéarisa les coordonnées internes en postulant des déplacements minimes. La version linéarisée de la coordonnée interne qt est notée St.

La SEP V peut être développée en série de Taylor au voisinage de son minimum en termes de St. Le troisième terme (le hessien de V), évalué au minimum, est une matrice 'F dérivée de la force. Dans l'approximation harmonique, la série de Taylor s'achève avec ce terme. Le deuxième terme, contenant les dérivées premières, est nul car évalué au minimum de V. Le premier terme peut être inclus dans le zéro d'énergie. D'où :

2V \approx \sum_{s,t=1}^{3N-6} F_{st} S_s\, S_t .

L'énergie cinétique de vibration classique a la forme :

2T = \sum_{s,t=1}^{3N-6} g_{st}(\mathbf{q}) \dot{S}_s\dot{S}_t ,

où gst est un élément du tenseur métrique des coordonnées internes (curvilignes). Les points indiquent les dérivées temporelles. L'évaluation du tenseur métrique g au minimum q0 de V donne la matrice définie positivement et symétrique G = g(q0)-1. On peut alors résoudre les deux problèmes matriciels suivants simultanément :

\mathbf{L}^\mathrm{T} \mathbf{F} \mathbf{L} =\boldsymbol{\Phi} \quad \mathrm{et}\quad \mathbf{L}^\mathrm{T} \mathbf{G}^{-1} \mathbf{L} = \mathbf{E},

car ils sont équivalents au problème aux valeurs propres généralisé :

\mathbf{G} \mathbf{F} \mathbf{L} = \mathbf{L} \boldsymbol{\Phi},

où \boldsymbol{\Phi}=\operatorname{diag}(f_1,\ldots, f_{3N-6}) et \mathbf{E}\, est la matrice unitaire. La matrice L-1 contient les coordonnées normales Qk dans ses lignes :

Q_k = \sum_{t=1}^{3N-6} (\mathbf{L}^{-1})_{kt} S_t , \quad k=1,\ldots, 3N-6. \,

En raison de la forme du problème aux valeurs propres généralisé, la méthode est appelée méthode GF, appellation à laquelle on accole parfois le nom de son inventeur : la méthode GF de Wilson.

On introduit les vecteurs

\mathbf{s} = \operatorname{col}(S_1,\ldots, S_{3N-6}) \quad\mathrm{et}\quad \mathbf{Q} = \operatorname{col}(Q_1,\ldots, Q_{3N-6}),

satisfaisant la relation :

\mathbf{s} = \mathbf{L} \mathbf{Q}.

En introduisant les résultats de l'équation aux valeurs propres, l'énergie de la molécule E = T + V (dans l'approximation harmonique) devient :

2E = \dot{\mathbf{s}}^\mathrm{T} \mathbf{G}^{-1}\dot{\mathbf{s}}+ \mathbf{s}^\mathrm{T}\mathbf{F}\mathbf{s}
= \dot{\mathbf{Q}}^\mathrm{T} \; \left( \mathbf{L}^\mathrm{T} \mathbf{G}^{-1} \mathbf{L}\right) \; \dot{\mathbf{Q}}+ \mathbf{Q}^\mathrm{T} \left( \mathbf{L}^\mathrm{T}\mathbf{F}\mathbf{L}\right)\; \mathbf{Q}
= \dot{\mathbf{Q}}^\mathrm{T}\dot{\mathbf{Q}} + \mathbf{Q}^\mathrm{T}\boldsymbol{\Phi}\mathbf{Q} = \sum_{t=1}^{3N-6} \big( \dot{Q}_t^2 + f_t Q_t^2 \big).

Le lagrangien L = T - V est :

L = \frac{1}{2} \sum_{t=1}^{3N-6} \big( \dot{Q}_t^2 - f_t Q_t^2 \big).

Les équations de Lagrange sont identiques aux équations de Newton :

\ddot{Q}_t + f_t \,Q_t = 0

pour un ensemble d'oscillateurs harmoniques non couplés. Ces équations différentielles ordinaires du second ordre sont facilement solvables, et donnent Qt en fonction du temps (voir l'article sur l'oscillateur harmonique).

Coordonnées normales en termes de coordonnées cartésiennes de déplacement

Les coordonnées normales sont parfois exprimées comme des combinaisons linéaires des coordonnées cartésiennes de déplacement. Soit RA le vecteur position du noyau A et RA0 la position à l'équilibre correspondante. \mathbf{d}_A \equiv \mathbf{R}_A -\mathbf{R}_A^0 est par définition la coordonnée cartésienne de déplacement du noyau A. La linéarisation de Wilson des coordonnées curvilignes internes qt exprime la coordonnée St en termes de coordonnées de déplacement :

S_t =\sum_{A=1}^N \sum_{i=1}^3 s^t_{Ai} \, d_{Ai}= \sum_{A=1}^N \mathbf{s}^t_{A} \cdot \mathbf{d}_{A}, \quad \mathrm{pour}\quad t = 1,\ldots,3N-6,

où sAt est le vecteur s de Wilson. Si l'on introduit s^t_{Ai} dans la matrice 3N-6 x 3N B, cette équation devient en notation matricielle :

\mathbf{s} = \mathbf{B} \mathbf{d}.

La forme réelle des éléments de matrice de B peut être relativement compliquée. Pour un angle de torsion en particulier, impliquant 4 atomes, on a besoin d'une algèbre vectorielle fastidieuse afin d'en déduire les valeurs correspondantes des s^t_{Ai}. On pourra, pour plus d'informations sur cette méthode connue sous le nom de méthode des vecteurs s de Wilson se reporter à l'ouvrage de Wilson et collaborateurs, ou à l'article vibration moléculaire. À ce stade,

\mathbf{Q} = \mathbf{L}^{-1} \mathbf{s} = \mathbf{L}^{-1} \mathbf{B} \mathbf{d} \equiv \mathbf{D} \mathbf{d}.

Sous forme de somme :

Q_k = \sum_{A=1}^N \sum_{i=1}^3 D^k_{Ai}\, d_{Ai} \quad \mathrm{pour}\quad k=1,\ldots, 3N-6.

D est ici une matrice 3N-6 x 3N donnée par la linéarisation des coordonnées internes q (procédé algébrique) et la solution des équations GF de Wilson (procédé numérique).

Relations aux conditions d'Eckart

Article principal : Conditions d'Eckart.

A partir de l'invariance des coordonnées internes St par rapport à toute rotation globale et translation de la molécule, on en déduit que les coordonnées linéarisées stA présentent les mêmes propriétés.On peut montrer que cela implique que les six conditions suivantes sont satisfaites par les coordonnées internes :

\sum_{A=1}^N \mathbf{s}^t_{A} = 0\quad\mathrm{et}\quad \sum_{A=1}^N \mathbf{R}^0_A\times \mathbf{s}^t_A= 0, \quad t=1,\ldots,3N-6.

Ces conditions découlent des conditions d'Eckart qui sont valables pour les vecteurs déplacement :

\sum_{A=1}^N M_A\; \mathbf{d}_{A} = 0 \quad\mathrm{et}\quad \sum_{A=1}^N M_A\; \mathbf{R}^0_{A} \times \mathbf{d}_{A} = 0.

Notes

  1. ↑ De manière générale, la SEP décrit les variations d'énergie induites par les différentes interactions entre atomes. Les positions atomiques y sont localisées dans des « vallées » où ces énergies sont minimales, et d'où il peut être difficile de s'extraire sauf apport extérieur d'énergie (excitations électromagnétiques ou thermiques, etc.).
  2. ↑ E. B. Wilson, Jr. Some Mathematical Methods for the Study of Molecular Vibrations, J. Chem. Phys. vol. 9, pp. 76-84 (1941)

Références

  • (en) E. B. Wilson, J. C. Decius, et P. C. Cross, Molecular Vibrations, McGraw-Hill, New York, 1955 (réimprimé par Dover, 1980).
  • (en) D. Papoušek et M. R. Aliev, Molecular Vibrational-Rotational Spectra Elsevier, Amsterdam, 1982.
  • (en) S. Califano, Vibrational States, Wiley, London, 1976.


  • (en) Cet article est partiellement ou en totalité issu d’une traduction de l’article de Wikipédia en anglais intitulé « GF method ».
  • Portail de la chimie Portail de la chimie
  • Portail de la physique Portail de la physique
Ce document provient de « M%C3%A9thode GF ».

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Méthode GF de Wilson de Wikipédia en français (auteurs)

Игры ⚽ Нужно решить контрольную?

Regardez d'autres dictionnaires:

  • Méthode FG — Méthode GF La méthode GF de Wilson, parfois appelée méthode FG, est une méthode de mécanique classique d obtention de certaines coordonnées internes pour une molécule semi rigide vibrante, les coordonnées normales Qk. Les coordonnées normales… …   Wikipédia en Français

  • Méthode GF — La méthode GF de Wilson, parfois appelée méthode FG, est une méthode de mécanique classique d obtention de certaines coordonnées internes pour une molécule semi rigide vibrante, les coordonnées normales Qk. Les coordonnées normales décomposent… …   Wikipédia en Français

  • WILSON (E.) — Pendant un demi siècle, Wilson a dominé intellectuellement la scène littéraire aux États Unis, interprétant pour ses contemporains les multiples aspects des grands courants artistiques, sociaux et politiques qui ont modifié le visage de… …   Encyclopédie Universelle

  • Méthode composite de chimie quantique — Les méthodes composites de chimie quantique sont des méthodes de chimie numérique dont le but est d obtenir une grande précision par combinaison des résultats de plusieurs calculs. Cette combinaison est constituée à partir de méthodes d un haut… …   Wikipédia en Français

  • Wilson-Bappu-Effekt — Der Wilson Bappu Effekt ist ein physikalisch astronomischer Effekt, der zur Distanzbestimmung von Sternen benutzt wird. Im Jahre 1957 berichteten Olin C. Wilson und M. K. Vainu Bappu die bemerkenswerte Korrelation zwischen der Breite der… …   Deutsch Wikipedia

  • Kenneth Wilson — Kenneth Geddes Wilson (* 8. Juni 1936 in Waltham, Massachusetts) ist ein US amerikanischer Physiker und Nobelpreisträger. Wilson, der Sohn von Edgar Bright Wilson, studierte als Putnam Fellow an der Universität Harvard und promovierte 1961 am… …   Deutsch Wikipedia

  • R.A. Wilson — Robert Anton Wilson (eigentlich Robert Edward Wilson; * 18. Januar 1932 in Brooklyn; † 11. Januar 2007 in Capitola bei Santa Cruz, Kalifornien) war ein US amerikanischer Bestsellerautor, Philosoph und Anarchist. Berühmt wurde er unter anderem… …   Deutsch Wikipedia

  • R. A. Wilson — Robert Anton Wilson (eigentlich Robert Edward Wilson; * 18. Januar 1932 in Brooklyn; † 11. Januar 2007 in Capitola bei Santa Cruz, Kalifornien) war ein US amerikanischer Bestsellerautor, Philosoph und Anarchist. Berühmt wurde er unter anderem… …   Deutsch Wikipedia

  • Robert Anton Wilson — (1977) Robert Anton Wilson (eigentlich Robert Edward Wilson; * 18. Januar 1932 in Brooklyn; † 11. Januar 2007 in Capitola bei Santa Cruz, Kalifornien) war ein US amerikanischer Bestsellerautor, Philosoph und Anarchist. Berühmt wurde er unter… …   Deutsch Wikipedia

  • Albright et Wilson — La compagnie Albright et Wilson (A W), fondée en 1856 au Royaume Uni par deux entrepreneurs quakers ; Arthur Albright et John Edward Wilson. Elle fabriquait du chlorate de potassium et du phosphore blanc pour l industrie des allumettes.… …   Wikipédia en Français

Share the article and excerpts

Direct link
… Do a right-click on the link above
and select “Copy Link”