An approximate method for calculating mean amplitudes of vibration from spectroscopic data—II

An approximate method for calculating mean amplitudes of vibration from spectroscopic data—II

Spectrochlmicn Act& 1050, pp. 341 to 348. Pergmon Press Ltd. Printed in Northern Irelmd An approximate method for calculating mean amplitudes of...

401KB Sizes 4 Downloads 77 Views

Spectrochlmicn

Act& 1050, pp. 341 to 348. Pergmon

Press Ltd.

Printed

in Northern

Irelmd

An approximate method for calculating mean amplitudes of vibration from spectroscopic data-II Application to simple polyatomic molecules S. J. CYVIN Institute

of Theoretical

Chemistry, Trondheim,

(Received

Technical Norway

21 December

University

of Norway,

1958)

Abstract-The approximate method for calculating mean-square amplitudes of vibration which was described in the first article of this series, has been applied to some special cases. The meansquare amplitudes have been expressed in terms of the atomic maasea and force constants for the following molecular structures: (i) linear eymmetrical XY,, ,(ii) linear symmetrical X,Y,, (iii) tetrahedral X, and (iv) plane symmetrical XY, molecules. IIltrodllction

IN the first article

of this series [l] (hereafter referred to as Paper I) an approximate method for computing the mean-square amplitudes of vibration in polyatomic molecules has been developed. It is an advantage of the method, that in several cases the mean-square amplitudes for the considered molecular model may be expressed explicitly in terms of the atomic masses and the force constants of the molecule. In the present work, this sort of expressions are evaluated for some symmetrical molecules including the plane symmetrical XY, molecular model as the most complicated case. According to equations (16) and (20) of Paper I, the matrices G, GFG, F-1 and. (FG F)-l are to be evaluated. However, the first term of the mean-square amplitude expression is always equal to a(,ux + ,a~), where px and ,uv denote the inverse masses of the adjacent atoms attached to the considered interatomic distance [2]. This fact may serve as a good check if the evaluation is performed by means of the G matrix. The present treatments of the special molecular models are in correspondence with the following scheme: (i) setting up the G and F matrices; (ii) determination of the corresponding inverse matrices; (iii) performing the matrix multiplications GFG and F-lG-lF-l; (’iv ) ex p ressing the interatomic displacements in terms of the chosen internal co-ordinates; and (v) determination of the terms of the meansquare amplitude formulae.

Linear symmetrical XY2 molecules . The simple case of linear symmetrical XY, molecules is taken as the first illustration of the method, although for this case the rigorous method for small harmonic vibrations [3] may also lead to explicit expressions [4]. [l] S. J. CYVIN, Spectrochim. Acta. 66 (1959). [2] E. B. WILSON, JR., J. Chena. Phys. 7, 1047 (1939); 9, 76 (1941); Y. MORINO, K. KUCHI~SU, A. TAKAEASHI and K. MAEDA, J. Chem. Phvs. 21, 1927 (1953). [3] Y. MORINO, K. KUCHITSU and T. SHINIANOUCHI, J. Client. Phys. 20, 726 (1952). [4] S. J. CYVIN, J. Chem. Phys. 30, 337 (1959).

341

S. J. CYVIN

Since we are interested in the interatomic displacements to a first approximation, only the non-degenerate vibrational modes need to be taken into account. The evalustion will be worked out in terms of the valence force co-ordinates rr, (r2, In these co-ordii.e. the deviations from equilibrium for the two X-Y distances. nates, the G and F matrices are reading*

k is the force constant for the X-Y constant. The corresponding inverse

stretching matrices

and k’ the bond-bond sre easily found to be

interaction

1 F-l=iqk +k')(lc - k')--k'/(k + Ic')(k k') k/(k +k')(lc - k') 1 [

G-l

-_

(Px + PYw+ux

PXIPPX

+ PY)PP

+ PY)PY

(Ps + PYWPUS

[

The G FG and F-l G-l F-l matrices GFG=

a

are specified b

[ 1

as follows,

F-lG-lF-1

;

+ PY)PY

=

[ I A

B

A

a

a = Q(pus+ pY)2 + px21 - 2k'CPx + PY)PX b = --BWps + /-+)ps + JWPS + PY)’ + ,d’l A = [(k2 + k'2)(px + PY) - 'L~Qu,ll(~ + k'12(h- W2(?4c + IUP)/+ B = [(k2 + V2)ps

- 21ck’(,ux + ,uJl/(k

+ ~%‘)~(k - 1~‘)~(2/4~ + ,uy),uy

The expression for the mean-square amplitude of the X-Y distance is obtained immediately from the diagonal elements of the evaluated matrices. The Y . . . Y interatomic distance deviations are expressed by TY...Y

=

1'1 +

r2

Hence the matrix Ui which enters into equation (20) of Paper I is the column matrix (1, l} in the present case, and the expressions which enter into the meansquare amplitude formula are consequently of the form 2(z + y). Here z denotes the diagonal and 2~ the off-diagonal element of one of the matrices G, G FG, F-l and (FG F)-I. The resulting formulae for the mean-square amplitudes (~2) are as follows : G&Y u;...y=

= dlux + yk(k 2apy

+ IUY) + Id)-l(k

hmk + W/JY2 + (k - wq-4~ + pyJ21 - l/)-l - $3[(k + ii’)-2/.+-l + (k - Id)-2(2px

- 2/9(k + k’)py2

* Throughout this paper omitted. ,ux and ,uy denote

+ 2y(k

+ k/)-l

- 2d(k + Ic’)-2py-1

the elements beneath the main diagonals the inverse masses of the X and Y atoms,

342

+ py)-l]

of symmetrical respectively.

matrices

will

be

in

approximate

method

for

calculating

mean

amplitudes

of vibration

from

Linear symmetrical X,Y, molecules

spectroscopic

data-II

_

It is profitable to utilize the molecular symmetry when working out the expressions in this case, since it is somewhat more complicated than the previous one. As before, only the non-degenerate vibrations will be taken into account. The inverse kinetic energy matrix and the potential energy matrix in terms of

Fig.

1. Notation used for the linear symmetrical X,Y, molecule. displacements from the equdibrium distances. Symmetry s, = 2-qr, + r*) Syrnm. species A,,: i s, = cl! $Symm. species A,,: s, = 2-q?., - T*)

The symbols co-ordinates:

denote

the

symmetry co-ordinates shall be denoted by theGermanletters 8 and $j,respectively. The applied symmetry co-ordinates are specified by Fig. 1 with the adherent text. In terms of these co-ordinates the energy matrices are reading Px +

PY

--2$is

0

211s

0

6= i

Ps f

k, + k,’ ; 5= PYI

2&k,’

0

k2

0

1

k, - k,’ I

1

The potential energy matrix, has been expressed by means of the force constants attached to the valence force co-ordinates, not the symmetry co-ordinates. k, is the force constant for the X-Y stretching and k’, for the interaction between the two X-Y bond stretchings. k, is the force constant for the X-X stretching and k’, for the X-Y and X-X bond interactions. By inversion one obtains

L

UPS + PY)l

L

KQJ

where

~1 = b/[(k, + k,'P, - 2(kz')21 K2 = (k, + W/W& + k,‘P, - 2(b,‘)21 K’ = -2%‘/[(kl + k,‘P, - 2(k2’121 K2

=

By matrix

l/(k,

- k,‘)

multiplication

@go

f

it is found

b”

81;

8-16-l8-1=

[^’

c

0

B

0 A2

343

1

s. J. h-VIN

a1

=

b

= 2(k, + k,’ + 2k, - 4k2’),ux2

(k,

+

c

= --2Vc,

a2

=

4

= C(kz - k,‘)2ps

B C

= !d(Jc, + h’ - 2bY2px = 2-*{[(k, + k,‘)k, - ,&

A2

=

(k,

h')(px

+

+

-

pd2

k,')(p,

k,')(px

+

+ +

2k2/4s2

py)px

--4k,'(px

-

+

PUY)PX

2%px2

+

%Vps

+ k,‘F,

-

W4)21-2(pxp~)-1

+

-

/-da

+ W)2,411k~

+ (k, + k~‘)2/4Hk~ 4 k,‘)k, - 2(k,‘)21-2(/w~)-1 + k,’ + 2k, - %‘)k,‘lpx - (k, + Wchs} x [(k,

(k,

-

A set of

k,')-2(px

+

+ h’)k,

-

2(k,‘)21-2(luxruy)-1

PI-)-~

representative

four

PY)PX

interatomic

distance

TX-Y

= 2-ybs,

1’x-x

= s,

TX . ..p

=

Let one of the matrices Q, $g%, matrices are of the typical form

is

given

by

+ S,)

2-y31

+

rp.. . p = 2q

deviations

8,)

+

82

f s,

3-1 and (gSi’j)-l

z

0

Y

0

be denoted

by

8.

All these

i 1 Xl

8=

52

For computing the expressions

the mean-square vl@Vl, F,SjV,,

amplitudes according to equation (20) of Paper v&V2 and i&V4 are needed, where

v, = (2-1, 0, 2-t);

v,

= (0,

v, = (2-5, 1, 2-i};

v, = (26, 1, O}

I,

1, O}

It is found = *(xl

V&V2

= $(x1 + x2) + y + 2%; V&V4

The mean-square &-Y

= a(px + SyVMk, -

t.&-=

+x21;

= y = 2x, + y + 2%

amplitudes

may now be evaluated.

+ PUP) - f&h(cLx

+ PY)~ + ~SPX~ -

+ k,‘)k,

%$,‘pxlW, = Zaps -

V2~V2

w4

- W2’)23-1

+ k,‘P,

-

+ (k, - W-l) - 3Wh2px + (k2’)2(clx + PI-) 2(k,‘)21-2(px/+r)-1 + (k, - k,‘r2(px + prF1}

2B(k, + k,’ + 2k2 -

WlXk, + k,’ -

~~‘CPU, + PYJPUXI

4k2’),ux2 + y(k, + k,‘)[(k,

2k,‘)21ux + (k, + k~‘)2puyl[(k~ 344

+ k,‘W,

+ k,‘)k, -

- 2(k2’)2]-1

2(k,‘)21-2(pxpy)-’

An

approximate

&...Y

=

4~s

method

for

+

-

1~1-1

calculating PP~(Ps~

+ h4?(k, + k’,) +

mean +

amplitudes

py2)

of vibration

2k’,~xpu-k-

-

from

spectroscopic

dat+II

+ k,,+i2 - 2k’dps - py)~xl

- 4k’d& + Ic’,P, - 2Wz)21-1 + (k, - k’,Y} - MUk, + k’J2(,e + PY) + kz2/+z+ W, + k’,Vws + W2Ppx + PY) - -Yk, + k’,P’,Plus + PY) - 6k,k’,pxl[@, + k’,)k, - 2W,)21-2(pxpyY + (k, - k’,F2(ps + ruyW k,

= 2wY - Wkl + k’dpY2 + y& + k’, + 2k2 - 4k’,)[(k, + k’,P, - 2W,)21-1 - -kY& + k’A2(px + puy)+ 4bi2px + 4(k, + k’,)k,px + 4(k’A2(4,e + PY) - 4(k, + k’,)k’,(2ps + puy)- l%k’,pu,lW, + k’,)k, - W’,)21-2(/wy)-1

~"Y...Y

Tetrahedral X, molecules The development for this type of molecules turns out to be less complicated than the previous case due to the high molecular symmetry. In fact, the treatment of tetrahedral X, molecules may also be carried through by the rigorous method.

Fig.

2. Notation

symm. symm.

species

used for the tetrahedrnl from equilibrium. A, :

species

E:

symm.species The symmetry be written

P, :

co-ordinates

4 s,, 1&, s,, s,, s,,

X, molecule. Symmetry = 6-+1 = 124(2r, = f(f.2 = 2-qr, = 2-q?, = 2-k(Ta

are given

The symbols denote displacements co-ordinates: + rz + ~3 + ~4 + rs + ~6) - T2 - 9.3 + 25-4 - 7.5 - r,) r3 + rs - To) - +J - TJ - ?-J

by Fig. 2 and the energy

matrices

may

8 = diag (4~~~ psj ius, 2px, 2px, 2~~) 5 = diag (&, g2, K,, K,, K,, K3) Let US assume that the potential field is given by the following force constants: k for the bond stretching, k’ for the interaction between a pair of adjacent bonds and k” for a pair of opposite bonds. Then we have K, = k + 4k’ + k”, K, = k .- 2k’ + k” and K, = k - k”. All the matrices to be evaluated are diagonal and obtained straightforward. 8-l 3-l OgS 5-16-1g-1 3

= diag (&+-l, = diag (K,-l,

,ux-i, ~~-1, +ps-l, +,us-l, &,~~-l) K2-l, K2-l, KS-l, KS-l, KS-l)

= diag (16RlpUX2, R2px2, . . . 4K,pX2, . . , . . ) = diag ( $K,-2px-1, K2-2px-1, . . , -&K3-2ps-1, ..,..) 345

S. J. CYVIN

The X-X

displacements

may be represented ‘)‘z= 6-*s,

yielding

the following

result

t&a = 2u/& - /?(jK, + iK,

+ 3-f&,

by

$ 2-Q5,,

for th’e mean-square

+ 2K,)pg+

amplitude.

y(&-l+

iKz-l

--6(&K,-a

+ iK,l)

+ iK,-2

+ $K,-2)/&

Fig. 3. Notation used for the plane symmetrical XY 3 molecule. The symbols denote the displacements from the equilibrium distances and angles. The equilibrium X-Y &stance is denoted by R. Symmetry co-ordinates: Symm. species A,‘: 4 = 3-)(71 -I- ?a + pa) St,, = 6-4(27, - r2 - T-J S,, = 6-?R(2a, - a, - as) Symm. species E’: = 2-q, - r3) 23i = 2-6R(a, - a$)

Plane symmetrical XY3 molecules In this case only the in-plane vibrations need to be considered. For. the introduction of symmetry co-ordinates, see Fig. 3. In terms of these co-ordinates, the energy matrices may be written

‘PY

0

%3/b

+

2PY)

0

0

0

$3"pu,

0

0

0

0

%3/b

6=

+

2clY)

3(3/b

+

2clY)

:3$% %(3Px f2PY)

L -K,

0

0

K,l? iv=

CD

0

O-

0

0

0

0

K,

r CD

346

--

An

approximate

If the 2f’R2

method

potential 2

aiaj

+

for

c&lculating

fun&ion 2gR

2

is given

riai

i-zj k + 2k’, K, = k - k’, @ = brevity be given by specifying way, s,:

= l//k;

6,; 3,: g,',

= (3~x + 2,+-)/6(3clx = l/R,; = K,/(K,aJ - l?);

By matrix

a2

=

multiplication

+

iK2(3f-4~

b =WbX2 c =

23”@2

mean

+

2g’R

f

amplitudes

of vibration

from

spectroscopic

d&-II

as 2 V = k C ri2 + 2k’ .z.ri~, +fR2 .x1. (riaj

q-,)2

+

z CC>+

(i, j = 1,‘;: 3), one ha; X, =

- f’ izd I’ = g - g’. The inverse matrices will for the different non-vanishing elements in the following

+ PY)PY;

6: @Ii 3,;

= (3k + 2h.)/2(3k = --3%&(3~u, + = Q/(K',@ - I?)

3;;

= --~/(K,Q

+.P~)P~ PUP)PY

- r2)

one obtains

+ YQx2 + 2d2

21~~)~

+8@(3~, + 3@)(3px

+ Wv~x + :3*r(3px

+ 2py)px

+ 2py)px + 2puy)px

+ uwpx

+ 2py)2 + 9/4x2]

A, = K1-2,uy-1

A2 = M3Q2

+ r2)(3px

+ 2py) +~3Wpx](~,o

- ry(3px

+ pu,)-lpp-~

B = E(K22 + 3rwpx + 2~~) + 3*~c2rpxw2~ - r2)-y3pu, C = -WW2Q + r2bx+ WC, + 3wy3px + 2py)] x (K,@ - ry(3px The interatomic

If one of the matrices

distance

deviations

are represented

by

rxmy = 3-*Sl + 6-*2S,, )'y...y =s, - 2-*S2, + S-5!& 8, 636, 3-l and (565)-1 is denoted x1

0 22

xi=

0

0

o-

2

0

0

Y

0

0

x2

2 Y-

347

by

+ E"p)-lpy-l + py)-lpy-l

S. J. Cm

it is found that the expressions +x1 + $x2 and x1 + $x2 + +y - 3-&z are to be respectively. used for computing the X-Y and Y . . . Y mean-square amplitudes, The resulting formulae are given below : U2X-Y

= 4pu~ +

PY)

+

--B&=&~Y~

W2(3p~

+

+ %%)pxl + M&-l - 4wiI-2~~-1 + w3P + rv3~ +

x U"Y...Y

=

2aruY

J72)-2(3Px

-

BUG

+

+ y[lil-l

+*(I<,

+

+

&iW2

+ 3-*W2

WC2

+

+

+

t&x2

21.4~) + 3*.w7~~1

PY)-lPY-l)

34)

+ 36rlpu,2

+ 34 + 3*2r)pi,+

3412clx

+ 34

+

+

+ 24w,+ - r3-11

3’u3pu,

(K24

-5%~)~

9(Jff22

+ 3*r)r(%

+

- ry-11

- 6{1i,-2~~-1

!V2)pu,

+ P~)I(K~+ - JJ2r2(3px + P~)-~P~-~>

Conclusion As a conclusion it may be stated that the general formulae turn out to be rather complicated when the considered molecular models are not quite simple. But individual molecules of higher complexity may be treated ?tumericaEly by the described method.

348