ORCA 5 ,Yes

hebrewsnabla · 发表于 Post on 2021-7-12 23:10:10

biogon 发表于 2021-7-12 22:52
那我用ano-rcc基组如果想做类似stable=opt之类的东西是不是就没辙了，睾斯完全算不动

pyscf好像有

喵星大佬 · 发表于 Post on 2021-7-12 23:17:38

biogon 发表于 2021-7-12 22:52
那我用ano-rcc基组如果想做类似stable=opt之类的东西是不是就没辙了，睾斯完全算不动

在openmolcas里面算好用mokit传给高斯做stable，那玩意算ano之类的基组应该是没有更好的了

biogon · 发表于 Post on 2021-7-13 08:44:32

本帖最后由 biogon 于 2021-7-13 08:48 编辑

喵星大佬发表于 2021-7-12 23:17
在openmolcas里面算好用mokit传给高斯做stable，那玩意算ano之类的基组应该是没有更好的了

ano基组即便是2zeta的算二三十个原子的体系单点高斯都算不动

biogon · 发表于 Post on 2021-7-13 08:46:06

hebrewsnabla 发表于 2021-7-12 23:10
pyscf好像有

pyscf有的话那也不错

biogon · 发表于 Post on 2021-7-13 08:47:56

wzkchem5 发表于 2021-7-12 23:07
可以用小一些的基组做stable=opt，再读取这个波函数做ANO-RCC的计算，然后做stable。。。
BDF的stable功 ...

高斯的问题是在于ano一般也单点就能算的动几个原子，更别说干别的事了

hdhxx123 · 发表于 Post on 2021-7-13 09:08:34

wzkchem5 发表于 2021-7-12 22:05
你可以算一下看看，我没用DIIS算。导致这个问题的关键在于f8的电子构型对格点特别敏感，其他f电子构型不 ...

测试结果：1.如果是这种猜测，则Eu即f6的电子构型也应该有这样的情况，但实际上DIIS和TRAH均在默认格点下正常迅速收敛。不知道应当如何证明格点敏感的猜测？
2.
实测：

Eu：默认格点
TRAH：-708.784257738945 1 minutes 5 seconds 322 msec
DIIS：-708.784257738238 38 seconds 320 msec
Tb：refgrid
TRAH:-822.886863903288 20 minutes 22 seconds 210 msec
DIIS:-822.886863902504 24 minutes 32 seconds 67 msec

复制代码

处理器：Xeon E3-1505M v5 四核。相同格点下应该都是收敛到相同解

wzkchem5 · 发表于 Post on 2021-7-13 15:35:16

hdhxx123 发表于 2021-7-13 02:08
测试结果：1.如果是这种猜测，则Eu即f6的电子构型也应该有这样的情况，但实际上DIIS和TRAH均在默认格点下 ...

收敛的难易程度方面，不能简单看particle-hole symmetry。举个极端的例子，假设你算Ce3+，除了4f以外所有电子都冻结，且初猜是Ce4+的轨道，也就是球对称的，那么你只迭代一步就收敛了，哪怕格点破坏了4f的简并性也是如此。但是如果你算Yb3+，初猜是Yb2+，那么需要迭代很多步，因为你的基组里的f基函数不止一套。另外4f8有一个特点，就是其中双占的那个f轨道很高，容易和5d混合，加剧对称性破缺，这也是4f6不具备的特点。
或许你可以试一下其他价态的稀土，看看4f6是不是普遍比4f8好收敛

hdhxx123 · 发表于 Post on 2021-7-13 16:10:41

wzkchem5 发表于 2021-7-13 15:35
收敛的难易程度方面，不能简单看particle-hole symmetry。举个极端的例子，假设你算Ce3+，除了4f以外所有 ...

学习了，谢谢老师

hdhxx123 · 发表于 Post on 2021-7-13 22:19:51

wzkchem5 发表于 2021-7-13 15:35
收敛的难易程度方面，不能简单看particle-hole symmetry。举个极端的例子，假设你算Ce3+，除了4f以外所有 ...

老师您这里说的双占f轨道能量很高，我在高斯里计算得到的Tb.wfn用Multiwfn打开发现有两个轨道能量是相对低的，但是其中一个要低得多，不知道是什么原因？是双占轨道的alpha和beta电子能量不一样吗？想咨询一下老师

wzkchem5 · 发表于 Post on 2021-7-14 00:22:07

hdhxx123 发表于 2021-7-13 15:19
老师您这里说的双占f轨道能量很高，我在高斯里计算得到的Tb.wfn用Multiwfn打开发现有两个轨道能量是相对 ...

你看的是哪两个轨道？
UKS计算，alpha轨道和beta轨道能量不一样。像这个体系，因为有大量的alpha电子，自旋极化很严重，所以alpha轨道和对应的beta轨道的能量完全没有可比性

hdhxx123 · 发表于 Post on 2021-7-14 07:30:26

wzkchem5 发表于 2021-7-14 00:22
你看的是哪两个轨道？
UKS计算，alpha轨道和beta轨道能量不一样。像这个体系，因为有大量的alpha电子， ...

看的是alpha和beta轨道的最低轨道能量，学习了，谢谢老师

wzkchem5 · 发表于 Post on 2021-7-24 02:11:15

ORCA 5.0.1已发布！
https://orcaforum.kofo.mpg.de/viewtopic.php?f=8&t=7785
修复了ORCA 5.0.0的相当一部分已知bug（https://orcaforum.kofo.mpg.de/viewforum.php?f=51），欢迎试用。视接下来用户报告的bug数量而定，预期会择机发布ORCA 5.0.2

sobereva · 发表于 Post on 2021-7-25 01:07:15

前面我提到的MD会崩溃的问题实测5.0.1已经解决了，而且经测试也没有ORCA 4.2跑AIMD到了2000步左右之后每一步越来越慢的问题

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] ORCA 5 ,Yes

评分 Rate

评分 Rate

评分 Rate