Что программа умеет считать 6 страница

Эта часть диалогового меню служит для задания требуемой точности расчета SCF волновой функции и максимального количества итераций для достижения этой точности.

• Convergence (Limit Предел сходимости) служит для того, чтобы остановить SCF процедуру, когда разница в энергиях между двумя последующими итерациями становится меньше заданной величины. Для ab initio расчетов это значение, как правило, выбирается равным 0.00001 ккал/моль, а используемый параметр может лежать в интервале от 1 до 0.00000001 ккал/моль. Тем не менее, если этот параметр задать меньше 10**-10, то сходимости можно не достигнуть, так как сам метод Хартри-Фока дает соразмерную ошибку. При поиске переходных состояний рекомендуется использовать более жесткие критерии сходимости.

• Iteration Limit (Предельное количество итераций) определяет максимальное количество итераций в SCF расчете. Расчет останавливается, в случае если программа выполнила заданное количество итераций и при этом не достигла сходимости. В этом случае результат расчета может быть не верен, так как энергия может быть далека от истинной, либо она осциллировала в ходе расчета. 50 итераций представляется разумным для большинства случаев, тогда как задание большего количества (скажем, до 200) может быть оправданным при расчетах переходных состояний. В случае, когда расчет не может сойтись и при большем количестве итераций, то это, как правило, означает, что дальнейшее увеличение количества итераций не может изменить, так как система не сходится. Лишь в некоторых случаях Ц Accelerate Convergence (Алгоритм ускорения сходимости) может исправить ситуацию.

Выбор этого параметра убыстряет сходимость SCF расчетов. При этом HyperChem включает процедуру, известную как «Прямое инвертирование подпространства итераций» (Direct Inversion of Iterative Subspase, DIIS) (Более подробно смотри HyperChem Computational Chemistry), которая может потребовать большего объема памяти. Включение этой опции может увеличить время, затрачиваемое на одну итерацию из-за того, что в этом подходе матрица Фока вычисляется как линейная комбинация текущей матрицы Фока и матриц Фока от предыдущих итераций. Как правило, такой подход уменьшает общее количество требуемых итераций.

Spin pairing (Спиновое состояние).Возможно выбрать два метода расчета спиновых состояний молекул. Первый – неограниченный метод Хартри-Фока (Unrestricted Hartree-Fock method, UHF) и ограниченный метод Хартри-Фока (Restricted Hartree-Fock method, RHF).

•UHF рассматривает спин-орбитали с различным пространственным распределением для a - и b- орбиталей. Этот метод применяется при изучении систем, как с открытыми, так и с закрытыми электронными оболочками. Так, для последних он хорошо описывает реакции диссоциации. Однако, из-за удвоения количества орбиталей, время расчета этим методом увеличивается вдвое. У этого метода существуют и другие ограничения, связанные с его основами.

•В RHF считается, что электроны с различным спином занимают одинаковые, в смысле пространственного распределения, орбитали. При этом неспаренные электроны тоже могут занимать отдельные орбитали. Этот метод применяется как для открытых, так и для закрытых электронных оболочек.

Configuration Interaction (Конфигурационное взаимодействие). Эта опция используется для активации расчета конфигурационных взаимодействий и открывает соответствующее диалоговое окно. Такой подход необходимо применять при расчетах УФ и оптических спектров в видимом диапазоне. Выбор этой опции существенно увеличивает время расчетов.

Кнопка OK служит для сохранения выбранных параметров и закрытия диалогового окна. Кнопка Cancel приводит к закрытию диалогового окна без сохранения выбранных параметров. Как уже говорилось ранее, в программе HyperChem существует возможность определять различные базисные наборы для разных частей рассчитываемой системы.

Необходимо отметить, что HyperChemне поддерживает ограниченный метод Хартри-Фока для систем с открытыми оболочками (ROHF) на ab initio уровне.

Диалоговое окно параметров высокого уровня метода ab initio (Ab Initio Advanced options)<o:p></o:p>

Это диалоговое окно служит для более точной настройки параметров ab initio расчетов.

Integral Format

• Regular определяет использование обычного формата для записи двухэлектронных интегралов. HyperChem использует 16 байт для записи каждого из интегралов. Первые 8 байт хранят 4 индекса интеграла, а последние 4 – его значение. Программа сохраняет эти величины только в том случае, если абсолютное значение интеграла больше или равно параметру Cutoff. В противоположном случае значение этого интеграла приравнивается к 0. Двухэлектронные интегралы и их индексы сохраняются на диске без модификации при выборе опции Regular и могут быть записаны в .log-файл при соответствующем выборе параметра QuantumPrintLevel меню StartLog.

• Raffenetti определяет использование формата Раффенети (R.C. Raffenetti, Chem.Phys.Lett., 20, 335 (1973)), который позволяет более просто формировать матрицу Фока в ходе SCF расчета. Этот формат, как правило, требует больше памяти и больше дискового пространства, однако позволяет повышать скорость расчета. Этот формат нельзя использовать при проведении MP2-расчетов.

Параметр Cutoff позволяет сохранять на диске только те интегралы, абсолютное значение которых равно или превышает задаваемый параметр. По умолчанию он равен 10**-10 Хартри. Этот параметр контролирует осуществление SCF-итераций, точность волновых функций и энергии, так как он уменьшает количество рассчитываемых двухэлектронных интегралов.

<![if !supportEmptyParas]><![endif]> <o:p></o:p>Параметр Buffer size (Размер буфера) определяет размер операционной памяти (в словах двойной точности, 8 байт для одного слова), которая требуется для хранения двухэлектронных интегралов до того, как записать их во временный на жестком диске (выбор этого диска может быть сделан в меню File/Preferences/Path). Более большой размер буфера способен уменьшить расчетное время из-за того, что программа будет обращаться к диску более редко. Если этот буфер будет достаточно велик, то HyperChem не будет обращаться к диску вовсе. Необходимо отметить, что в случае, когда программа останавливается из-за сбоя компьютера, либо по другой некорректной причине, этот временный файл остается на диске. Для того чтобы освободить это место, его необходимо удалять вручную.

Direct SCF calculation вычисляет двухэлектронные интегралы на каждой итерации, а не один раз перед SCF шагом, как если бы это было без нее. Такой расчет, безусловно, гораздо медленнее, но он позволяет не использовать дисковое пространство и операционную память под большое количество интегралов. Включение этой опции требуется при расчетах больших систем на компьютерах с маленьким диском и памятью.

Ghost-atoms Control (Использование атомов-призраков). Эта опция позволяет вводить центры, которым приписываются базисные функции тех или иных атомов, при этом в систему не вводятся ни дополнительные ядра, ни дополнительные электроны. Эта опция позволяет вводить только дополнительные базисные функции, центрированные в любой точке пространства. Необходимо помнить, что использование таких «виртуальных» атомов может приводить к неким артефактам. Например, так как существуют базисные волновые функции, будут существовать и мулликеновские заряды, которые будет центрированы на соответствующих центрах. Для того чтобы использовать эту возможность, необходимо в меню Select, предварительно сделав соответствующее выделение, выбрать пункт Name selection и, выбрав пункт Other, приписать ему имя ghost-atoms. После этого в диалоговом меню Ab Initio Method / Advanced Options становится активным пункт Ghost-atoms Control. Такой подход бывает эффективным при описании ряда неординарных химических связей. Эта опция активна только для расчетов одной точки.

MO initial guess (Параметр стартового заселения МО) определяет стартовое заполнение коэффициентов молекулярных орбиталей при помощи диагонализации остового гамильтониана. При выборе параметра Projected Huckel, эти параметры определяются по методу Хюккеля. Аналогично определяются и первоначальные коэффициенты при выборе параметров Projected CNDO(методом CNDO) и Projected INDO(методом INDO).

Number of d Orbitals (Количество d-орбиталей). Этот параметр определяет вид d-орбиталей, используемых в расчете. Выбор пяти(five) орбиталей соответствует расчету с использованием эрмитовых орбиталей (d 0, d 1, d –1, d 2, d –2), а выбор шести(six) – соответствует расчету с использованием d-орбиталей в декартовом представлении (d _xx_, d _yy, d _zz, d _xy, d _xz, d _yz). <![if !supportEmptyParas]><![endif]> <o:p></o:p>

<![if !supportEmptyParas]><![endif]> <o:p></o:p>С другими особенностями программы HyperChem читатель может познакомиться, изучив файлы Referenc.pdf, CDK.pdf и GetStart.pdf.

Меню Compute определяет вид расчета. Это может быть расчет в одной точке (Single Point), оптимизация геометрии (Geometry Optimization), молекулярная динамика (Molecular Dynamics), ланжевеновская динамика (Langevin Dynamics), расчеты методом Моне-Карло, а также некоторые другие виды расчетов. Далее по тексту книги нас в основном будут интересовать результаты расчетов методом молекулярной динамики, которые представлены файлами с динамическим кино. Для того чтобы просмотреть результаты молекулярно-динамического моделирования, необходимо прочитать соответствующий файл (в тексте они будут указываться), а затем, выбрав в меню Compute пункт Molecular Dynamics, отметить Ц Playback и нажать Proceed. Для увеличения скорости проигрывания файлов (когда это необходимо, в зависимости от сложности задачи и мощности компьютера) можно увеличить time steps (по умолчанию стоит 1, при увеличении скорость увеличивается кратно количеству шагов).

<![if !supportEmptyParas]><![endif]> <o:p></o:p>Как уже говорилось выше, подробная инструкция о пользовании программой HyperChem дается в файлах CDK.pdf, GetStart.pdf, Referenc.pdf. Но научиться с ней работать не составляет большого труда, так как программа сконструирована очень удобно и предназначена для широкого круга пользователей с разным уровнем квалификации.

Литература<o:p></o:p>

1. W. J. Hehre, R. F. Stewart, and J. A. Pople, J. Chem. Phys., 51, 2657 (1969); J. B. Collins, P. V. Schleyer, J. S. Binkley, and J. A. Pople, J. Chem. Phys., 64, 5142 (1976); R. F. Stewart, J. Chem. Phys., 52, 431 (1970).<o:p></o:p>

2. J. S. Binkley, J. A. Pople, and W. J. Hehre, J. Am. Chem. Soc., 102, 939 (1980); M. S. Gordon, J. S. Bin-kley, J. A. Pople, W. J. Pietro, and W. J. Hehre, J. Am. Chem. Soc., 104, 2797 (1982); W. J. Pietro, M. M. Francl, W. J. Hehre, D. J. Defrees, J. A. Pople, and J. S Binkley, J. Am. Chem. Soc., 104, 5039 (1982).<o:p></o:p>

3. W. J. Hehre, R. Ditchfield, and J. A. Pople, J. Chem. Phys., 56, 2257 (1972); P. C. Hariharan and J. A. Pople, Theor. Chim. Acta, 28, 213 (1973); M. S. Gordon, Chem. Phys. Lett., 76, 163 (1980).<o:p></o:p>

4. R. Krishnan, J. S. Kinkley, R. Seeger, and J. A. Pople, J. Chem. Phys., 72, 650 (1980); A. D. McLean and G. S. Chandler, J. Chem. Phys., 72, 5639 (1980).<o:p></o:p>

T.H. Dunning and P.J. Hay, in Modern Theoretical Chemistry, Plenum, New York, 1976.<o:p></o:p>

5. Кларк Т., Компьютерная химия, М., «Мир» (1990).