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

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

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

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

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

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

• <SPAN style="COLOR: navy">UHF</SPAN> рассматривает спин-орбитали с различным пространственным распределением для a - и b- орбиталей. Этот метод применяется при изучении систем, как с открытыми, так и с закрытыми электронными оболочками. Так, для последних он хорошо описывает реакции диссоциации. Однако, из-за удвоения количества орбиталей, время расчета этим методом увеличивается вдвое. У этого метода существуют и другие ограничения, связанные с его основами.

• В <SPAN style="COLOR: navy">RHF</SPAN> считается, что электроны с различным спином занимают одинаковые, в смысле пространственного распределения, орбитали. При этом неспаренные электроны тоже могут занимать отдельные орбитали. Этот метод применяется как для открытых, так и для закрытых электронных оболочек.

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

Кнопка <SPAN style="COLOR: blue">OK</SPAN> служит для сохранения выбранных параметров и закрытия диалогового окна. Кнопка <SPAN style="COLOR: blue">Cancel</SPAN> приводит к закрытию диалогового окна без сохранения выбранных параметров. Как уже говорилось ранее, в программе HyperChem существует возможность определять различные базисные наборы для разных частей рассчитываемой системы.

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

<SPAN style="FONT-SIZE: 14pt; COLOR: #993300; mso-bidi-font-size: 12.0pt">Диалоговое окно параметров высокого уровня метода ab initio (Ab Initio Advanced options)<o:p></o:p></SPAN>

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

<SPAN style="COLOR: navy">Integral Format</SPAN>

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

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

Параметр <SPAN style="COLOR: navy">Cutoff</SPAN> <SPAN style="COLOR: blue"> </SPAN>позволяет сохранять на диске только те интегралы, абсолютное значение которых равно или превышает задаваемый параметр. По умолчанию он равен 10**-10 Хартри. Этот параметр контролирует осуществление SCF -итераций, точность волновых функций и энергии, так как он уменьшает количество рассчитываемых двухэлектронных интегралов.

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

<SPAN style="mso-spacerun: yes"> </SPAN> <SPAN style="COLOR: navy">Direct SCF</SPAN> <SPAN style="COLOR: navy"> calculation </SPAN> вычисляет двухэлектронные интегралы на каждой итерации, а не один раз перед SCF шагом, как если бы это было без нее. Такой расчет, безусловно, гораздо медленнее, но он позволяет не использовать дисковое пространство и операционную память под большое количество интегралов. Включение этой опции требуется при расчетах больших систем на компьютерах с маленьким диском и памятью.

<SPAN style="COLOR: navy">Ghost-atoms Control </SPAN>(<SPAN style="COLOR: navy">Использование атомов-призраков</SPAN>). Эта опция позволяет вводить центры, которым приписываются базисные функции тех или иных атомов, при этом в систему не вводятся ни дополнительные ядра, ни дополнительные электроны. Эта опция позволяет вводить только дополнительные базисные функции, центрированные в любой точке пространства. Необходимо помнить, что использование таких «виртуальных» атомов может приводить к неким артефактам. Например, так как существуют базисные волновые функции, будут существовать и мулликеновские заряды, которые будет центрированы на соответствующих центрах. Для того чтобы использовать эту возможность, необходимо в меню <SPAN style="COLOR: blue">Select</SPAN>, предварительно сделав соответствующее выделение, выбрать пункт <SPAN style="COLOR: blue">Name selection</SPAN> и, выбрав пункт <SPAN style="COLOR: blue">Other</SPAN>, приписать ему имя ghost-atoms. После этого в диалоговом меню <SPAN lang=EN-US style="COLOR: maroon; mso-ansi-language: EN-US">Ab Initio Method</SPAN><SPAN lang=EN-US style="mso-ansi-language: EN-US"> /</SPAN> <SPAN lang=EN-US style="mso-ansi-language: EN-US"> <SPAN style="COLOR: navy">Advanced Options</SPAN> </SPAN>становится<SPAN style="mso-ansi-language: EN-US"> </SPAN>активным<SPAN style="mso-ansi-language: EN-US"> </SPAN>пункт<SPAN style="mso-ansi-language: EN-US"> <SPAN lang=EN-US style="COLOR: navy">Ghost-atoms Control</SPAN> <SPAN lang=EN-US style="COLOR: navy">.</SPAN><SPAN lang=EN-US> </SPAN></SPAN>Такой подход бывает эффективным при описании ряда неординарных химических связей. Эта опция активна только для расчетов одной точки.

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

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

Кнопка <SPAN style="COLOR: blue">OK</SPAN> служит для сохранения выбранных параметров и закрытия диалогового окна. Кнопка <SPAN style="COLOR: blue">Cancel</SPAN> приводит к закрытию диалогового окна без сохранения выбранных параметров.

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

Меню <SPAN style="COLOR: blue">Compute</SPAN> определяет вид расчета. Это может быть расчет в одной точке (<SPAN style="COLOR: blue">Singl</SPAN><SPAN lang=EN-US style="COLOR: blue; mso-ansi-language: EN-US">e</SPAN><SPAN style="COLOR: blue"> Point</SPAN>), оптимизация геометрии (<SPAN style="COLOR: blue">Geometry Optimization</SPAN>), молекулярная динамика (<SPAN style="COLOR: blue">Molecular Dynamics</SPAN>), ланжевеновская динамика (<SPAN style="COLOR: blue">Langevin Dynamics</SPAN>), расчеты методом Моне-Карло, а также некоторые другие виды расчетов. Далее по тексту книги нас в основном будут интересовать результаты расчетов методом молекулярной динамики, которые представлены файлами с динамическим кино. Для того чтобы просмотреть результаты молекулярно-динамического моделирования, необходимо прочитать соответствующий файл (в тексте они будут указываться), а затем, выбрав в меню <SPAN style="COLOR: blue">Compute</SPAN> пункт <SPAN style="COLOR: blue">Molecular Dynamics</SPAN>, отметить Ц <SPAN style="COLOR: blue">Playback</SPAN> и нажать <SPAN style="COLOR: blue">Proceed</SPAN>. Для увеличения скорости проигрывания файлов (когда это необходимо, в зависимости от сложности задачи и мощности компьютера) можно увеличить <SPAN style="COLOR: blue">time steps</SPAN> (по умолчанию стоит 1, при увеличении скорость увеличивается кратно количеству шагов).

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

 

<SPAN style="FONT-SIZE: 16pt; COLOR: #993300; mso-bidi-font-size: 12.0pt">Литература</SPAN><SPAN lang=EN-US style="FONT-SIZE: 16pt; COLOR: #993300; mso-bidi-font-size: 12.0pt; mso-ansi-language: EN-US"><o:p></o:p></SPAN>

<SPAN lang=EN-US style="mso-ansi-language: EN-US">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></SPAN>

<SPAN lang=EN-US style="mso-ansi-language: EN-US">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></SPAN>

<SPAN lang=EN-US style="mso-ansi-language: EN-US">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. Act a, 28, 213 (1973); M. S. Gordon, Chem. Phys. Lett., 76, 163 (1980).<o:p></o:p></SPAN>

<SPAN lang=EN-US style="mso-ansi-language: EN-US">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></SPAN>

<SPAN lang=EN-US style="mso-ansi-language: EN-US">T.H. Dunning and P.J. Hay, in Modern Theoretical Chemistr y, Plenum, New York, 1976.<o:p></o:p></SPAN>

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