Выбор типа расчёта

В GAMESS (US) тип расчёта (то есть, его конечную цель расчёта и алгоритм пути к ней) определяет параметр RUNTYP. Всего в GAMESS (US) доступно 34 значения. О каждом из них них можно почитать в официальном мануале. Здесь же будут рассмотрены только некоторые основные значения.

Важно

Все типы расчётов можно разделить на две группы: градиентные и неградиентные.

К градиентным относятся расчёты, задаваемые значениями RUNTYP=GRADIENT, HESSIAN, OPTIMIZE, SADPOINT, GLOBOP, IRC, GRADEXTR, DRC и RAMAN. Остальные относятся к неградиентным.

Градиентные расчёты наиболее эффективны, когда присутствует возможность вычислять градиент аналитически, что доступно далеко не для всех методов. Численное вычисление градиентов в GAMESS (US) доступно для всех методов и активируется значением NUMGRD=.TRUE. в группе $CONTRL.

=ENERGY – рассчёт полной энергии системы, геометрия частицы при этом не изменяется (установлено по умолчанию).

Примечание

Расчёты, которые не предполагают изменение геометрии системы, также известны как single point расчёты.

=GRADIENT – то же, что и RUNTYP=ENERGY, но дополнительно рассчитывается градиент.

=HESSIAN – то же, что и RUNTYP=GRADIENT, но дополнительно рассчитывается гессиан и проводится анализ колебаний в гармоническом приближении. Для аналитического расчёта гессиана в FMO расчётах применяется другое значение – RUNTYP=FMOHESS.

=OPTIMIZE – градиентная оптимизация геометрии системы в минимум энергии.

=TRUDGE – неградиентная оптимизация геометрии системы в минимум энергии. Данное значение также может быть использовано для неградиентной оптимизации базисного набора.

=GLOBOP – оптимизация геометрии системы в минимум энергии методом Монте-Карло.

=SADPOINT – градиентная оптимизация геометрии системы в седловую точку.

Значения RUNTYP=IRC и RUNTYP=DRC отвечают за поиск пути наименьшей энергии реакции, разница между ними заключается в том, что при значении RUNTYP=IRC кинетическая энергия на всех этапах приравнивается к нулю, а при RUNTYP=DRC сумма кинетической и потенциальной энергий остаётся постоянной на протяжении всего расчёта.