Conditions d'Eckart - Définition

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

Énergie de vibration

L'énergie de vibration d'une molécule peut être écrite en termes de coordonnées par rapport au repère d'Eckart comme :

 2T_\mathrm{vib} = \sum_{A=1}^N M_A \dot{\mathbf{R}}_A\cdot \dot{\mathbf{R}}_A = \sum_{A=1}^N M_A \dot{\mathbf{d}}_A\cdot \dot{\mathbf{d}}_A.

Le repère d'Eckart n'étant pas galiléen, l'énergie cinétique totale comprend aussi les termes de rotation et de Coriolis, que l'on considère comme implicites par la suite. L'énergie de vibration est écrite en termes de coordonnées de déplacment, qui sont linéairement dépendantes car influées par les 6 modes externes, nuls, i.e. les dA satisfont à 6 relations linéaires. Il est possible d'écrire l'énergie de vibration uniquement en termes de modes internes qr (r =1, ..., 3N-6) comme il va être montré par la suite. On écrit les différents modes en termes de déplacements

 \begin{align} q_r = \sum_{Aj} d_{Aj}& \big( q^A_{rj} \big) \\ s_i = \sum_{Aj} d_{Aj}& \big( M_A \delta_{ij} \big) =0 \\ s_{i+3} = \sum_{Aj} d_{Aj}& \big( M_A \sum_k \epsilon_{ikj} R^0_{Ak} \big)=0 \\ \end{align}

Les expressions entre parenthèses définissent une matrice B liant les modes internes et externes aux déplacements. la matrice B peut être divisée en une partie interne (3N-6 x 3N) et une partie externe (6 x 3N)

 \mathbf{v}\equiv  \begin{pmatrix} q_1 \\ \vdots \\ \vdots \\ q_{3N-6} \\ 0 \\ \vdots \\ 0\\ \end{pmatrix} = \begin{pmatrix} \mathbf{B}^\mathrm{int} \\ \cdots \\ \mathbf{B}^\mathrm{ext} \\ \end{pmatrix} \mathbf{d} \equiv \mathbf{B} \mathbf{d}.

On définit la matrice M par :

 \mathbf{M} \equiv \operatorname{diag}(\mathbf{M}_1, \mathbf{M}_2, \ldots,\mathbf{M}_N) \quad\textrm{et}\quad \mathbf{M}_A\equiv  \operatorname{diag}(M_A, M_A, M_A)

et à partir des relations données précédemment on obtient les relations matricielles :

 \mathbf{B}^\mathrm{ext} \mathbf{M}^{-1} (\mathbf{B}^\mathrm{ext})^\mathrm{T} = \operatorname{diag}(N_1,\ldots, N_6) \equiv\mathbf{N},

et

 \mathbf{B}^\mathrm{int} \mathbf{M}^{-1} (\mathbf{B}^\mathrm{ext})^\mathrm{T} = \mathbf{0}.

On définit

 \mathbf{G} \equiv \mathbf{B}^\mathrm{int} \mathbf{M}^{-1} (\mathbf{B}^\mathrm{int})^\mathrm{T}.

En utilisant les règles de multiplication des matrices par blocs, on montre que :

 (\mathbf{B}^\mathrm{T})^{-1} \mathbf{M} \mathbf{B}^{-1} = \begin{pmatrix} \mathbf{G}^{-1} && \mathbf{0} \\ \mathbf{0}      && \mathbf{N}^{-1} \end{pmatrix},

G-1 est de dimension (3N-6 x 3N-6) et N-1 de dimension (6 x 6). Lénergie cinétique devient :

 2T_\mathrm{vib} =  \dot{\mathbf{d}}^\mathrm{T}  \mathbf{M} \dot{\mathbf{d}} = \dot{\mathbf{v}}^\mathrm{T}\; (\mathbf{B}^\mathrm{T})^{-1} \mathbf{M} \mathbf{B}^{-1}\; \dot{\mathbf{v}} = \sum_{r, r'=1}^{3N-6} (G^{-1})_{r r'} \dot{q}_r \dot{q}_{r'}

où les six derniers composants sont nuls. Cette forme de l'énergie cinétique de vibration fait partie de la méthode GF de Wilson. Il est intéressant de signaler que l'énergie potentielle dans l'approximation harmonique peut être écrite de la façon suivante :

 2V_\mathrm{harm} =  \mathbf{d}^\mathrm{T} \mathbf{H} \mathbf{d} = \mathbf{v}^\mathrm{T} (\mathbf{B}^\mathrm{T})^{-1} \mathbf{H} \mathbf{B}^{-1} \mathbf{v} = \sum_{r, r'=1}^{3N-6} F_{r r'} q_r q_{r'},

H est le hessien du potentiel au minimum et F, défini par cette équation, est la matrice F de la méthode GF.

Page générée en 0.076 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