Не получается запустить CCSD(T) расчёт

[Vit Nhoc]

Я запускаю расчёт оксида никеля с водородом (изучаем каталитические свойства оксида никеля). Расчёты на DFT (PBE), MP2 дают разную дурь. Попробовал запустить CCSD(T)/cc-pVTZ, после 50 вторичных итераций выдаётся "*MAX. CYCLES*".
Входной файл:

%MEM=2000MW
#P CCSD(T)/cc-pVTZ SCF(XQC) OPT(Z-Matrix)

NiO

0,3
O
Ni 1 R12
Ni 1 R13 2 A213
O 2 R24 1 A124 3 D3124
H 2 R25 4 A425 3 -D3425
H 5 R56 2 A256 1 D1256
Variables:
R12 = 1.78256368
R13 = 1.69813043
A213 = 75.50025433
R24 = 1.78256166
A124 = 100.91079339
D3124 = 0.00155249
R25 = 1.59057116
A425 = 127.79696805
D3425 = 160.01706249
R56 = 0.85884182
A256 = 74.34054345
D1256 = 102.54483074


Можете что-нибудь посоветовать? В частности, как в Гауссиане сделать выдачу тау-критерия?

[Гесс]

попробуйте CCSD(T)=Maxcyc=500/cc-pVTZ

что вы имеете в виду под тау-критерием?

[Vit Nhoc]

попробуйте CCSD(T)=Maxcyc=500/cc-pVTZ

что вы имеете в виду под тау-критерием?

Это мне знакомый квантовик сказал (он считает на Molpro), тау-критерий служит индикатором адекватности CC-расчёта.
Делать больше циклов не имеет смысла - я посмотрел сходимость CC, там энергия болтается вверх-вниз.

[YuraM]

Насчет количества итераций - попробуйте ХФ волновую функцию посчитать с SCF=tight. Также проверьте стабильность волновой функции stable=opt.

Проверить адекватность - два критерия:
1) T1 diagnostic:

T1Diag
Computes the T1 diagnostic of T. J. Lee and coworkers [Lee89] (CCSD only).

2) Печатайте наибольшую амплитуду в CCSD.

А у нас в ИГХТУ кто-то занимается каталитическими свойствами оксида никеля?

[Vit Nhoc]

А у нас в ИГХТУ кто-то занимается каталитическими свойствами оксида никеля?

[cut]Как я понял, на кафедре физхимии занимаются скорее скелетным никелем, но он всегда покрывается оксидной плёнкой[/cut]

[Гесс]

Я прогнал ваш инпут с CCSD(T)=Maxcyc=500/cc-pVDZ (даблзета ибо мне не хочется ждать триплзета). Оно делает по 81 итерации и переходит на следующий шаг. Когда сделает первый шаг оптимизации - убью и выложу сюда.
Там прослеживается довольно стабильное снижение Delta, и когда оно достигает 1*10^-7 процедура считается сошедшейся.

Я не уверен насколько поведение компактной фазы может быть описано столь малыми кластерами.

[Vit Nhoc]

Гесс,
А мне только что сказали что такую задачу нереально решать на CCSD(T), поскольку она мультиконфигурационная.
Сказали, нужен CASSCF или аналоги, но на персональном компьютере просчитать это нереально - слишком много неспаренных электронов, которые надо включить в активное пространство.

[Гесс]

Гесс,
А мне только что сказали что такую задачу нереально решать на CCSD(T), поскольку она мультиконфигурационная.
Сказали, нужен CASSCF или аналоги, но на персональном компьютере просчитать это нереально - слишком много неспаренных электронов, которые надо включить в активное пространство.

это вполне может быть. что там с активным пространством я сказать не возьмусь (и уж тем более не глядя) но кас-расчет болеее требователен к памяти чем к числу ядер. Светит вам путь в Firefly или на самый худой конец в Орку. Гауссиан для этих целей не подойдет.

[Shorku]

Структурка довольно странная сама по себе. Но если это зачем-то надо считать, попробуйте нарисовать диаграмму МО и оценить, единственным ли образом можно раскидать электроны, или нет.

[Vit Nhoc]

Я прогнал ваш инпут с CCSD(T)=Maxcyc=500/cc-pVDZ (даблзета ибо мне не хочется ждать триплзета). Оно делает по 81 итерации и переходит на следующий шаг. Когда сделает первый шаг оптимизации - убью и выложу сюда.
Там прослеживается довольно стабильное снижение Delta, и когда оно достигает 1*10^-7 процедура считается сошедшейся.

Я не уверен насколько поведение компактной фазы может быть описано столь малыми кластерами.

Ну как, посчиталось?

[Гесс]

Я приношу свои извинения. Я решил что вы отказались от использования Каплдкласера для этой задачи и что результат не представляет для вас интереса. Расчет умер существенно позже и я не искал как такое можно лечить.
[cut]WARNING| EIGENVALUE 12 TOO LARGE. REPLACED BY 25.000000
WARNING| EIGENVALUE 11 TOO LARGE. REPLACED BY 25.000000
WARNING| EIGENVALUE 10 TOO LARGE. REPLACED BY 25.000000
WARNING| EIGENVALUE 9 TOO LARGE. REPLACED BY 25.000000
WARNING| EIGENVALUE 8 TOO LARGE. REPLACED BY 25.000000
HESSIAN DOES NOT HAVE THE DESIRED LOCAL STRUCTURE
TAKING SIMPLE RFO STEP
SEARCHING FOR LAMDA THAT MINIMIZES ALONG ALL MODES


*****************************************
*** UNABLE TO DETERMINE LAMDA IN FmD114 **
*****************************************


Error termination via Lnk1e in /usr/local/g09/l114.exe at Tue Jul 5 19:35:41 2016.[/cut]

VitCC500.out.txt