计算化学公社

标题: 求助：ORCA对铀酰配合物进行结构优化基组选择 [打印本页]

作者
Author: suyanxu 时间: 2024-3-10 21:08
标题: 求助：ORCA对铀酰配合物进行结构优化基组选择
目前做的体系原子数在20-60个之间，除铀外还含有C、H、O、N，我目前对铀原子使用SDD赝势，基组选择SARC-ZORA-TZVP，其他原子使用def2-SVP基组。输入文件如下
! PBE0 D3 opt SARC/J noautostart miniprint nopop
%maxcore  2000
%pal nprocs  30 end
%cpcm smd true
  SMDsolvent "water"
end
%basis
  NewGTO U "SARC-ZORA-TZVP" end
  NewECP U "SDD" end
  NewGTO H "def2-SVP" end
  NewGTO N "def2-SVP" end
  NewGTO O "def2-SVP" end
  NewGTO C "def2-SVP" end
end
* xyz 2 1
原子坐标略
*
请问各位老师我这样做合理吗？

作者
Author: wzkchem5 时间: 2024-3-10 21:42
SARC系列基组是配合显式相对论方法（ZORA、DKH2）用的，不能搭配赝势使用。如果计算量允许，尽量用当前基组搭配ZORA。
SMD溶剂模型做结构优化，容易因为数值误差不收敛或有虚频（不过如果实际上收敛了且没有额外虚频，那么结果仍然是可以用的），一般建议用CPCM。
基组有更简单的写法，第一行写def2-SVP，然后删掉所有的NewGTO xxx "def2-SVP" end，和你现在写的效果是一样的。不过你现在的写法也是正确的。
其他应该问题不大

作者
Author: suyanxu 时间: 2024-3-10 22:34

wzkchem5 发表于 2024-3-10 21:42
SARC系列基组是配合显式相对论方法（ZORA、DKH2）用的，不能搭配赝势使用。如果计算量允许，尽量用当前基组 ...

老师请问用当前基组搭配ZORA是指使用ZORA-def2-TZVP基组吗？我试了一下发现报错There are no main basis functions on atom number 43 (U )。如果想使用赝势的话用什么基组比较好？还请老师指点

作者
Author: wzkchem5 时间: 2024-3-10 22:50

suyanxu 发表于 2024-3-10 15:34
老师请问用当前基组搭配ZORA是指使用ZORA-def2-TZVP基组吗？我试了一下发现报错There are no main basi ...

不是，是用SARC-ZORA-TZVP，ZORA关键词是另写的。SARC-ZORA-TZVP、ZORA-def2-TZVP里的ZORA是基组名字的一部分，本身并没有指定用ZORA方法的作用，需要另写"ZORA"这个关键词，说明书里有很多示例。
铀比较重，赝势可能未必可靠，个人觉得还是用ZORA比较好

作者
Author: suyanxu 时间: 2024-3-10 23:19

wzkchem5 发表于 2024-3-10 22:50
不是，是用SARC-ZORA-TZVP，ZORA关键词是另写的。SARC-ZORA-TZVP、ZORA-def2-TZVP里的ZORA是基组名字的一 ...

老师输入文件这样写对吗
! PBE0 D3 ZORA SARC/J opt CPCM(water) noautostart miniprint nopop
%maxcore  2000
%pal nprocs  30 end
%basis
  NewGTO U "SARC-ZORA-TZVP" end
  NewGTO H "def2-SVP" end
  NewGTO N "def2-SVP" end
  NewGTO O "def2-SVP" end
  NewGTO C "def2-SVP" end
end

作者
Author: wzkchem5 时间: 2024-3-10 23:31

suyanxu 发表于 2024-3-10 16:19
老师输入文件这样写对吗
! PBE0 D3 ZORA SARC/J opt CPCM(water) noautostart miniprint nopop
%maxcor ...

我觉得没问题

作者
Author: suyanxu 时间: 2024-3-10 23:41

wzkchem5 发表于 2024-3-10 23:31
我觉得没问题

感谢老师抽出时间解答

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