ВУЗ:
Составители:
88
H,2,R4,1,A4,3,D4
Variables:
R2=1.3
R3=1.1
R4=0.95
A3=120.
A4=120.1
D4=90.
--Link1--
%Chk=formp
%NoSave
#T RHF/6-31G(d) IRC(RCFC) Guess=Read Geom=AllCheck Test
Для расчета пути следования реакции (зависимость энергии вдоль пути
реакции) используется опция IRC. Такой расчет начинается в седловой точке
и следует вдоль пути в обоих направлениях из переходного состояния к
оптимизации геометрии молекулярной системы в каждой точке вдоль пути.
Окончательный расчет IRC соединяет два минимума на поверхности
потенциальной энергии по
пути, который проходит через переходное
состояние между ними. Процедура для выполнения расчетов IRC
выполняется в следующей последовательности:
• Оптимизация начальной структуры перехода (выполняется один из
расчетов TS, QST2, QST3, Path):
• Запуск расчета частот для оптимизации структура перехода. Этот
расчет необходимо выполнить по нескольким причинам:
• Чтобы убедиться, что седловая точка найдена
правильно;
• Для определения нулевой энергии для перехода структуры;
• Для расчета силовых констант, необходимых для расчета
IRC.
• Выполнение расчета IRC (с ключевыми словами IRC).
В процессе выполнения расчета с целью использования результатов
оптимизации (первый шаг) используются для записи %Chk=FileName, и
Guess=Read, Geom=AllCheck для чтения этих результатов. IRC расчет может
использоваться совместно
с опциями:
Phase=(N1, N2 [,N3 [,N4]]);
Forward следование вдоль пути в прямом направлении;
Reverse следование вдоль пути в обратном направлении;
MaxPoints=N количество рассчитываемых точек вдоль пути реакции. По
умолчанию N=10;
StepSize=N размер шага в единицах 0.01 Бор. По умолчанию 10.
93
0.4573695180D+03 0.1403732281D-01
0.1039486850D+03 0.6884262226D-01
0.2921015530D+02 0.2321844432D+00
0.9286662960D+01 0.4679413484D+00
0.3163926960D+01 0.3623119853D+00
Atom C1 Shell 2 SP 3 bf 2 - 5 -0.087148798088
1.236442079086 0.000000000000
0.7868272350D+01 -0.1193324198D+00 0.6899906659D-01
0.1881288540D+01 -0.1608541517D+00 0.3164239610D+00
0.5442492580D+00 0.1143456438D+01 0.7443082909D+00
Atom C1 Shell 3 SP 1 bf 6 - 9 -0.087148798088
1.236442079086 0.000000000000
0.1687144782D+00 0.1000000000D+01 0.1000000000D+01
Atom C1 Shell 4 SP 1 bf 10 - 13 -0.087148798088
1.236442079086 0.000000000000
0.4380000000D-01 0.1000000000D+01 0.1000000000D+01
Atom C1 Shell 5 D 1 bf 14 - 19 -0.087148798088
1.236442079086 0.000000000000
0.8000000000D+00 0.1000000000D+01
…
Строка
Atom C1 Shell 1 S 6 bf 1 - 1 -0.087148798088 1.236442079086 0.000000000000
означает: Тип атома (Atom C1); Тип оболочки (Shell 1 S); Количество функций
(6); Координаты X, Y, Z (-0.087148798088 1.236442079086 0.000000000000).
Затем следует информация об экспоненциальном коэффициенте и
коэффициенте контракции
0.3047524880D+04 0.1834737132D-02
0.3163926960D+01 0.3623119853D+00
Для 2p оболочки волновая функция будет равна
.
Функция имеет вид, где нормирующий коэффициент.
