Квантово-химические расчеты в программе Gaussian. Серба П.В - 36 стр.

UptoLike

56
-------------
1 8 -0.000068691 0.000000000
-0.000051574
2 1 -0.000574392 0.000000000
0.000836549
3 1 0.000643083 0.000000000
-0.000784975
------------------------------------------------------
-------------
Cartesian Forces: Max 0.000836549 RMS
0.000479219
GradGradGradGradGradGradGradGradGradGradGradGradGradGra
dGradGradGradGrad
Berny optimization.
Internal Forces: Max 0.001030037 RMS
0.000903979
Search for a local minimum.
Завершение второго шага итерационного процесса:
Step number 2 out of a maximum of 20
Третий и последующие шаги итерационного процесса:
All quantities printed in internal units (Hartrees-
Bohrs-Radians)
Update second derivatives …
…......................................................
...........
Max 0.000012189 RMS
0.000010809
Search for a local minimum.
Step number 3 out of a maximum of 20
Последний шаг итерационного процесса:
All quantities printed in internal units (Hartrees-
Bohrs-Radians)
Update second derivatives using D2CorX and points 1
2 3
Trust test= 9.77D-01 RLast= 5.66D-03 DXMaxT set to
3.00D-01
The second derivative matrix:
R1 R2 A1
R1 0.59670
R2 0.01989 0.59670
A1 0.03340 0.03340 0.17798
Eigenvalues --- 0.17296 0.57681 0.62161
65
Следующее задание иллюстрирует модификацию внутренних координат:
# HF/6-31G(d) Opt=ModRedun Test
Opt job
0,1
C1 0.000 0.000 0.000
C2 0.000 0.000 1.505
O3 1.047 0.000 -0.651
H4 -1.000 -0.006 -0.484
H5 -0.735 0.755 1.898
H6 -0.295 -1.024 1.866
O7 1.242 0.364 2.065
H8 1.938 -0.001 1.499
3 8 Добавляет связи атомов H
2 1 3 Добавляет угол C-C-O.
3.2. Анализ заселенностей по Малликену
Вероятность электронного распределения в молекуле определяется
функцией
, причем нормировка требует, чтобы , где
полное число электронов. Для однодетерминантной волновой функции, в
которой молекулярные орбитали представлены в виде линейной комбинации
базисных функций , функция плотности вероятности дается
выражением
, (3.1)
где
-элементы матрицы плотности.
Согласно Mалликену, анализ заселенности можно произвести путем
интегрирования выражения (3.1), что ведет к формуле
, (3.2)
где
-матрица интегралов перекрывания на базисных функциях, которые
являются нормированными, то есть
. Диагональные члены
характеризуют полную заселенность орбитали
(net population). Сумма
недиагональных компонент выражения (3.2), и , равных
по величине, называется заселенностью перекрывания (overlap population)
Заселенность перекрывания связана с двумя базисными функциями
и
, которые могут быть как на одном и том же атоме, так и на двух
различных атомах. Теперь полный электронный заряд можно представить в