As calculated, the overlap integrals represent the overlap of atomic orbitals that are aligned along the z-axis. In general, this will not be the case, and the diatomic overlap integral matrix must be rotated in order to represent the actual orientation used. The rotation operation can be constructed in a completely general way, or it can be constructed using trigonometric functions, as shown here.
The rotation matrices are well known, but the method by which they are constructed is by no means simple. Consider the general case of one atom at the origin, and a second atom at some point p(x,y,z). Let the angle from the second atom to the first atom to
the z z axis be b, and the angle of the projection of the second atom onto the xy plane to the y axis be a. Then, the p-orbital
transform is as shown in the following table:
Angular Dependence of the p-orbitals | |||
p(x) | p(y) | p(z) | |
π+ | cos a cos b | -sin a | cos a sin b |
π- | sin a cos b | cos a | sin a sin b |
σ | -sin b | 0 | cos b |
Mathematical Tools for use with Spherical Harmonics |
|
When working with the angular components, it is normally easier to use the Cartesian symbols x, y, and z instead of the trigonometric forms. This avoids the potentially confusing expressions that involve the angles a, b, θ, and φ.
Consider the rotation of d(x2-y2). The rotation components involved are: R(x) = (cosa.cosb)x + (sina.cosb)y + (-sinb)z, and R(y) = (-sina)x + (cosa)y .
R|d(x2-y2)> | = | (cos2a.cos2b)x2 + (2cosa.cosb.sina.cosb)xy + |
(sin2a.cos2b)y2 + (-2sina.cosb.sinb)yz + (sin2b)z2 + | ||
(-sin2a)x2 + (2sina.cosa)xy + (-cos2a)y2 |
In this expression, the integral over all odd terms vanish, therefore:
<d(x2-y2)|R|d(x2-y2)> | = | <x2(cos2a.cos2b - sin2a)x2> + | |
<x2(sin2a.cos2b - cos2a)y2> + | |||
<x2(sin2b)z2> + | |||
<y2(-cos2a.cos2b + sin2a)x2> + | |||
<y2(-sin2a.cos2b + cos2a)y2> + | |||
<y2(-sin2b)z2> |
cos2a.cos2b - sin2a -sin2a.cos2b + cos2a
or(cos2b + 1).(2cos2a - 1)
sin2a.cos2b - cos2a + sin2b - cos2a.cos2b + sin2a - sin2b
or-(cos2b + 1).(2cos2a - 1)
The rotation matrix element is thus:<d(x2 - y2)|R|d(x2 - y2)> = <x4 -2x2y2 +y4>½(cos2b + 1).(2cos2a - 1)
As a result of the fact that the angular components of the atomic orbitals are normalized, i.e., <x4 - 2x2y2 + y4>=1:<d(x2 - y2)|R|d(x2 - y2)> = (2cos2a - 1)cos2b + ½(2cos2a - 1)sin2b
which is the form used in the calculation of the d-rotation matrix. The full set of d-orbital rotation matrix elements is:
Angular Dependence of the d-orbitals |
|||||
|
δ+ ≡ d(x2 - y2) |
π+ ≡ d(xz) |
σ ≡ d(2z2-x2-y2) |
π- ≡ d(yz) |
δ- ≡ d(xy) |
d(x2-y2) |
(2cos2a - 1)cos2b + ½(2cos2a - 1)sin2b |
-cosa.sinb.cosb |
√¾.sin2b |
-sina.sinb.cosb |
2.sina.cosa.cos2b + sina.cosa.sin2b |
d(xz) |
(2cos2a - 1)sinb.cosb |
cosa(2cos2b - 1) |
-√3.sinb.cosb |
sina (2cos2b - 1) |
2.sina.cosa.sinb.cosb |
d(2z2-x2-y2) |
√¾.(2cos2a - 1)sin2b |
√3.cosa.sinb.cosb |
cos2b - ½sin2b |
√3.sina.sinb.cosb |
√3.sina.cosa.sin2b |
d(yz) |
-2sina.cosa.sinb |
-sina.cosb |
0 |
cosa.cosb |
(2cos2a - 1)sinb |
d(xy) |
-2sina.cosa.cosb |
sina.sinb |
0 |
-cosa.sinb |
(2cos2a - 1)cosb |