计算化学公社

标题: 求助orca优化Cu16Pd1簇结构输入文件写法 [打印本页]

作者
Author: Huyang 时间: 2023-4-22 21:22
标题: 求助orca优化Cu16Pd1簇结构输入文件写法
orca新手，用的orca版本5.0.4，请问一下大家，如果优化339个原子的Cu16Pd簇用下面的

! B97-3c opt noautostart miniprint nopop
%maxcore  5000
%pal nprocs 13 end
* xyz 0 1

坐标
*

还是
! B3LYP D3 def2-TZVP(-f) def2/J RIJCOSX opt noautostart miniprint nopop
%maxcore  x
%pal nprocs y end
* xyz 0 1

坐标
*

比较好？还有就是def2应该是对Pd自动加了赝式基组把，如果我对铜也加赝式，下面的写法对吗？

! B3LYP D3 def2-TZVP(-f) def2/J RIJCOSX opt noautostart miniprint nopop
%maxcore  x
%pal nprocs y end
%basis
ECP{LanL=[Cu]} end
end
* xyz 0 1

坐标
*

作者
Author: wzkchem5 时间: 2023-4-22 23:35
按http://sobereva.com/272选泛函
对铜用赝势的话，需要用NewGTO和NewECP关键字，具体写法在手册里搜

作者
Author: dnlx 时间: 2023-4-22 23:49
体系太大，应该从b97-3c 开始试，而且怀疑超算几天能算出来。
另外按Grimme 组去年DFT protocol文章，B3LYP对于b97-3c 的优化结果并没有优势

作者
Author: sobereva 时间: 2023-4-23 08:02
B3LYP算这种体系没任何优点
ORCA靠RIJ，用B97-3c纯泛函计算对这种大体系远比用RIJCOSX的B3LYP快，结果也不比用B3LYP差，输入文件写起来还省事

作者
Author: Huyang 时间: 2023-4-23 11:43

wzkchem5 发表于 2023-4-22 23:35
按http://sobereva.com/272选泛函
对铜用赝势的话，需要用NewGTO和NewECP关键字，具体写法在手册里搜

多谢老师，不知道为啥orca官网上不去，特殊手段也不行。我看到有Gaussian Vs ORCA里面写到ECP时有以下关键词，然后我就用了第一个。
ECP{LanL=[Rh]}[/b]
ECP{LanL=[Rb-Ag]}
ECP{LanL=[Rh],def2-TZVP/J}
ECP{LanL=[Rh],LanLDZ,def2-TZVP/J}

我看到要对某些原子指定基组的时候用到%basis NewGTO 元素 “基组” end这种写法。
ECP这里我不太懂，下面这样写可以吗？
! B3LYP D3 def2-TZVP(-f) def2/J RIJCOSX opt noautostart miniprint nopop
%maxcore x
%pal nprocs y end
%basis
NewECP Cu "lanl2" end (还需要在这下面加一行写NewGTO吗？）
NewGTO Cu "lanl2dz" end(这一行需要吗？）
end
* xyz 0 1

坐标
*

作者
Author: Huyang 时间: 2023-4-23 11:44

dnlx 发表于 2023-4-22 23:49
体系太大，应该从b97-3c 开始试，而且怀疑超算几天能算出来。
另外按Grimme 组去年DFT protocol文章，B3LY ...

好的，多谢老师，我用B97-3C跑起来了。条件限制只能16核，每个核分配大概6G，不知道能不能跑出来
作者
Author: Huyang 时间: 2023-4-23 11:45

sobereva 发表于 2023-4-23 08:02
B3LYP算这种体系没任何优点
ORCA靠RIJ，用B97-3c纯泛函计算对这种大体系远比用RIJCOSX的B3LYP快，结果也不 ...

多谢sob老师，用B97-3c跑起来了。
作者
Author: wzkchem5 时间: 2023-4-23 15:07

Huyang 发表于 2023-4-23 04:43
多谢老师，不知道为啥orca官网上不去，特殊手段也不行。我看到有Gaussian Vs ORCA里面写到ECP时有以下关 ...

需要NewGTO
作者
Author: Huyang 时间: 2023-5-9 20:02
感谢大家，我把配体的一些柔性基团用H代替，用B97-3C，36核，每个核5G，跑了2天12h正常收敛。

欢迎光临计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3