ORCA 5 ,Yes

wzkchem5 · 发表于 Post on 2021-7-8 22:39:34

hdhxx123 发表于 2021-7-8 15:25
老师，想请问一下Tcut和Thresh在orca4.2.1中的设置与scf级别的设置是怎么样的呢？是orca5中有改变吗？我 ...

应该没有大变化。你自己用orca跑一下就知道了，输入文件不写miniprint，就会打出tcut和thresh值

hdhxx123 · 发表于 Post on 2021-7-8 22:47:35

wzkchem5 发表于 2021-7-8 22:39
应该没有大变化。你自己用orca跑一下就知道了，输入文件不写miniprint，就会打出tcut和thresh值

好的，谢谢老师

hdhxx123 · 发表于 Post on 2021-7-8 23:30:02

wzkchem5 发表于 2021-7-8 20:53
注意一个问题：你把thresh和tcut设得比默认值还松（使用verytightscf时， thresh默认为1e-12，tcut默认为 ...

我之前不清楚这个事儿然后orca4.2.1手册里也没具体写，也一直带着miniprint，结果以为这个设置是比较好的设置。跑了才发现像做opt时verytightscf默认就是这么严，其实这两个值一般就应该是Thresh 1e-12 TCut 1e-14；如果不收敛应该设的更小，谢谢老师提醒了！

wzkchem5 · 发表于 Post on 2021-7-9 00:56:14

hdhxx123 发表于 2021-7-8 16:30
我之前不清楚这个事儿然后orca4.2.1手册里也没具体写，也一直带着miniprint，结果以为这个设置是比较好的 ...

个人建议平时不要带miniprint和nopop。。。如果计算成功，你基本上可以用grep读你需要的结果，程序多输出一些内容也完全不影响；如果计算失败，那么程序多打出来一些信息反而有助于你找问题在哪。尤其是population analysis，有时候计算收敛到了错误的波函数上，看布居分析结果就能看出来，如果没有布居分析结果，就得导出molden、下到本地然后用可视化软件看，会麻烦不少

biogon · 发表于 Post on 2021-7-10 15:40:58

wzkchem5 发表于 2021-7-9 00:56
个人建议平时不要带miniprint和nopop。。。如果计算成功，你基本上可以用grep读你需要的结果，程序多输出 ...

orca未来有支持更高级的相对论哈密顿比如DKH4 X2C之类的计划吗？
哦对了现在有没有能在相对论哈密顿下做波函数优化的程序，bdf发现似乎是没有？

wzkchem5 · 发表于 Post on 2021-7-10 16:06:46

biogon 发表于 2021-7-10 08:40
orca未来有支持更高级的相对论哈密顿比如DKH4 X2C之类的计划吗？
哦对了现在有没有能在相对论哈密顿下做 ...

你指的波函数优化是stable=opt吗？BDF本来有这个功能，但是有bug，而且只能用MO基算法，还没有实现AO基算法。
orca有做sf-X2C的打算，但是近期做的概率不太大

hdhxx123 · 发表于 Post on 2021-7-11 19:48:45

wzkchem5 发表于 2021-7-8 20:53
注意一个问题：你把thresh和tcut设得比默认值还松（使用verytightscf时， thresh默认为1e-12，tcut默认为 ...

老师，我测试了一下，TRAH的效果确实在Tb体系来说远不如纯DIIS。

wzkchem5 · 发表于 Post on 2021-7-12 15:57:56

hdhxx123 发表于 2021-7-11 12:48
老师，我测试了一下，TRAH的效果确实在Tb体系来说远不如纯DIIS。

TRAH的迭代次数多于DIIS是正常的，但是至于TRAH对于这个例子为什么不收敛，我还得看一下

wzkchem5 · 发表于 Post on 2021-7-12 20:55:35

hdhxx123 发表于 2021-7-11 12:48
老师，我测试了一下，TRAH的效果确实在Tb体系来说远不如纯DIIS。

这个算例TRAH不收敛是因为格点不够大。而且很有意思的是连DefGrid3都不够大，需要用RefGrid才能收敛。DIIS用默认格点能收敛应该只是巧合，你也看到了DIIS收敛也很困难。用RefGrid的情况下，TRAH用了265步收敛，收敛能量-822.886863903289。
这个体系之所以对格点要求这么高，应该是因为f轨道的简并性导致的。Tb3+是4f8，自旋向下的那个f电子可以占据任一个f轨道，理论上占据哪个f轨道得到的能量都是一样的才对，但是如果格点的角向部分不够大，这些态的能量会有细微的区别，导致程序不知道该收敛到哪里去。

hdhxx123 · 发表于 Post on 2021-7-12 21:31:43

本帖最后由 hdhxx123 于 2021-7-12 21:53 编辑

wzkchem5 发表于 2021-7-12 20:55
这个算例TRAH不收敛是因为格点不够大。而且很有意思的是连DefGrid3都不够大，需要用RefGrid才能收敛。DII ...

老师，那diis在默认格点下收敛的波函数与RefGrid下TRAH收敛的波函数有明显区别吗？耗时相差多少？（应该不是巧合，我其他十五种稀土离子在默认格点下diis都能正常收敛）
另：diis同TRAH在相同格点下收敛的能量相差多少。

hdhxx123 · 发表于 Post on 2021-7-12 21:54:55

wzkchem5 发表于 2021-7-12 20:55
这个算例TRAH不收敛是因为格点不够大。而且很有意思的是连DefGrid3都不够大，需要用RefGrid才能收敛。DII ...

另：发现一个问题，在关键词里同时写noTRAH和TRAH它并不会报错，而计算时会没有[F,P]和Damp两行并且出现收敛会问题，建议优化

wzkchem5 · 发表于 Post on 2021-7-12 22:05:20

hdhxx123 发表于 2021-7-12 14:31
老师，那diis在默认格点下收敛的波函数与RefGrid下TRAH收敛的波函数有明显区别吗？耗时相差多少？（应该 ...

你可以算一下看看，我没用DIIS算。导致这个问题的关键在于f8的电子构型对格点特别敏感，其他f电子构型不一定有这个问题，所以其他稀土能收敛，不代表Tb3+能收敛不是巧合。

wzkchem5 · 发表于 Post on 2021-7-12 22:25:58

hdhxx123 发表于 2021-7-12 14:54
另：发现一个问题，在关键词里同时写noTRAH和TRAH它并不会报错，而计算时会没有[F,P]和Damp两行并且出现 ...

感谢反馈，检查了一下代码，这是因为TRAH会关掉DIIS，而NoTRAH不会打开DIIS。
我改了一下代码，从ORCA 5.0.1开始，同时写TRAH和NoTRAH会报错

biogon · 发表于 Post on 2021-7-12 22:52:14

wzkchem5 发表于 2021-7-10 16:06
你指的波函数优化是stable=opt吗？BDF本来有这个功能，但是有bug，而且只能用MO基算法，还没有实现AO基算 ...

那我用ano-rcc基组如果想做类似stable=opt之类的东西是不是就没辙了，睾斯完全算不动

wzkchem5 · 发表于 Post on 2021-7-12 23:07:02

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

可以用小一些的基组做stable=opt，再读取这个波函数做ANO-RCC的计算，然后做stable。。。
BDF的stable功能也有问题，但是可以用spin-conserving U-TDDFT没有负、虚、复激发能这个标准来代替（这个条件严格等价于没有任何的instability），多费大概一倍计算量而已

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

[ORCA] ORCA 5 ,Yes

评分 Rate

评分 Rate

浏览过的版块