orca ccsd diis

[Гесс]

Доброе время суток!
Делается оптимизация масюсенького кластера в каплкластере.

Очень быстро наступает:

------------------
CLOSED SHELL GUESS
------------------

Initial guess performed in 0.013 sec
E(0) ... -301.329105739
E(MP2) ... -1.148269331
Initial E(tot) ... -302.477375070
<T|T> ... 0.801826836
Number of pairs included ... 136
Total number of pairs ... 136

------------------------------------------------
RHF COUPLED CLUSTER ITERATIONS
------------------------------------------------

Number of amplitudes to be optimized ... 1054592

Iter E(tot) E(Corr) Delta-E Residual Time
0 -302.477375070 -1.148269331 0.000000000 0.019094053 1.62
*** Turning on DIIS ***
1 -302.349626800 -1.020521061 0.127748269 0.014299186 1.78
2 -302.435376948 -1.106271209 -0.085750147 0.002997234 2.20
3 -302.463759150 -1.134653411 -0.028382202 0.004127612 2.40
4 -302.475770174 -1.146664435 -0.012011024 0.002898209 2.60
5 -302.484828961 -1.155723222 -0.009058787 0.002777329 3.32
6 -302.495874760 -1.166769021 -0.011045799 0.002082174 3.04
7 -302.502926165 -1.173820426 -0.007051405 0.001761169 3.20
8 -302.507674554 -1.178568815 -0.004748389 0.001262132 3.16
9 -302.512595544 -1.183489805 -0.004920990 0.001035519 3.20
10 -302.514768322 -1.185662583 -0.002172777 0.000792015 3.15
11 -302.517541548 -1.188435809 -0.002773226 0.000732229 3.14
12 -302.518244686 -1.189138947 -0.000703138 0.000652266 3.19
13 -302.519405844 -1.190300105 -0.001161158 0.000603942 3.17
14 -302.520045250 -1.190939511 -0.000639406 0.000606257 3.16
15 -302.520764205 -1.191658466 -0.000718955 0.000663196 3.19
16 -302.521204583 -1.192098844 -0.000440377 0.000701938 3.20
17 -302.521710370 -1.192604631 -0.000505787 0.000741500 3.19
18 -302.522052761 -1.192947022 -0.000342391 0.000744690 3.19
19 -302.522360299 -1.193254560 -0.000307538 0.000747105 3.19
20 -302.522530428 -1.193424689 -0.000170129 0.000731934 3.73
21 -302.522637033 -1.193531294 -0.000106604 0.000735727 3.19
22 -302.522646161 -1.193540422 -0.000009128 0.000742025 3.42
23 -302.522634886 -1.193529147 0.000011275 0.000749547 3.49
24 -302.522665800 -1.193560061 -0.000030914 0.000756972 3.20
25 -302.522616861 -1.193511122 0.000048939 0.000754013 3.21
26 -302.522670365 -1.193564627 -0.000053505 0.000758382 3.18
27 -302.522734542 -1.193628803 -0.000064176 0.000762716 3.17
28 -302.522782688 -1.193676949 -0.000048146 0.000766500 3.48
29 -302.522782426 -1.193676687 0.000000261 0.000767619 3.20
30 -302.522808648 -1.193702909 -0.000026222 0.000770229 3.18
31 -302.522826193 -1.193720454 -0.000017544 0.000772116 3.19
32 -302.522846809 -1.193741070 -0.000020616 0.000773935 3.19
33 -302.522896521 -1.193790782 -0.000049712 0.000778171 3.18
34 -302.522930251 -1.193824512 -0.000033729 0.000781019 3.16
--- The Coupled-Cluster iterations have NOT converged ---
************************************************************
* Shut down parallel processing *
************************************************************
-------
TIMINGS
-------

Расчет на этом не умирает (он умирает позже и причины его умирания мне пока непонятны), но собственно хотелось бы решить эту несходимость.
Она не реагирует на

%scf
DIISMaxit=100
Maxiter=200
end

(Собственно и не должно реагировать, блок "SCF ITERATIONS" с горем пополам сходится перед этим)
Равно как и на NODIIS SOSCF в ключах

Ошибка одинакова как при использовании перегруженного инпута

! CCSD(T) cc-pvdz cc-pvdz/J RIJCOSX TightSCF VeryTightSCF
! NODIIS SOSCF NumFreq NUMGRAD Opt TightOpt Grid6
%pal nprocs 32 end
%geom Maxiter=100 end
%scf
#DIISMaxit=100
Maxiter=200
end

так и при бедненькой версии

! CCSD cc-pvdz VeryTightSCF TightOpt
! KDIIS Grid6
%pal nprocs 32 end
%geom Maxiter=100 end
%scf
#DIISMaxit=100
Maxiter=200
end
* xyz 0 1
Координаты
*

Собственно вопрос - как можно увеличить число итераций в данном блоке? либо как в него впихнуть альтернативные DIISу методы?, ну либо что еще можно сделать?

[kbob]

Настройки Coupled Cluster расчета задаются в %mdci блоке:
Вот те, что могут относится к вашему случаю с дефотными значениями:

Код: Выделить всё

MaxIter 35  # Max. number of iterations
MaxDIIS 7   # Max number of DIIS vectors to be stored

отказаться от DIIS тут видимо нельзя.
Но советую отказаться от RIJCOSX, из-за него обычно сильно ухудшается сходимость.
Если вам нужно уменьшить время расчета, то лучше использовать LPNO приближениие:

Код: Выделить всё

! RHF LPNO-CCSD cc-pVTZ cc-pVTZ/C TightSCF

[Гесс]

Спасибо! Убежал пробовать!

[alien308]

Здесь была глупость.

[Гесс]

Здесь была глупость.

любая глупость может быть полезна.

Еще раз большое спасибо kbobу, конкретно данная проблема успешно решена. Расчет валится в других местах с какими то совсем невнятными причинами (или пишет об ошибках но все равно шурует вперед и вперед, короче загадочно), но об этом если не решу сам - создам тему через пару дней.
LPNO мне не подойдет, ибо в норме моя задача должна считаться в CCSD(T), а вот с (T) орчье ЛПНО пока не дружит.
Был пойман странный идиотизм - вроде ничего кроме mdci блока не менял, но внезапно полезли проблемы со сходимостью в SCF блоке, причем в одном случае проблема рещилась с горем пополам увеличением числа итераций, а другая система янтеллеровская и чето теперь я ейный SCF не могу свести даже с применением NRSCF (там оно опять же странно рушится, подробнее распишу позже если не поборю).

[endlesslake]

в третьей орке появилось DLPNO-CCSD(T), правда только closed shell.

[Гесс]

в третьей орке появилось DLPNO-CCSD(T), правда только closed shell.

задача большей частью в open-shelle.
но всеравно спасибо, даст бог где-нибудь пригодится

[kbob]

На странице 66 доки для 3.0.0 написано:

There are a few things to be noticed about (D)LPNO methods:
• The LPNO methods obligatorily make use of the RI approximation. Hence, a correlation fit set must be
provided.
• For DLPNO-CCSD presently only the closed-shell variants are available.
• Besides the closed-shell version which uses a RHF or RKS reference determinant there is an open-shell version
for high-spin open-shell molecules (see original paper) using an UHF or UKS reference determinant build
from quasi-restricted orbitals (QROs, see section 7.8.3). Since the results of the current open-shell version are
slightly less accurate than that of the closed-shell version it is mandatory to specify if you want to use the
closed-shell or open-shell version for calculations of closed-shell systems, i.e. always put the “RHF” (“RKS”)
or “UHF” (“UKS”) keyword in the simple keyword line. Open-shell systems can be of course only treated by
the open-shell version. Do not mix results of the closed- and open-shell version (e.g. if you calculate
reaction energies of a reaction in which both closed- and open-shell molecules take part, you should use the
open-shell version throughout).
• Triples will become available in the next ORCA release.
• F12-variants will become available in the foreseeable future.
• Parallelization is done.

ну а если у вас что-то валится то прейдите на версию 3.0.0, а если всеравно валится задавайте вопросы либо здесь либо на форуме ORCA.

[Гесс]

ну а если у вас что-то валится то прейдите на версию 3.0.0, а если всеравно валится задавайте вопросы либо здесь либо на форуме ORCA.

В свете моих "успехов" в установке и распараллеливании КХ-программ под линуксом это звучит как проклятье.
Но это уже поставлено в мои задачи.
Возможно проблема решится с другой стороны - у меня уже исчезла задача оптимизации Open-shell систем в CCSD(T), так что возможно жить станет легче. Посмотрим.

[alien308]

Здесь была глупость.

любая глупость может быть полезна.

Извиняюсь, ввёл в заблуждение. Глупость была моя. Ответил не вполне правильно не прочитав ответа уважаемого kbob.

[Гесс]

Да я понял что у вас была самокритика, просто возможно в написанной вами глупости было чтото что мне либо кому нибудь другому бы пригодилось ))).
Ладно, проехали ))

[VTur]

Возможно проблема решится с другой стороны - у меня уже исчезла задача оптимизации Open-shell систем в CCSD(T), так что возможно жить станет легче. Посмотрим.

Гесс, единствено, что хочется заметить, такие методы удовлетворительно работают в больших и сверхбольших базисах. Не много молекул можно ими посчитать. Конешно, можно посчитать и опубликовать расчет в существенно меньшем базис, но стоит ли?

[ Post made via Android ]

[Гесс]

Гесс, единствено, что хочется заметить, такие методы удовлетворительно работают в больших и сверхбольших базисах. Не много молекул можно ими посчитать. Конешно, можно посчитать и опубликовать расчет в существенно меньшем базис, но стоит ли?

Для клоздшелла уже есть оптимизации в cc-pvdz и aug-cc-pvtz, которые мало отличаются. Это было сделано в молпре, и там есть странные фигни которые мне ненравятся. И вот пришел я в орку и тут все тоже плохо. Короче куда ни кинь всюду нет маланцу счастья.

[VTur]

Если в расчете ccsd есть странности, уберите букву s и посмотрите, может они исчезнут?

[ Post made via Android ]

[Гесс]

Если в расчете ccsd есть странности, уберите букву s и посмотрите, может они исчезнут?

тут много разных букв можно убрать и добавить:
CCD, OOCCD, pCCSD, токо для каждого надо разбираться стоит ли мне в него лезть и не получу ли я там дурь на выходе.

[VTur]

Возьмите CCD, уберете однократные возбуждения. Может дело в них

[ Post made via Android ]

[Гесс]

Возвращаясь к нашим баранам:
Делается клоздшеллнутое DLPNO/CCSD(T)

! DLPNO-CCSD(T) cc-pvdz cc-pVDZ/C Opt Freq
! KDIIS VeryTightSCF TightOpt Grid6 Numgrad
%pal nprocs 16 end
%geom Maxiter=100 end
%scf
DIISMaxit=100
Maxiter=500
end
%mdci
Maxiter=500
end
* xyz 0 1

RHF COUPLED CLUSTER ITERATIONS

сводятся довольно быстро

36 -302.512588577 -1.186549863 0.000000067 0.000008386 5.22
--- The Coupled-Cluster iterations have converged ---

и мы дорываемся до

ORCA NUMERICAL GRADIENT
(16-process run)

которые шли шли шли
<<<Energy at the input geometry>>>
<<<Displacing 7/coordinate 1 (+)>>>
[два десятка однотипных строк выкинуты]
<<<Displacing 8/coordinate 2 (-)>>>

An error has occured in the MDCI module
CALLING COMMAND: /Адрес/Prog/orca_3_0_1_linux_x86-64/orca_mdci CCSD.proc7.mdciinp.tmp>CCSD.proc7.lastcpt
RETURN CODE : 256
Check for MDCI-logfiles
ORCA_NUMFREQ: ORCA finished with an error in the energy calculation
Please check CCSD.lastint.dpl10 for details
Please check CCSD.lastscf.dpl10 for details
Please check CCSD.lastcpt.dpl10 for details

и еще несколько аналогичных ошибок

причем этого CCSD.proc7.lastcpt я потом вообще не нашел
а CCSD.lastcpt.dpl10 for details заканчивается тем что

499 -302.503399735 -1.172381691 0.000000003 0.000025085 33.57
--- The Coupled-Cluster iterations have NOT converged ---

Ну я конечно добавлю ему еще итераций, но у меня была мысль как либо по бороться за сходимость внутри %mdci блока
Но смотрю я мануал 3.0.1 в раздел 6.8.4 Use of the MDCI Module и вижу там фигу. Может ли мне в этом помочь Level shift? Другие опции?

Уходить с SD(T) не поборовшись пока не хочется.

[Гесс]

Насколько невежливо поднять "Max. element of the residual vector for convergence check" с дефолтного STol 1e-5 до 6e-5 ?
Во всех неуспевших сконвергнуться за 500 итераций эта величина достигла 6,2е-5 или менее но снижалась очень вяло (линейно на последних 300 шагах, со скоростью примерно 1e-6 на 100 шагов)

[kbob]

В orca нет аналитического градиента для СС методов, а на вопрос когда будет отвечают, что не скоро:
https://cec.mpg.de/forum/viewtopic.php?f=10&t=116
С другой стороны на форуме поднималась тема про пакеты с аналитическим градиентом в СС, и таких пакетов оказалось предостаточно.
viewtopic.php?f=71&t=100543
Про имплементацию IRC вопрос разработчикам ORCA был задан, но ответа пока нет.
https://cec.mpg.de/forum/viewtopic.php?f=10&t=487
На мой вопрос о том будут ли реализованы методы оптимизации пути реакции NEB и ZTS был дан очень интересный ответ:
https://cec.mpg.de/forum/viewtopic.php?f=10&t=101
было сказано, что Chemshell совместно с ORCA может производить такую оптимизацию, так что в ORCA ничего реализовавать и не надо. Возможно тоже самое относится и к IRC.
С другой стороны в ORCA можно просканировать ППЭ вдоль нормальной координаты (это примерно то-же самое)

6.14.8.8 Potential Energy Surface Scans along Normal Coordinates
The ground and excited state potential energy surfaces can also be mapped as a function of normal coordinates. The
normal mode trajectory run is invoked by the keyword !MTR. In addition several parameters have to be specified in
the block %mtr. The following example illustrates the use:
First you run a frequency job:

P.S. Я лично для расчета IRC пользую Nwchem или Psi4.

[kbob]

Насколько невежливо поднять "Max. element of the residual vector for convergence check" с дефолтного STol 1e-5 до 6e-5 ?
Во всех неуспевших сконвергнуться за 500 итераций эта величина достигла 6,2е-5 или менее но снижалась очень вяло (линейно на последних 300 шагах, со скоростью примерно 1e-6 на 100 шагов)

Если у вас за 500 итераций не сходится то нужно искать причину этого. Может быть у вас базис неадекватен задаче или 1-конфигурационный подход в принципе не применим.