detcas - онлайн у хмарі

ПРОГРАМА:

ІМ'Я

detcas - Взаємодія детермінанта конфігурації / Повний активний
Програма космічного самоузгодженого поля

ОПИС

Програма detcas виконує самоузгоджене поле повного активного простору (CASSCF)
оптимізація молекулярних орбіталей за допомогою двоетапної процедури в поєднанні з
програма взаємодії визначника конфігурації DETCI. Програма досить проста і
в даний час використовує підхід Ньютона-Рафсона для оновлення орбіталей, використовуючи простий
наближений орбітальний гессенський. Зближення прискорюється за допомогою прямої інверсії Пулея
процедура ітераційного підпростору (DIIS). Код був написаний, щоб дозволити більш загальні
хвильові функції, які не обов’язково мають повну обробку CI активного простору. в
зокрема, будь-яка хвильова функція CI з обмеженим активним простором (RAS), підтримувана DETCI, може бути
використовується, що дозволяє хвильові функції MCSCF типу RASSCF.

Посилання

Приблизний орбітальний гессиан:

1. Г. Чабан, М. В. Шмідт, М. С. Гордон, Теор. Chim. Acta 97, 88-95 (1997).

Обмежений активний простір CI:

1. Алгоритми взаємодії конфігурації на основі детермінанта для повних і обмежених
Configuration Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen, and HJ Aa.
Jensen, J. Chem. фіз. 89, 2185 (1988).

Обмежений активний простір SCF:

1. П.-А. Malmqvist, A. Rendell, and BO Roos, J. Phys. хім. 94, 5477 (1990).

Програма DETCI:

1. CD Sherrill, Обчислювальні алгоритми для великомасштабного повного та багатопосилкового посилання
Хвильові функції конфігураційної взаємодії, докторська дисертація, Університет Джорджії, Афіни,
ГА, 1996.

2. CD Sherrill і HF Schaefer, Метод взаємодії конфігурації: прогрес
у висококорельованих підходах, Adv. Quantum Chem. 34, 143-269 (1999).

ФАЙЛИ ВИМАГАЄТЬСЯ

input.dat - вхідний файл
file78 - МО одноелектронні інтеграли (оператор fzc)
файл72 - МО двоелектронні інтеграли
file73 - одночастинкова матриця щільності МО
file74 - МО двочасткова матриця щільності
файл75 - М. О. Лагранж

ФАЙЛИ ОНОВЛЕНО

output.dat - вихідний файл
file14.dat - Запис енергій і орбітальних градієнтів

ВХІД ФОРМАТ

Доступні такі аргументи командного рядка:

-спокійно Це дає той самий результат, що і PRINT=0.

-o fname
Дає ім’я файлу для вихідного файлу. За замовчуванням – output.dat.

Додатковий вхід для цієї програми зчитується з файлу input.dat. Чим частіше використовується
ключові слова:

WFN = рядок
Підтримувані типи хвильових функцій: CASSCF і RASSCF.

КОНВЕРГЕНЦІЯ = ціле
Бажана конвергенція на орбітальному градієнті. Збіжність досягається при RMS
похибки орбітального градієнта менше 10**(-n). За замовчуванням 4 for
розрахунки енергії та 7 для градієнтів.

DOCC = цілочисельний_масив
Цей вектор дає кількість подвійно зайнятих орбіталей у кожній ірреп. існує
за замовчуванням немає.

SOCC = цілочисельний_масив
Цей вектор дає кількість однократно зайнятих орбіталей у кожній ірреп. існує
за замовчуванням немає.

ЕНЕРГІЯ_КОНВЕРГЕНЦІЯ = ціле
Бажана збіжність по повній енергії MCSCF. За замовчуванням 7.

FROZEN_DOCC = цілочисельний_масив
Кількість найнижчих енергетичних подвійно зайнятих орбіталей у кожній з незводних
подання, яке буде буквально заморожено (не оновлюється в MCSCF). The
Використовується бавовняне впорядкування незмінних уявлень. За замовчуванням є
нульовий вектор.

FROZEN_UOCC = цілочисельний_вектор
Кількість незайнятих орбіталей з найвищою енергією в кожному незводному представленні
який буде буквально заморожений (не оновлюється в MCSCF). За замовчуванням є нуль
вектор.

RESTRICTED_DOCC = цілочисельний_масив
Кількість найнижчих енергетичних подвійно зайнятих орбіталей у кожній з незводних
представлення, яке буде оптимізовано, але подвійно зайняте в MCSCF.
Ці орбіталі приходять після FROZEN_DOCC орбіталі. За замовчуванням є нуль
вектор.

RESTRICTED_UOCC = цілочисельний_масив
Кількість незайнятих орбіталей з найвищою енергією в кожному незводному представленні
які будуть оптимізовані, але не зайняті в MCSCF. Ці орбіталі приходять
перед FROZEN_UOCC орбіталі. За замовчуванням є нульовий вектор.

NCASITER = ціле
Максимальна кількість ітерацій для оптимізації орбіталей. Цей варіант Повинен be
зазначений in ПОВЕРНЕНО розділ of вхід, тому що це має бути видимим для
керуюча програма PSI. За замовчуванням 1.

ДРУК = ціле
Цей параметр визначає багатослівність виводу. Значення 1 або 2 визначає
мінімальний друк, значення 3 визначає детальний друк. Значення 4 або 5 є
використовується для налагодження. Не використовуйте рівень 5, якщо тестовий випадок не дуже маленький (наприклад
STO H2O CISD).

Рідше використовувані ключові слова:

DIIS_FREQ = ціле
Процедура екстраполяції DIIS буде спробована щоразу n ітерації. За замовчуванням
є 1.

DIIS_MAX_VECS = ціле
Максимальна кількість векторів підпростору для процедури DIIS. Після цього числа
векторів, старі вектори будуть вилучені з підпростору як
необхідно. За замовчуванням – 8.

DIIS_MIN_VECS = ціле
Мінімальна кількість векторів підпростору перед інтерполяцією DIIS може бути
виконано. За замовчуванням встановлено значення 2. Немає сенсу мати значення менше 2.

DIIS_START = ціле
Номер ітерації під час прямої інверсії ітеративного підпростору Пулея (DIIS)
процедуру прискорення конвергенції слід увімкнути. За замовчуванням 3.
До цієї ітерації вектори не додаються до підпростору DIIS. Перший
Крок DIIS не буде зроблено доти DIIS_MIN_VECS вектори знаходяться в підпросторі DIIS.

LAG_FILE = ціле
Файл (номер одиниці) для читання матриці Лагранжа. Значення за замовчуванням
на даний момент 75.

OPDM_FILE = ціле
Файл (номер одиниці) для зчитування одночасткової матриці щільності. Значення за замовчуванням
наразі 73.

SCALE_GRAD = boolean
Вказує, чи потрібно масштабувати орбітальний градієнт за наближеною (діагональною) орбіталлю
Гессенський. За замовчуванням є TRUE.

SCALE_STEP = реальний
Масштабний коефіцієнт для кроку орбітального обертання. За замовчуванням 1.0.

TPDM_FILE = ціле
Файл (номер одиниці) для зчитування двочасткової матриці щільності. Значення за замовчуванням
наразі 74.

LEVEL_SHIFT = boolean
Вказує, чи допускати зсув рівня матриці Гессе. Цим можна скористатися
щоб переконатися, що гессиан є позитивно визначеним для початкових ітерацій. Якщо
гессиан не є позитивно визначеним поблизу збіжності, то це може бути ознакою
чисельні нестабільності в MCSCF. За замовчуванням є TRUE.

SHIFT = реальний
Вказує, наскільки діагональні елементи гессена повинні бути зміщені для рівня
переміщення. За замовчуванням 0.01.

DETERM_MIN = реальний
Мінімально дозволене значення для визначника матриці Гессе при зсуві рівня
увімкнено. Це не використовується, якщо LEVEL_SHIFT=FALSE. За замовчуванням 0.00001.

STEP_MAX = реальний
Це максимально дозволений один орбітальний оберт. За замовчуванням 0.30.

FORCE_STEP = boolean
Це дозволяє користувачеві змінити обчислений крок і примусово виконати крок у a
конкретний напрям. Це може бути корисно, якщо ви намагаєтеся вигнати обчислення
від сідлової точки. За замовчуванням встановлено значення FALSE.

FORCE_PAIR = ціле
Це індекс для незалежної пари, яку потрібно повернути, якщо
FORCE_STEP=TRUE. Інакше це ігнорується. За замовчуванням немає.

FORCE_VALUE = реальний
Це значення орбітального обертання для обертання орбіталей, визначених як
FORCE_PAIR. Це використовується, лише якщо для FORCE_STEP встановлено значення TRUE. Значення за замовчуванням
0.0.

CHECK_HESSIAN = boolean
Обчисліть власні значення орбітального гессена. Це добре для перевірки
сідловидні рішення. Слід використовувати без зміни рівня. За замовчуванням є
ПОМИЛКОВИЙ.

EIGEN_VECTORS = boolean
Ви хочете також обчислити власні вектори орбітального гессеана. Це
ігнорується, якщо CHECK_HESSIAN=FALSE. За замовчуванням – FALSE.

8 травень, 1998 detcas(1)

Використовуйте detcas онлайн за допомогою служб onworks.net