计算化学公社

标题: orca收敛不出一个局域最小值 [打印本页]

作者
Author: dlldwtz 时间: 2024-7-30 10:04
标题: orca收敛不出一个局域最小值
对于图中的这个结构，试过不同的scf求解器、初猜、或者rohf都不收敛，orca6.0或者5.0.3/4都没成功，但是molpro用hf,so-sci可以收敛

作者
Author: zjxitcc 时间: 2024-7-30 11:28
本帖最后由 zjxitcc 于 2024-7-30 11:34 编辑

既然讨论的是SCF收敛性问题，不需要进行UCCSD(T)-F12计算，直接进行UHF计算即可。你可以晒出你能收敛的UHF结果（不收敛的就不要展示了），我可以用MOKIT给出非常简单、迅速、便捷、不需要很多人为尝试的解决办法（当然你自己用MOKIT也可以实现，只不过有的用户没用过MOKIT，思想里给自己设置了一个壁垒，因此我来展示例子和结果更为直观和震撼）。
PS：
（1）给出Molpro结果时应当指明是ROHF能量还是UHF能量，因为如果你直接做UCCSD(T)-F12计算，Molpro不支持UHF-UCCSD(T)，你只能进行ROHF计算，而这与UHF能量必然不相等，ROHF收敛性也与UHF是两个问题，没法混在一起比较。
（2）单参考计算不需要PModel关键词。
（3）你用的基组过于复杂，如果讨论UHF，我们可以先局限于cc-pVTZ和cc-pVTZ-PP，简单方便。
这三点的意思是这里面蕴含太多的其他因素，我们需要先控制变量。

作者
Author: dlldwtz 时间: 2024-7-30 14:15

zjxitcc 发表于 2024-7-30 11:28
既然讨论的是SCF收敛性问题，不需要进行UCCSD(T)-F12计算，直接进行UHF计算即可。你可以晒出你能收敛的UHF ...

好的，我把输入文件简化一下试一下。请问老师，orca如何将rohf轨道转为uhf轨道从而进行uccsd计算呢？

作者
Author: wzkchem5 时间: 2024-7-30 15:08
https://www.faccts.de/docs/orca/ ... tions-to-an-element

Note that the NewGTO keyword does not change the ECP for the given element - you must use NewECP or DelECP (see section Advanced Specification of Effective Core Potentials).

作者
Author: dlldwtz 时间: 2024-7-30 15:25

wzkchem5 发表于 2024-7-30 15:08
https://www.faccts.de/docs/orca/6.0/manual/contents/detailed/basisset.html#assigning-or-adding-basis ...

原来如此，谢谢老师！

作者
Author: dlldwtz 时间: 2024-7-30 15:58

wzkchem5 发表于 2024-7-30 15:08
https://www.faccts.de/docs/orca/6.0/manual/contents/detailed/basisset.html#assigning-or-adding-basis ...

老师您好，现在uhf可以收敛了，但是出现以下报错，请问是为什么呢？多个单点计算均出现这个问题

Initializing the integral package       ... done
Warning: reference  - re-canonicalizations have been set to INT 1 VIRT 1
Fock densities were made
Making Fock matrix and reference energy:

--------------------------
UNRESTRICTED FOCK OPERATOR
--------------------------

Recanonicalizing the internal orbitals
Recanonicalizing the virtual orbitals
Storing new orbitals

ORCA finished by error termination in MDCI
Calling Command: ../orca_6_0_0_shared_openmpi416_avx2/orca_mdci pw.mdciinp.tmp
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run

作者
Author: wzkchem5 时间: 2024-7-30 16:36

dlldwtz 发表于 2024-7-30 08:58
老师您好，现在uhf可以收敛了，但是出现以下报错，请问是为什么呢？多个单点计算均出现这个问题

Initi ...

检查默认的内存是否不足。如果不是内存问题，可能是bug，可以在orca论坛上报一下bug

作者
Author: dlldwtz 时间: 2024-7-30 16:46

wzkchem5 发表于 2024-7-30 16:36
检查默认的内存是否不足。如果不是内存问题，可能是bug，可以在orca论坛上报一下bug

老师好，经过我测试发现，将以下输入文件

! uhf ccsd(t)-f12/ri cc-pvdz-f12 cc-pvdz-f12-optri cc-pvtz/c verytightscf notrah

%basis
ecp "sk-mcdhf-rsc"
newgto br "cc-pvdz-pp-f12" end
newcabsgto br "cc-pvdz-pp-f12-optri" end
newauxcgto br "cc-pvtz-pp/c" end
end

%scf
convforced 0
stabperform true
end

* gzmt 0 2
h
o 1 0.9606
h 2 0.9606 1 104.729
br 3 3.0973 2 60.85 1 105.244
*

中的Br原子改为Cl原子并删除%basis块后可以正常结束uhf ccsd(t)-f12/ri的计算。我的赝势部分是否仍格式有问题？

作者
Author: wzkchem5 时间: 2024-7-30 17:15

dlldwtz 发表于 2024-7-30 09:46
老师好，经过我测试发现，将以下输入文件

! uhf ccsd(t)-f12/ri cc-pvdz-f12 cc-pvdz-f12-optri cc-pv ...

我觉得还是值得检查一下是不是默认的maxcore不够，毕竟Br的电子比Cl多，可能空轨道也多

作者
Author: zjxitcc 时间: 2024-7-30 19:57

dlldwtz 发表于 2024-7-30 14:15
好的，我把输入文件简化一下试一下。请问老师，orca如何将rohf轨道转为uhf轨道从而进行uccsd计算呢？

不是把ROHF轨道转为UHF轨道，是你要做UHF计算，不是做ROHF计算。

作者
Author: dlldwtz 时间: 2024-7-30 21:59

wzkchem5 发表于 2024-7-30 17:15
我觉得还是值得检查一下是不是默认的maxcore不够，毕竟Br的电子比Cl多，可能空轨道也多

老师好，我测试发现应该是6.0.0的bug，因为5.0.4可以正常运行。
另外想请教一个问题，对于含Br原子的体系，f12基组和cabs基组使用了带-pp的赝势基组，ri加速使用的auxc基组是否一定也要带-pp？至少程序不认为newauxcgto br "cc-pvtz/c"与之前的赝势基组冲突。我问这个问题的原因是ri-jk加速所支持的基组中不包含cc-pvtz-pp/jk而只有cc-pvtz/jk。谢谢老师！

作者
Author: wzkchem5 时间: 2024-7-30 23:48

dlldwtz 发表于 2024-7-30 14:59
老师好，我测试发现应该是6.0.0的bug，因为5.0.4可以正常运行。
另外想请教一个问题，对于含Br原子的体 ...

应该确实是bug，已经和相关负责人反映了，正在排查。
关于auxc基组不带pp是不是一定不行，我不是太清楚，但是带pp肯定更保险。程序不能保证把所有不合理的基组-辅助基组搭配都检查出来，因为基组搭配从合理到不合理是一个连续谱

欢迎光临计算化学公社 (http://bbs.keinsci.com/)