Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

QuantMech_Chem_p3_85

.pdf
Скачиваний:
150
Добавлен:
11.06.2015
Размер:
1.43 Mб
Скачать

Поэтому для основного (невырожденного) состояния эффект Штарка первого порядка отсутствует. Рассмотрим первое возбужденное (n = 2) состояние. Оно четырехкратно вырожденно. Нормированные волновые функции для разных значений орбитального момента | 2lm > имеют вид:

Ψ00(r,r1) :=

1

 

 

2

r

exp

 

r

i := −1

 

 

 

 

 

 

 

4 r1 2 π r1

 

r1

 

2 r1

 

 

 

 

 

 

Ψ10

(

r,r1

)

:=

 

 

 

 

 

1

 

 

 

r exp

 

 

r

 

 

 

(

θ

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2 r1

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 r1

 

 

2 π r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ11

(

r,r1

,θ,φ

)

:=

 

 

 

 

 

1

 

 

r exp

 

 

r

 

 

(

θ

)

exp

(

i φ

)

 

 

 

 

 

 

 

2

 

 

 

 

 

2 r1

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 r1

 

π r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ1_1

(

r,r1,θ,φ

)

:=

 

 

 

 

 

1

 

 

r exp

 

r

 

 

(

 

)

exp

(

i φ

)

 

 

 

 

 

 

2

 

 

 

2 r1

sin

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 r1

π r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оператор возмущения V является нечетным относительно отражения в начале координат. Поэтому отличными от нуля будут только матричные элементы для невозмущенных состояний с противоположными четно-

стями.

Следовательно,

единственный ненулевой матричный элемент

< 200 | V | 210> , в чем можно убедиться и непосредственным вычислени-

ем интегралов.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

)

 

π

 

 

 

(

 

 

)

 

(

)

 

2

(

)

 

U01

r1

,qε

:=

 

 

 

 

 

 

 

 

 

dθ dr

 

 

Ψ00(r,r1) V qε,r

 

Ψ10 r,r1

 

2π r sin θ

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U01(r1,qε) assume,r1 > 0

3 r1 qε

 

 

 

 

 

Согласно [1-3], при наличии вырождения поправки первого порядка к

энергии находятся как корни векового уравнения

 

 

 

 

 

 

 

 

 

 

E 3 r1 qε 0

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

0

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

3 r1 qε −

 

 

solve , E

 

 

 

 

 

 

 

 

 

 

0

 

0

 

E 0

 

 

3 r1 qε

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

E

 

 

 

 

 

3 r1 qε

 

 

Пусть Eo - энергия исходного вырожденного состояния. Тогда с учетом поправок исходный уровень с n = 2 расщепляется на три разных уровня. Поскольку E пропорциональны ε , то это - линейный эффект Штарка.

E1 = Eo 3

q

ε r1 E2 = Eo + 3

q

ε r1

E3 = E4 = Eo

61

Упражнения

1.Вычислить поправки первого и второго порядка по теории возмущений для первого возбужденного состояния осциллятора с потенциалом V(x) = a |x| . Оценить изменение частоты перехода (n=1--> n=0) для принятых выше параметров молекулы водорода и потенциала возмущения.

2.Для ангармонического осциллятора с потенциалом возмущения вида

V(α,β,x) = α x3 + β x4

проверить, что поправка первого и второго порядка к энергии n - го состояния гармонического осциллятора равна (см.[3-5])

E1(β) =

3 h2 β

n2

+ n +

1

E2(α) = 15 h2 α2 n2

+ n +

11

 

 

2 μ

2

ω

2

 

 

2

4 μ

3

ω

4

 

 

30

 

 

 

 

 

 

 

 

 

 

 

3. Показать, что в задаче об атоме водорода в электрическом поле

уровню

E1

соответствует

волновая функция (| 200 > + | 210 >) / 2 , а

уровню

E2

соответствует

функция (| 200 > - | 210 >) / 2 . Обе эти

функции являются гибридными (диагональными) sp - орбиталями.

4. Пространственный ротатор с моментом инерции J и дипольным моментом d, направленным вдоль оси ротатора, помещен в однородное электрическое поле ε, рассматриваемое как возмущение. Найти первые неисчезающие поправки к энергии основного уровня и показать (см[4]),

что поляризуемость основного состояния ротатора равна β = 2Jd2/3h2 .

5.В первом порядке теории возмущений найти поправку к уровням энергии водородоподобного иона при учете слагаемого -- p 4/8m3c2 в квазирелятивистской функции Гамильтона электрона в поле ядра (см[4]).

6.Cчитая, что гамильтониан спин-орбитального взаимодействия имеет

вид

 

h2

 

 

1

d

 

(L S)

,

 

V(r) =

 

 

 

 

r

 

U(r)

 

2

 

2

 

 

c

 

dr

 

 

 

 

 

2m

 

 

 

 

 

 

 

найти в первом порядке по теории возмущений энергию спин-орбиталь- ного взаимодействия Eсо = < nlm | V(r) | nlm > для водородоподобного иона. Вычислить эту энергию для электрона, находящегося на 2p-AO атома водорода (см.[4]).

7. Найти энергию основного состояния атома гелия в первом порядке теории возмущений, рассматривая взаимодействие электронов как малое возмущение. Сопоставить результат с вариационным расчетом и указать возможные причины различий (см.[3,4]).

62

Элементы теории групп и представлений групп Задача 11

Для правильного треугольника: а) указать все элементы симметрии и построить взаимно однозначное соответствие с элементоми группы перестановок ядер симметричной трехатомной молекулы; б) построить совокупность матриц третьего порядка, однозначно связанных с элементами группы перестановок и составить таблицу умножения этой группы; в) найти неприводимые представления и их характеры; г) построить базисные функции неприводимых представлений; д) найти матрицы прямого произведения представлений.

Решение

Элементы симметрии и перестановки

Обозначим вершины треугольника цифрами 1,2,3, а соответствующие оси симметрии буквами σ1,σ2,σ3 (Рис.1). Имеем 6 операций симметрии:

Е - единичная операция; σ1 - отражение в плоскости σ1;

σ2 - отражение в плоскости σ2; σ3 - отражение в плоскости σ3;

C3 - поворот вокруг оси, проходящей

 

через центр треугольника пер-

 

пендикулярно плоскости по

 

часовой стрелке на угол 1200;

 

С32 - поворот вокруг этой оси на 2400.

 

Легко убедиться в том, эти элементы

Рис. 1. Элементы симметрии

составляют группу C3v шестого

порядка (см. [7-9,14]).

Выполняя операции симметрии, обнаружим, что каждому элементу можно поставить в соответствие определенную перестановку индексов

1, 2 и 3, а именно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

1 2 3

 

С

3

1 2

3

 

С 2

1 2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2 3

 

 

2 3

1

 

3

 

 

3 1

2

 

 

 

 

 

 

 

 

 

 

σ

1

1 2 3

 

σ

2

1

2

3

 

σ

3

1

2

3

 

1 3 2

 

3

2

1

 

2

1

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученная группа перестановок обозначается Р3. Соответствие между элементами групп однозначно и не нарушается при выполнении опера-

ций умножения перестановок (справа налево).

Например,

 

σ3 = С3 σ1

1

2

3

1

2

3

=

1

2

3

↔ σ3

2

3

 

1

3

2

2

1

3

 

 

1

 

 

63

Матрицы и таблица умножения группы S3

Сопоставим каждой вершине {1, 2, 3} треугольника вектор по правилу:

 

1

 

 

0

 

 

0

 

X :=

0

 

X :=

1

 

X :=

0

 

1

 

 

2

 

 

3

 

 

 

0

 

 

0

 

 

1

 

Тогда каждая операция симметрии может быть представлена матрицей. Это представление является точным. Вид матриц легко получить, анализируя их результаты действия на векторы X1, X2 и X3 и сопоставляя с операциями симметрии.

 

1 0

0

 

0 0 1

 

0 1 0

E :=

0

1

0

 

A :=

 

1 0 0

 

B :=

0 0 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

1

 

 

0 1

0

 

 

1 0 0

 

 

1 0

0

 

0 0

1

 

0 1

0

C :=

0

0

1

 

D :=

 

0

1

0

 

F :=

1

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

0

 

 

1

0

0

 

 

0

0

1

 

Группа не является абелевой, поскольку произведение операций сим-

метрии некоммутативно. Например,

 

 

 

 

0

0

1

 

 

0

1

0

 

C A = 0

1

0

= D

, но

A C = 1

0

0

= F

 

 

 

 

 

 

 

 

 

1

0

0

 

 

0

0

1

 

Данная группа имеет следующие подгруппы: 1) тривиальная подгруппа, состоящая из единицы Е; 2) подгруппа, состоящая из элементов {E,A,B}; 3) подгруппа {E,C}; 4) подгруппа {E,D}; 5) подгруппа {E,F}. Порядки подгрупп равны 1, 3 и 2. Это соответствует теореме Лагранжа, согласно которой, порядок подгруппы конечной группы является делителем порядка группы.

Для нахождения таблицы умножения группы (квадрата Кэли) воспользуемся описанной ниже программной функцией MTable. Для ее реализации представим элементы группы симметрии как элементы одномерного символьного массива t и как элементы трехмерного массива вложенных (nested) матриц G. Функция вычисляет произведение матриц и сопоставляет его с каждым элементом группы S3 .

t0 := "E"

t1 := "A"

t2 := "B"

t3 := "C"

t4 := "D"

t5 := "F"

G0 := E

G1 := A

G2 := B

G3 := C

G4 := D

G5 := F

64

MTable(g) := for i 0 .. 5

T0,0 " S3 "

T0,i+1

concat("_" ,ti,"_" )

Ti+1,0

concat("_" ,ti,"_|" )

for j 0 .. 5

 

 

 

for

k 0 .. 5

 

 

 

 

 

Ti+1, j+1 tk if gi g j = gk

 

 

 

continue

T

{3,3}{3,3}

{3,3} G = {3,3}{3,3}

{3,3}

 

" S3 " "_E_" "_A_" "_B_" "_C_" "_D_" "_F_"

 

"_E_|"

"E"

"A"

"B"

"C"

"D"

"F"

 

 

 

 

"_A_|"

"A"

"B"

"E"

"F"

"C"

"D"

 

MTable(G) =

"_B_|"

"B"

"E"

"A"

"D"

"F"

"C"

 

 

"_C_|"

"C"

"D"

"F"

"E"

"A"

"B"

 

 

 

 

"_D_|"

"D"

"F"

"C"

"B"

"E"

"A"

 

 

"_F_|"

"F"

"C"

"D"

"A"

"B"

"E"

 

 

 

Матрицы поворотов и отражений

Построим еще одно представление группы C3v . На этот раз - двумерное как преобразование координат точки (x,y) на плоскости вокруг центра треугольника, оставляя z- координату неизменной и принимая равным единице расстояние от центра треугольника до вершины.

Такое представление состоит из матриц (2x2), вложенных в трехмерные матрицы. Их можно получать на основе матрицы преобразования координат вектора при повороте вокруг центра треугольника в плоскости, перпендикулярной оси Z, и матрицы отражений. Матрицы собственных и несобственных (с отражением в горизонтальной плоскости) вращений вектора по часовой стрелке показаны ниже (см. [7-9,14]).

 

 

 

cos(φ)

sin(φ)

0

(φ) :=

 

 

 

(

φ

)

(

)

 

R1

 

 

 

sin

 

cos φ

 

0

 

(

 

 

(

 

0

 

 

0

 

1

det

 

 

φ

))

simplify

1

 

R1

 

 

 

cos(φ)

sin(φ)

0

 

R_1

(φ) :=

 

(

φ

)

(

)

0

 

 

 

sin

 

cos φ

 

 

 

 

 

0

 

 

0

 

1

 

det(R_1(φ)) simplify

→ −1

65

Поворотам на 1200 и 2400 по часовой стрелке отвечают матрицы a и b. Отражения в плоскостях σ1 , σ2 и σ3 реализуются матрицей с и как комбинации d = c a и f = a c .

 

1 0

0

 

 

 

1

3

0

 

 

 

 

1 3 0

e :=

 

0

1

0

 

 

a :=

1

3

1 0

 

b :=

1

3

1 0

 

2

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

0

 

0 2

 

 

 

 

0

0 2

 

1 0

0

 

1

1 3 0

 

 

1

1

3 0

c :=

 

0

 

1

0

 

d :=

3 1 0

 

f :=

 

3

1 0

 

2

2

 

 

0

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

1

 

 

0

 

2

 

 

 

0

0 2

Ортогональность и характеры представлений

Следуя законам преобразований функции при действии на нее операторов симметрии, можно получить два одномерных представления группы C3v . Поскольку

Γ(E) ψ = E ψ

и Γ (σ1)2 ψ = E ψ ,

то Γ(E) = 1 и Γ(σ1) = ± 1

Аналогично,

Γ(σ2) = ± 1

и

Γ(σ3) = ± 1

Поскольку

C3 = σ1 σ2 , то

 

Γ(C3)= Γ(σ1) Γ(σ2) = 1

 

(C3)2 = C3 C3

и

Γ (C3)

2 = Γ(C3)Γ(C3)= 1

 

 

 

 

 

Поэтому одномерные представления имеют вид :

 

Γ1 := ( 1 1 1 1 1 1 )T

Γ2 := ( 1 1 1 1 1 1 )T

Первое из них является полносимметричным представлением A1, а в представлении Γ2 при отражениях функция должна изменить знак (A2).

В качестве третьего представления выберем двумерное представление Γ3 типа Е, полученное выше с помошью матриц вращения и отражения.

 

 

 

1

0

 

 

 

 

 

1

 

1

 

3

 

 

 

 

1

 

1 3

Γ

 

:=

0

1

 

Γ

 

:=

 

 

 

 

 

Γ

 

:=

 

 

 

 

 

 

 

2

3

 

 

2

3

 

30

 

 

31

 

 

1

 

32

 

 

1

 

 

 

1 0

 

 

 

 

 

1

 

1

3

 

 

 

 

1

 

1

3

Γ3

:=

 

 

 

Γ3

 

:=

 

 

 

 

 

Γ3

 

:=

 

 

 

 

 

4

2

3

 

5

2

3

 

 

3

0 1

 

 

 

 

 

1

 

 

 

 

1

Четвертое представление G - совокупность матриц E, A, ... F группы S 3.

66

Проверим ортогональность матричных элементов, например, для представления Γ3, выполняя суммирование по всем элементам группы. Условие ортогональности имеет вид (см. [7-9,14]):

Γ

α

(g1) Γ

β

(g)

jq

=

g

δ

αβ

δ

ij

δ

pq

 

 

ip

 

 

nα

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для унитарных вещественных матриц обратная равна транспонированной, поэтому, например, для подобных элементов условие ортогональности дает:

5

 

 

 

 

 

5

 

 

 

 

 

Γ

 

Γ

 

= 3

Γ

 

Γ

 

= 3

 

3g 0,0

 

3g

0,0

 

3g 0,1

 

3g

0,1

g = 0

 

 

 

 

 

g = 0

 

 

 

 

 

Эти значения равны порядку группы, деленному на размерность представления 6/2 = 3. Другие элементы представления Γ3 ортогональны, например,

5

 

 

 

 

 

5

 

 

 

 

 

Γ

 

Γ

 

= 0

Γ

 

Γ

 

= 0

 

3g 0,0

 

3g

1,0

 

3g 0,0

 

3g

1,1

g = 0

 

 

 

 

 

g = 0

 

 

 

 

 

В соответствии с определением, вычислим характеры представления как сумму диагональных элементов матриц, например, для представления G.

Char(g) :=

 

c1,0 concat("Ch" )

 

 

 

 

 

 

for

k 1 .. 6

 

 

 

 

 

 

 

 

 

 

c0

,k concat("_" ,tk1,"_" )

 

 

 

 

 

 

 

 

 

 

 

 

 

c1

,k tr(gk1) if

rows(gk1) > 1

 

 

 

 

 

 

 

c1

,k gk1

otherwise

 

 

 

 

 

c

 

 

 

 

 

 

 

Char(G) =

 

0

 

 

 

 

"_E_" "_A_" "_B_" "_C_" "_D_" "_F_"

 

"Ch"

 

3

0

0

1

1

 

 

 

 

1

Условие ортогональности вещественных характеров имеет вид ([9,14]):

χα(g) χβ(g) = g δαβ g

67

Проверим это условие для представления G.

5

 

χG := submatrix(Char(G) ,1 ,1 ,1 ,6)T

χGi χGi = 12

 

i = 0

Поскольку сумма квадратов характеров не равна порядку группы g = 6, то это представление является приводимым. Вычислим коэффициенты его разложения по характерам неприводимых представлений Γ1, Γ2, Γ3.

 

 

χ1

:= Γ1

 

χ2 := Γ2

 

 

 

χ3

:= submatrix(Char(Γ3),1 ,1 ,1 ,6)T

 

m1 := 61

5

 

 

m2 := 61

5

 

m3 := 61

5

χGi χ1i

χGi χ2i

χGi χ3i

 

i = 0

 

 

 

i = 0

 

 

i = 0

m1 = 1

 

 

m2 = 0

 

m3 = 1

Следовательно, представление G разлагается в сумму Γ1, Γ2, Γ3 с коэффициентами m1 = 1, m2 = 0 и m3 = 1.

G = m1 A1 + m2 A2 + m3 E

Базисные функции представления

Пусть задан набор из трех одинаковых по форме функций φ1(r1) , φ2(r2), φ3(r3), значения которых зависят только от расстояния точек до вершин треугольника. Такими функциями могут быть атомные орбитали, центрированные на атомах 1, 2 и 3. Операции симметрии группы С3v переводят эти функции друг в друга. Если представить эти три функции как элементы вектор-столбца, то закон преобразования можно определить следующим образом:

 

φ1

 

 

φ1

 

 

φ1

 

 

φ3

 

 

φ1

 

 

φ2

 

 

 

 

 

 

 

 

 

 

 

E

φ

2

 

φ

2

 

A

φ

2

 

φ

1

 

B

φ

2

 

φ

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ3

 

φ3

 

φ3

 

φ2

 

φ3

 

φ1

 

φ1

 

 

φ1

 

 

φ1

 

 

φ3

 

 

φ1

 

 

φ2

 

 

 

 

 

 

 

 

 

 

 

 

 

C

φ

2

 

φ

3

 

D

φ

2

 

φ

2

 

F

φ

2

 

φ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ3

 

φ2

 

φ3

 

φ1

 

φ3

 

φ3

68

Результаты преобразований R(g) выбранных функций φk и характеры неприводимых представлений группы C3v можно свести в таблицу.

Для конструирования базисных функций, преобразующихся по конкретным неприводимым представлениям группы воспользуемся аппаратом проекционных операторов. По определению проекционных опе-

раторов (см. [7-9,14,15])

 

nα

 

или

 

nα

 

Pαi,k = g Γα(g)k,i R(g)

Pα = g

χα(g) R(g)

 

 

 

g

 

 

 

g

где nα - размерность представления, g - порядок группы, R(g) - операция симметрии, Γα(g)k,i - матричный элемент транспонированной (ве-

щественной) матрицы представления Γα(g), χα(g) - ее характер, а сум-

мирование ведется по всем элементам группы. Для полносимметричного представления А1 оператор и результат проектирования имеют вид:

PA1(Φ) :=

1

χ1

E + χ1

A + χ1

 

B + χ1

C + χ1

D + χ1

F Φ

 

6

0

 

 

1

 

 

2

 

 

 

3

 

 

4

5

 

 

 

 

 

 

1

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3

φ1 +

 

 

φ2

+

 

φ3

 

 

 

 

 

φ

1

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

φ2

1

φ1 +

φ2

+

φ3

 

 

 

 

PA1

3

 

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ3

1

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

φ1 +

 

 

φ2

+

 

φ3

 

 

 

 

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

Предполагая ортонормированность функций φk можно получить нормированную полносимметричную функцию.

ψ

A1

= 1 (φ1 + φ2 + φ3)

 

3

69

Для представления А2 оператор проектирования имеет вид:

PA2(Φ) := 16 χ20 E + χ21 A + χ22 B + χ23 C + χ24 D + χ25 F Φ,

а результат его действия на пробные функции свидетельствует о том, что среди выбранных функций нет базисных функций представления А2 .

 

φ

1

 

 

0

 

 

 

 

 

 

 

 

 

 

PA2 φ2

 

0

 

 

φ3

 

 

0

 

 

 

 

 

 

Для двумерного представления Е нужно построить две ортогональные компоненты. Для этого построим проектор и подействуем им на выбранные пробные функции.

(

Φ

)

:=

2

 

χ3

E + χ3

A + χ3

B + χ3

C

+ χ3

 

D + χ3

F Φ

 

 

PE

 

 

 

 

 

 

 

6

 

0

 

 

 

1

 

2

 

 

2

1

 

 

3

1

 

 

4

5

 

 

 

 

 

 

 

 

 

 

 

 

 

φ1

φ2

φ3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ

 

 

 

3

3

3

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

φ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PE

 

 

 

 

φ1 +

 

 

φ2

 

φ3

 

 

 

 

 

 

 

 

3

 

3

3

 

 

 

 

 

 

 

 

 

φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

1

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

φ1

φ2 +

φ3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После нормировки первой и второй строки находим:

ψE

1

= 1 (2 φ1 − φ2 − φ3)

ψ2

= 1 (2 φ2 − φ1 − φ3)

 

6

 

6

Поскольку первая и вторая строка не ортогональны, то для нахождения второй базисной функции воспользуемся процедурой ортогонализации Шмидта. В качестве ψE2 выберем функцию

ψE2 = ψ2 + CE ψE1 ,

где константа CE должна быть выбрана из условия ортогональности

ψE1 ψE2 dτ = 0

Это условие выполняется, если

70

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]