计算化学公社

标题: 关于orca计算锕系元素的基组原则以及写法的求助，希望得到检查 [打印本页]

作者
Author: 奶啵啵 时间: 2023-4-13 10:35
标题: 关于orca计算锕系元素的基组原则以及写法的求助，希望得到检查
本帖最后由奶啵啵于 2023-4-13 17:31 编辑

已解决利用orca计算铀单点能，输入文件如下所示希望得到检查，计算的特别慢，是不是我的输入有问题
! wB97M-V def2-QZVPP def2/J RIJCOSX strongSCF noautostart miniprint nopop
%maxcore 1500
%pal nprocs 12 end
%cpcm
smd true
SMDsolvent "water"
end
%basis
NewGTO U

（省略）（SARC-ZORA)
end
NewAuxJGTO U "autoaux" end

end
* xyz 2 1

（省略结构信息）
*

作者
Author: 七尺贱 时间: 2023-4-13 11:09
有多少个原子呐，def2QZVPP本来就是超大基组，算的就慢，而且你每个核心给的内存太少了，一核不得给好几个G甚至十几个

作者
Author: 奶啵啵 时间: 2023-4-13 11:13

七尺贱发表于 2023-4-13 11:09
有多少个原子呐，def2QZVPP本来就是超大基组，算的就慢，而且你每个核心给的内存太少了，一核不得给好几个G ...

3个原子，哈哈哈啊哈，就是很慢，哎，然后out经常卡在一个地方不懂

作者
Author: wxsxjj 时间: 2023-4-13 13:29

奶啵啵发表于 2023-4-13 11:13
3个原子，哈哈哈啊哈，就是很慢，哎，然后out经常卡在一个地方不懂

不要小瞧计算量啊

作者
Author: Satoru 时间: 2023-4-13 13:45
cc-pVQZ-DK是给DKH2相对论计算设计的基组，不适合和赝势一同使用

作者
Author: 奶啵啵 时间: 2023-4-13 15:22

Satoru 发表于 2023-4-13 13:45
cc-pVQZ-DK是给DKH2相对论计算设计的基组，不适合和赝势一同使用

嗯。。。。。。就，想算铀，还想和QZVPP一个级别，嗯。。。。。推荐一下？

作者
Author: 奶啵啵 时间: 2023-4-13 17:31

Satoru 发表于 2023-4-13 13:45
cc-pVQZ-DK是给DKH2相对论计算设计的基组，不适合和赝势一同使用

已解决

作者
Author: sobereva 时间: 2023-4-13 17:45

奶啵啵发表于 2023-4-13 15:22
嗯。。。。。。就，想算铀，还想和QZVPP一个级别，嗯。。。。。推荐一下？

普通泛函对U用def-TZVPP足矣

作者
Author: 奶啵啵 时间: 2023-4-14 08:24

sobereva 发表于 2023-4-13 17:45
普通泛函对U用def-TZVPP足矣

好的老师，我刚才，用SARC-ZORA-TZVPP试了一下，（参考了文献）得到的能量数据相当离谱

作者
Author: 奶啵啵 时间: 2023-4-14 08:45

sobereva 发表于 2023-4-13 17:45
普通泛函对U用def-TZVPP足矣

刚才查了下手册第9章,tzvpp 是只到36号元素么

作者
Author: sobereva 时间: 2023-4-15 00:26

奶啵啵发表于 2023-4-14 08:45
刚才查了下手册第9章,tzvpp 是只到36号元素么

完整的def系列的基组定义去https://www.cosmologic-services.de/basis-sets/basissets.php看

作者
Author: sobereva 时间: 2023-4-15 00:26

奶啵啵发表于 2023-4-14 08:24
好的老师，我刚才，用SARC-ZORA-TZVPP试了一下，（参考了文献）得到的能量数据相当离谱

要对全电子相对论计算有基本的概念，看比如
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

作者
Author: 奶啵啵 时间: 2023-4-16 19:53

sobereva 发表于 2023-4-15 00:26
完整的def系列的基组定义去https://www.cosmologic-services.de/basis-sets/basissets.php看

好的好的

作者
Author: 奶啵啵 时间: 2023-4-16 20:08

sobereva 发表于 2023-4-15 00:26
要对全电子相对论计算有基本的概念，看比如
使用Gaussian做镧系金属配合物的量子化学计算
http://sober ...

补一个小问题，看到def系列基组的输出文件只有Gaussian、Turbomole、molpro、Dalton这四种，Orca是不支持么？

作者
Author: zjxitcc 时间: 2023-4-16 20:16

奶啵啵发表于 2023-4-16 20:08
补一个小问题，看到def系列基组的输出文件只有Gaussian、Turbomole、molpro、Dalton这四种，Orca是不支持 ...

“def系列基组的输出文件”是什么意思？

作者
Author: sobereva 时间: 2023-4-16 23:31

奶啵啵发表于 2023-4-16 20:08
补一个小问题，看到def系列基组的输出文件只有Gaussian、Turbomole、molpro、Dalton这四种，Orca是不支持 ...

你自己手动转成ORCA的格式就完了
任何基于Gaussian函数的基组，在ORCA里都能自定义使用

作者
Author: 奶啵啵 时间: 2023-4-17 08:54
本帖最后由奶啵啵于 2023-4-17 08:56 编辑

zjxitcc 发表于 2023-4-16 20:16
“def系列基组的输出文件”是什么意思？

不好意思，没描述清楚，就是在网站上没看到orca格式的输出文件

作者
Author: 奶啵啵 时间: 2023-4-17 17:15

sobereva 发表于 2023-4-16 23:31
你自己手动转成ORCA的格式就完了
任何基于Gaussian函数的基组，在ORCA里都能自定义使用

在补一个小问题
基组定义NewGTO U "cc-pVQZ-PP" 、（根据文献）
赝势选择Stuttgart赝势NewECP U "SDD"
最后的单点能看起来很可靠（和之前结构优化数据差距不大，不像以前那么离谱）
这样合理么？

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