GF method: Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Paul Wormer
(Got some valid criticism on WP version, tried to answer that)
imported>Paul Wormer
(Some further clarification of the method (after WP comments))
Line 1: Line 1:
{{subpages}}
{{subpages}}


In [[chemistry]] and [[molecular physics]], Wilson's '''GF method''', sometimes referred to as the '''FG method''', is a classical mechanical  method to obtain certain ''internal coordinates'' for a vibrating [[semi-rigid molecule]], the so-called ''normal coordinates'' ''Q''<sub>k</sub> (also known as normal modes). Simultaneously, the vibrational energies associated with these normal modes are obtained.
In [[chemistry]] and [[molecular physics]], Wilson's '''GF method''', sometimes referred to as the '''FG method''', is a classical mechanical  method to obtain certain ''internal coordinates'' for a vibrating [[semi-rigid molecule]], the so-called ''normal coordinates'' ''Q''<sub>k</sub> (also known as normal modes). Simultaneously, the vibrational energies associated with these normal modes are obtained.
   
   
Normal coordinates decouple the classical vibrational motions of the molecule making it possible to compute the vibrational amplitudes of atoms in a molecule.  In Wilson's GF method it is assumed that the molecular kinetic energy consists only of  vibrations of the atoms, ''i.e.,''  the overall rotational and translational energies of the molecule are ignored. Normal coordinates are usually defined  in a classical mechanical context, but they also appear in the quantum mechanical descriptions of the  vibrational motions of molecules and the Coriolis coupling between rotations and vibrations.
Normal coordinates decouple the classical vibrational motions of the molecule making it possible to compute the vibrational amplitudes of atoms in a molecule.  In Wilson's GF method it is assumed that the molecular kinetic energy consists only of  vibrations of the atoms, ''i.e.,''  the overall rotational and translational energies of the molecule are ignored. Normal coordinates are usually defined  in a classical mechanical context, but they also appear in the quantum mechanical descriptions of the  vibrational motions of molecules and the Coriolis coupling between rotations and vibrations.
Line 14: Line 14:


The atoms in a molecule are bound by a [[potential energy surface]] (PES) (or a [[force field (chemistry)|force field]]) which is a function of 3''N''-6 coordinates.
The atoms in a molecule are bound by a [[potential energy surface]] (PES) (or a [[force field (chemistry)|force field]]) which is a function of 3''N''-6 coordinates.
The internal degrees of freedom ''q''<sub>1</sub>, ...,  ''q''<sub>3''N''-6</sub> describing the PES in an optimum way are often non-linear; they are for instance ''valence coordinates'', such as bending and  torsion angles and bond stretches. It is possible to write the quantum mechanical kinetic energy operator for such [[curvilinear coordinates]], but it is hard to formulate a general theory applicable to any molecule. This is why Wilson<ref>E. B. Wilson, Jr.
The internal degrees of freedom ''q''<sub>1</sub>, ...,  ''q''<sub>3''N''-6</sub> describing the PES in an optimum way are often non-linear; they are for instance ''valence coordinates'', such as bending and  torsion angles and bond stretches. It is possible to write the quantum mechanical kinetic energy operator for such [[curvilinear coordinates]], but it is hard to formulate a general theory applicable to any molecule. This is why Wilson<ref name="Wilson41"/> linearized the internal coordinates by assuming small displacements. The linearized version of  the internal coordinate ''q''<sub>''t''</sub> is denoted by ''S''<sub>''t''</sub>.  
''Some Mathematical Methods for the Study of Molecular Vibrations'', J. Chem. Phys. vol. '''9''', pp. 76-84 (1941)</ref> linearized the internal coordinates by assuming small displacements. The linearized version of  the internal coordinate ''q''<sub>t</sub> is denoted by ''S''<sub>t</sub>.  


The PES ''V'' can be Taylor expanded around its minimum in terms of the ''S''<sub>t</sub>.  The third term (the [[Hessian matrix|Hessian]] of ''V'') evaluated in the minimum is a force derivative matrix '''F'''. In the harmonic approximation the Taylor series is ended after this term. The second term, containing first derivatives, is zero because it is evaluated in the miminum of ''V''. The first term can be included in the zero of energy.
The PES ''V'' can be [[Taylor series|Taylor-expanded]] around its minimum in terms of the ''S''<sub>''t''</sub>.  The third term (the [[Hessian matrix|Hessian]] of ''V'') evaluated in the minimum is a force derivative matrix '''F'''. In the ''harmonic approximation''&mdash;on which the method is based&mdash;the Taylor series is ended after this term. The second term, containing first derivatives, is zero because it is evaluated in the miminum of ''V''. The first term can be included in the zero of energy.
Thus,
Thus,
:<math> 2V \approx \sum_{s,t=1}^{3N-6} F_{st} S_s\, S_t </math>.
:<math> 2V \approx \sum_{s,t=1}^{3N-6} F_{st} S_s\, S_t </math>.
Line 23: Line 22:
The classical vibrational kinetic energy has the form:
The classical vibrational kinetic energy has the form:
:<math> 2T = \sum_{s,t=1}^{3N-6} g_{st}(\mathbf{q})  \dot{S}_s\dot{S}_t ,</math>
:<math> 2T = \sum_{s,t=1}^{3N-6} g_{st}(\mathbf{q})  \dot{S}_s\dot{S}_t ,</math>
where ''g''<sub>''st''</sub> is an element of the metric tensor of the internal (curvilinear) coordinates. The dots indicate [[time derivative]]s.  Evaluation of the metric tensor '''g''' in the minimum '''q'''<sup>0</sup> of ''V'' gives  
where ''g''<sub>''st''</sub> is an element of the [[metric tensor]] of the internal (curvilinear) coordinates. The dots indicate [[time derivative]]s.  Evaluation of the metric tensor '''g''' in the minimum '''q'''<sup>0</sup> of ''V'' gives  
the positive definite and symmetric matrix '''G''' = '''g'''('''q'''<sup>0</sup>)'''<sup>-1</sup>.  
the [[positive define matrix|positive definite]] and [[symmetric matrix]] '''G''' = '''g'''('''q'''<sup>0</sup>)'''<sup>-1</sup>.  
One can solve the following two matrix problems simultaneously
One can solve the following two matrix problems simultaneously
:<math> \mathbf{L}^\mathrm{T} \mathbf{F} \mathbf{L} =\boldsymbol{\Phi}
:<math> \mathbf{L}^\mathrm{T} \mathbf{F} \mathbf{L} =\boldsymbol{\Phi}
\quad \mathrm{and}\quad \mathbf{L}^\mathrm{T} \mathbf{G}^{-1} \mathbf{L} = \mathbf{E},
\quad \mathrm{and}\quad \mathbf{L}^\mathrm{T} \mathbf{G}^{-1} \mathbf{L} = \mathbf{E},
</math>
</math>
since they are equivalent to the generalized eigenvalue problem
since they are equivalent to the [[generalized eigenvalue problem]]
:<math>
:<math>
\mathbf{G} \mathbf{F} \mathbf{L} = \mathbf{L} \boldsymbol{\Phi},
\mathbf{G} \mathbf{F} \mathbf{L} = \mathbf{L} \boldsymbol{\Phi},
Line 38: Line 37:
often with the name of its originator attached to it: '''Wilson's GF method'''. By matrix transposition in both sides of the equation and using the fact
often with the name of its originator attached to it: '''Wilson's GF method'''. By matrix transposition in both sides of the equation and using the fact
that both '''G''' and '''F''' are symmetric matrices, as are diagonal matrices, one can recast this equation into a very similar one  for '''FG''' . This is why the method
that both '''G''' and '''F''' are symmetric matrices, as are diagonal matrices, one can recast this equation into a very similar one  for '''FG''' . This is why the method
is also referred to as '''Wilson's FG method'''.
is also referred to as '''Wilson's FG method'''. In [[numerical mathematics]] several algorithms are known to solve  generalized eigenvalue problems. The solution yields the matrix '''L''' (its inverse yields the normal modes ''Q''<sub>''k''</sub>)  and the diagonal matrix '''&Phi;''' which has the vibrational [[force constant]]s on the diagonal.  


We  introduce  the vectors
We  introduce  the vectors
Line 65: Line 64:
  = \sum_{t=1}^{3N-6} \big( \dot{Q}_t^2 +  f_t Q_t^2 \big).
  = \sum_{t=1}^{3N-6} \big( \dot{Q}_t^2 +  f_t Q_t^2 \big).
</math>
</math>
Solution of the GF equations yields ''Q''<sub>''t''</sub> and ''f''<sub>''t''</sub>. In the harmonic approximation the time dependence of the [[harmonic oscillator]]s is given by sine-type functions, so that the time derivatives of the ''Q''<sub>''t''</sub>'s are simple, and hence the  3''N''&minus;6 vibrational (kinetic plus potential) energies are known once the GF equations are solved.
The Lagrangian ''L'' = ''T'' - ''V'' is
The Lagrangian ''L'' = ''T'' - ''V'' is
:<math>  
:<math>  
Line 74: Line 75:
</math>
</math>
for a set of uncoupled
for a set of uncoupled
harmonic oscillators. These ordinary  second-order differential equations are easily solved yielding ''Q''<sub>''t''</sub> as function of time, see the article on  [[harmonic oscillator]]s.
harmonic oscillators. These ordinary  second-order differential equations are easily solved yielding ''Q''<sub>''t''</sub> as sine-type functions of time, see the article on  [[harmonic oscillator]]s.


==Normal coordinates in terms of Cartesian displacement coordinates ==
==Normal coordinates in terms of Cartesian displacement coordinates ==

Revision as of 05:03, 20 January 2009

This article is developing and not approved.
Main Article
Discussion
Related Articles  [?]
Bibliography  [?]
External Links  [?]
Citable Version  [?]
 
This editable Main Article is under development and subject to a disclaimer.

In chemistry and molecular physics, Wilson's GF method, sometimes referred to as the FG method, is a classical mechanical method to obtain certain internal coordinates for a vibrating semi-rigid molecule, the so-called normal coordinates Qk (also known as normal modes). Simultaneously, the vibrational energies associated with these normal modes are obtained.

Normal coordinates decouple the classical vibrational motions of the molecule making it possible to compute the vibrational amplitudes of atoms in a molecule. In Wilson's GF method it is assumed that the molecular kinetic energy consists only of vibrations of the atoms, i.e., the overall rotational and translational energies of the molecule are ignored. Normal coordinates are usually defined in a classical mechanical context, but they also appear in the quantum mechanical descriptions of the vibrational motions of molecules and the Coriolis coupling between rotations and vibrations.

It follows from application of the Eckart conditions that the matrix G−1 gives the kinetic energy in terms of arbitrary linear internal coordinates, while F represents the harmonic potential energy in terms of these internal coordinates. Solution of the GF equations gives the matrix that transforms from general internal coordinates to the special set of normal coordinates. The associated eigen energies of the vibrations are obtained along the way.

The method is called after the two matrices F and G that enter the final equations. The name of the most important of its early investigators[1] E. Bright Wilson, Jr. (1908–1992) is often attached to the method.

The GF method

A non-linear molecule consisting of N atoms has 3N-6 internal degrees of freedom, because positioning a molecule in three-dimensional space requires three degrees of freedom and the description of its orientation in space requires another three degree of freedom. These degrees of freedom must be subtracted from the 3N degrees of freedom of a system of N particles.

The atoms in a molecule are bound by a potential energy surface (PES) (or a force field) which is a function of 3N-6 coordinates. The internal degrees of freedom q1, ..., q3N-6 describing the PES in an optimum way are often non-linear; they are for instance valence coordinates, such as bending and torsion angles and bond stretches. It is possible to write the quantum mechanical kinetic energy operator for such curvilinear coordinates, but it is hard to formulate a general theory applicable to any molecule. This is why Wilson[1] linearized the internal coordinates by assuming small displacements. The linearized version of the internal coordinate qt is denoted by St.

The PES V can be Taylor-expanded around its minimum in terms of the St. The third term (the Hessian of V) evaluated in the minimum is a force derivative matrix F. In the harmonic approximation—on which the method is based—the Taylor series is ended after this term. The second term, containing first derivatives, is zero because it is evaluated in the miminum of V. The first term can be included in the zero of energy. Thus,

.

The classical vibrational kinetic energy has the form:

where gst is an element of the metric tensor of the internal (curvilinear) coordinates. The dots indicate time derivatives. Evaluation of the metric tensor g in the minimum q0 of V gives the positive definite and symmetric matrix G = g(q0)-1. One can solve the following two matrix problems simultaneously

since they are equivalent to the generalized eigenvalue problem

where and is the unit matrix. The matrix L-1 contains the normal coordinates Qk in its rows:

Because of the form of the generalized eigenvalue problem, the method is called the GF method, often with the name of its originator attached to it: Wilson's GF method. By matrix transposition in both sides of the equation and using the fact that both G and F are symmetric matrices, as are diagonal matrices, one can recast this equation into a very similar one for FG . This is why the method is also referred to as Wilson's FG method. In numerical mathematics several algorithms are known to solve generalized eigenvalue problems. The solution yields the matrix L (its inverse yields the normal modes Qk) and the diagonal matrix Φ which has the vibrational force constants on the diagonal.

We introduce the vectors

which satisfy the relation

Upon use of the results of the generalized eigenvalue equation the energy E = T + V (in the harmonic approximation) of the molecule becomes,

Solution of the GF equations yields Qt and ft. In the harmonic approximation the time dependence of the harmonic oscillators is given by sine-type functions, so that the time derivatives of the Qt's are simple, and hence the 3N−6 vibrational (kinetic plus potential) energies are known once the GF equations are solved.

The Lagrangian L = T - V is

The corresponding Lagrange equations are identical to the Newton equations

for a set of uncoupled harmonic oscillators. These ordinary second-order differential equations are easily solved yielding Qt as sine-type functions of time, see the article on harmonic oscillators.

Normal coordinates in terms of Cartesian displacement coordinates

Often the normal coordinates are expressed as linear combinations of Cartesian displacement coordinates. Let RA be the position vector of nucleus A and RA0 the corresponding equilibrium position. Then is by definition the Cartesian displacement coordinate of nucleus A. Wilson's linearizing of the internal curvilinear coordinates qt expresses the coordinate St in terms of the displacement coordinates

where sAt is known as a Wilson s-vector. If we put the into a 3N-6 x 3N matrix B, this equation becomes in matrix language

The actual form of the matrix elements of B can be fairly complicated. Especially for a torsion angle, which involves 4 atoms, it requires tedious vector algebra to derive the corresponding values of the . See for more details on this method, known as the Wilson s-vector method, the book by Wilson et al., or molecular vibration. Now,

In summation language:

Here D is a 3N-6 x 3N matrix which is given by (i) the linearization of the internal coordinates q (an algebraic process) and (ii) solution of Wilson's GF equations (a numeric process).

Relation with Eckart conditions

From the invariance of the internal coordinates St under overall rotation and translation of the molecule, follows the same for the linearized coordinates stA. It can be shown that this implies that the following 6 conditions are satisfied by the internal coordinates,

These conditions follow from the Eckart conditions that hold for the displacement vectors,

See this article for more details.

References

Cited reference

  1. 1.0 1.1 E. B. Wilson, Jr. Some Mathematical Methods for the Study of Molecular Vibrations, J. Chem. Phys. vol. 9, pp. 76-84 (1941)

Further references

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