cscf - Интернет в облаке

ПРОГРАММА:

ИМЯ

cscf - решает уравнения Хартри-Фока

ОПИСАНИЕ

Программа cscf выполняет итерационную процедуру решения уравнений Хартри-Фока.

Эта программа ограничена симметрией D2h и ее подгруппами, а также орбитальными заполнениями
должны быть целыми числами. Таким образом, некоторые чистые состояния углового момента, полученные из
частичное заполнение вырожденных орбиталей не может быть получено с помощью настоящих кодов. Для
Например, состояние 2PIu (дублет PI u) линейного ONO, полученное из самой низкой энергии
линейная (pi u) 1 конфигурация может быть вычислена только как 2B2u (дублет B2u) или 2B3u
(дублет B 3u) компонента состояния 2PIu (дублет PI u), и результирующая пространственная
волновая функция не будет иметь симметрии PI. Однако в определенном смысле это желательно,
поскольку энергия будет непрерывной функцией угла изгиба. Расчет энергии
изогнутых конфигураций как 2B2u (дублет B 2u) или 2B3u (дублет B 3u) и выполнение чистого 2PIu
(дублет PI u) состояние при линейной геометрии приводит к выраженному разрыву.

По большей части триплетные состояния, возникающие в результате двукратного заполнения дважды вырожденного
орбитальное, такое как состояние 3A2 (триплет A 2), возникающее в результате (e ') 2 или (e ") 2
конфигурации в симметрии D3h, или 3SIGMAg (триплет SIGMA g) состояние (pi g) 2 или (pi
u) 2 в симметрии Dinfh (D бесконечность h), будет иметь правильную пространственную симметрию.
Синглетные состояния, возникающие в результате тех же электронных конфигураций, по своей природе
многоконфигурация и, как таковая, плохо представлены одной конфигурацией
волновые функции.

Ссылки

PK-файловый метод:

1. RC Raffenetti, Chem. Phys. Lett. 20 (1973) 335.

Молекулярная симметрия и расчеты HF с закрытой оболочкой:

1. M.Dupuis, HFKing, Int. J. Quant. Chem. 11 (1977) 613.

DIIS для закрытой оболочки:

1. P. Pulay, Chem. физ. лат. 73 (1980) 393.

2. P. Pulay, J. Comp. хим. 3 (1982) 556.

Коэффициенты связи (альфа и бета) для открытой оболочки:

1. CCJ Roothaan, Rev. Mod. Phys. 32 (1960) 179.

Демпфирование:

1. Д. Р. Хартри, "Расчет атомных структур" (Wiley: Нью-Йорк), 1957.

2. MC Zerner, M. Hehenberger, Chem. Phys. Lett. 62 (1979) 550.

Смещение уровня:

1. VR Saunders и IH Hillier, Int. J. Quant. Chem. 7 (1973) 699.

ОБЪЕДИНЕНИЕ CSCF

Для сложных открытых корпусов рекомендуется использовать соответствующий закрытый корпус.
вычисление должно быть выполнено первым (добавить или удалить лишний электрон) и что этот вектор SCF затем
использоваться в качестве предположения для желаемой волновой функции открытой оболочки. Для случаев TCSCF всегда
разумно сначала запустить закрытую оболочку (или, возможно, соответствующий триплет) SCF, а затем использовать
это предположение для TCSCF.

Для систем с открытой оболочкой рекомендуется значение сдвига уровня от 0.5 до 3.0. Начните с
высокое значение (2.0 - 3.0) для первого расчета SCF, а затем уменьшите его (до 0.5 - 1.0) для
последующие прогоны, в которых в качестве отправной точки используется сходящийся вектор SCF.

Чрезвычайно важно отметить, что эта версия кода больше не поддерживает
ОТКРЫТЫЙ. Одна обязательно используют новый ключевые слова СПРАВКА и МУЛЬТП в указывать напишите of SCF
необходимо.

ВХОД ФОРМАТ

The cscf программа выполняет поиск по пути ключевого слова по умолчанию (сначала SCF , а затем DEFAULT)
для следующих ключевых слов:

LABEL = string
Это строка символов, которая будет включена в вывод. Эта строка не используется
по программе. По умолчанию нет.

ВФН = string
Это тот тип волновой функции, который в конечном итоге желателен. По умолчанию SCF.

ОТКРЫТЫЙ is нет дольше поддержал

СПРАВКА = string
Это определяет тип вычисления SCF, который нужно выполнить. Это может быть один из РГНФ
(для синглета с закрытой оболочкой), РОХФ (для расчета ограниченной открытой оболочки), UHF
(для неограниченного расчета открытой оболочки), ДВАКОН (для двух конфигураций
синглет), или SPECIAL. Если SPECIAL дается, то альфа- и бета-связь
коэффициенты должны быть заданы с АЛЬФА и BETA ключевые слова. По умолчанию РГНФ.

MULTP = целое
Задает кратность молекулы. По умолчанию синглет.

ЗАРЯД = целое
Задает заряд молекулы. По умолчанию 0.

ДОКК = целочисленный_вектор
Это дает количество дважды занятых орбиталей в каждой неприводимой
представление. По умолчанию нет. Если это не указано, CSCF попытается
угадайте профессии, используя основной гамильтониан.

SOCC = целочисленный_вектор
Это дает количество однократно занятых орбиталей в каждой неприводимой
представление. По умолчанию нет.

ДЕРТИП = string
Это определяет порядок производной, которая должна в конечном итоге быть сделана. Это используется
не провела обыск SCF программа, чтобы определить, должны ли быть записаны определенные файлы, а также
используется для определения сходимости волновой функции по умолчанию. По умолчанию
ПЕРВЫЙ.

МАКСИТЕР = целое
Это дает максимальное количество итераций. По умолчанию - 40.

КОНВЕРГЕНЦИЯ = целое
Это определяет, насколько сильно волновая функция будет сходиться. Конвергенция
определяется путем сравнения среднеквадратичного изменения матрицы плотности ("дельта P") с
данное значение. Критерий сходимости - 10 ** (-целое). По умолчанию 7, если
и ДЕРТИП = NONE и ВФН = SCF даны и 10 в противном случае.

СДВИГ УРОВНЯ = реальные
Это определяет сдвиг уровня. По умолчанию 1.

НЕПОСРЕДСТВЕННЫЙ = логический
Указывает, следует ли выполнять вычисление SCF с помощью прямого интегрального метода. В
по умолчанию - ложь.

PRINT_MOS = логический
Указывает, печатать ли молекулярные орбитали или нет. По умолчанию - false.

Также существует большое количество менее часто используемых входных параметров. Если вы этого не сделаете
понять, что означают следующие параметры, затем убедитесь, что они не отображаются в вашем
Вход. Значения по умолчанию будут работать в подавляющем большинстве случаев. Они указаны
со следующими ключевыми словами:

DELETE_INTS = логический
Файлы интегралов будут удалены, если ВФН = SCF и ДЕРТИП = ПЕРВЫЙ or ДЕРТИП = NONE.
Если вы хотите сохранить файлы интегралов, установите DELETE_INTS = ложь. По умолчанию
правда.

ИЗМЕНЕНИЕ ПОРЯДКА = string
Параметр управляет переупорядочением молекулярных орбиталей. Если установлено на ДО то
угаданные орбитали из файла контрольной точки будут переупорядочены. Если установлено на ПОСЛЕ, сошлись
орбитали будут переупорядочены перед записью в файл контрольной точки. В любом
, признали МОРДЕР необходимо указать параметр, чтобы указать карту переупорядочения. По умолчанию
не менять порядок орбиталей.

МОРДЕР = целочисленный_вектор
Это определяет вектор переупорядочения молекулярных орбиталей. Он будет использоваться только в том случае, если
ИЗМЕНЕНИЕ ПОРЯДКА установлен. Этот вектор отображает каждую орбиталь на ее новый индекс, например МОРДЕР = (0
2 1) указывает, что после переупорядочения орбитали 1 и 2 поменяются местами. Ранг
этот вектор совпадает с количеством МО. Индексы в порядке Питцера.
(упорядочено по симметрии, затем по энергии в каждом блоке симметрии), base-0. CSCF
скорее всего потерпит неудачу, если данный MOORDER смешивает орбитали из разных репсов. Там
не по умолчанию.

АЛЬФА = реальный_вектор
If ОТКРЫТЫЙ = SPECIAL, то этот параметр дает коэффициенты альфа-связи.
Количество элементов в этом векторе - MM (MM + 1) / 2, где MM - количество
неприводимые представления, содержащие однократно заполненные молекулярные орбитали. Там
не по умолчанию.

BETA = реальный_вектор
If ОТКРЫТЫЙ = SPECIAL, то этот параметр дает коэффициенты бета-связи.
Количество элементов в этом векторе - MM (MM + 1) / 2, где MM - количество
неприводимые представления, содержащие однократно заполненные молекулярные орбитали. Там
не по умолчанию.

УГАДАЙ = string
Эта опция определяет тип начального предположения о собственном векторе, который будет использовать CSCF.
Единственные допустимые варианты на данный момент: (1) УГАДАЙ = CORE, что заставляет его использовать
базовый собственный вектор гамильтониана для начала расчета; (2) УГАДАЙ = AUTO , которые
приводит к попытке использовать вектор MO в файле контрольной точки или прибегает к
основное предположение, если в этом файле нет собственного вектора. По умолчанию, если AUTO.

ПЕЧАТЬ = целое
Это вариант печати. По умолчанию - 0.

MO_OUT = логический
Распечатывает орбитали с симметрией и заполнениями в конце
расчет. По умолчанию верно.

ПОВОРОТ = логический
Молекулярные орбитали не будут повернуты, если это неверно. Только вращение
влияет на виртуальные орбитали для систем с открытыми оболочками. Этот параметр должен быть истинным
для коррелированных градиентов, и он должен быть ложным для вторых и более высоких производных.
По умолчанию false, если ВФН = SCF и верно в противном случае.

CHECK_ROT = логический
Проверьте вращение молекулярной орбиты, описанное выше, чтобы убедиться, что нет столбцов с
матрица собственных векторов SCF меняются местами при вращении. Не действует, если ПОВОРОТ =
ложный. По умолчанию это правда.

CHECK_MO_ORTHOGONALITY = логический
Проверьте, ортонормированы ли молекулярные орбитали. Полезно только для отладки. В
по умолчанию - ложь.

ДИИС = логический
Это определяет, будет ли использоваться diis. По умолчанию это правда.

DIISTSTART = целое
Это дает первую итерацию, для которой будет использоваться DIIS. По умолчанию - 0.

НДИИС = целое
Это дает количество матриц ошибок для использования в процедуре diis. По умолчанию
равно 6 для закрытой оболочки, 4 для открытой оболочки и 3 для tcscf.

ДИИСДАМП = реальные
Это дает коэффициент демпфирования для процедуры diis. По умолчанию для
закрытая оболочка, 0.02 для открытой оболочки и 0.01 для tcscf.

ИНКР = реальные
Это используется в tcscf, чтобы определить, как часто пересчитываются коэффициенты ci.
Небольшое число (~ 0.25) заставит их пересчитывать почти каждый scf
итерация. По умолчанию 0.25.

ДИН_АКК = логический
При выполнении прямого scf это указывает, будут ли отсечки динамической интегральной точности
будет использоваться. По умолчанию true (используйте динамические отсечки). Начальные итерации
выполняется с интегралами с точностью до шести знаков. После того, как плотность сведена к
Выполнено 10 ^ -5 или 30 итераций, используется полная интегральная точность. Если scf
проблемы сходимости возникают при отключении динамической отсечки путем установки этого
значение переменной false может помочь.

ORTHOG_ONLY = логический
Иногда в CASSCF или других схемах, отличных от HF / KS для орбитальной оптимизации, это
полезно для реортогонализации МО из других геометрий для текущей геометрии, чтобы
их можно использовать в качестве первоначального предположения для новых МО. Это может быть выполнено
запуск CSCF с ORTHOG_ONLY = правда. После ортогонализации орбиталей
программа завершится без выполнения вычисления SCF. Это ключевое слово будет
игнорируется, если в файле контрольной точки нет предыдущих орбиталей. По умолчанию правда
if ВФН = ДЕТКАС.

30 Май, 1991 cscf(1)

Используйте cscf онлайн с помощью сервисов onworks.net