计算化学公社

标题: 用高斯优化硫化铜团簇的设置求助 [打印本页]

作者
Author: luffys 时间: 2021-10-10 23:35
标题: 用高斯优化硫化铜团簇的设置求助
各位老师好，我想计算CuS的激发态性质和实验上对比，感觉用第一性原理软件计算周期结构的计算量非常大，所以想先研究CuS小团簇。
我先用MS建了个CuS周期性结构，扩超胞后挖了一个团簇，再用氢饱和。结构图我贴在下面了（图分别是用MS显示的团簇示意图；以及GV的示意图，发光的是固定的原子）

我想用高斯先优化下面的结构，下面是具体的输入，请问这样的设置合理吗？赝势和基组的选择是依据文献中计算Cd[size=6.029pt]10Se[size=6.029pt]6S[size=6.029pt]14团簇的设置。
[size=6.029pt]

先谢谢各位老师了！

%chk=1.chk
#p opt b3lyp/lanl2dz  SCF=(Fermi) int=(ultrafine,acc2e=12)

CuS_Cluster

0 1
Cu             0 1.09600000 1.89800000 6.46000000
Cu             0 1.09600000 1.89800000 1.76000000
Cu             0 2.19200000 0.00000000 4.11000000
S             0 0.00000000 0.00000000 7.16700000
S             0 0.00000000 0.00000000 1.05400000
S             -1 0.00000000 0.00000000 9.27400000
S             0 1.09600000 1.89800000 4.11000000
S             0 0.00000000 3.79700000 1.05400000
Cu             0 4.38400000 0.00000000 1.76000000
S             -1 3.28800000 -1.89800000 1.05400000
S             0 4.38400000 0.00000000 4.11000000
S             0 3.28800000 1.89800000 1.05400000
S             -1 6.57600000 0.00000000 1.05400000
Cu             0 -2.19200000 -3.79700000 1.76000000
S             -1 -3.28800000 -5.69500000 1.05400000
S             0 -2.19200000 -3.79700000 4.11000000
Cu             0 -2.19200000 0.00000000 6.46000000
Cu             0 -2.19200000 0.00000000 1.76000000
Cu             0 -1.09600000 -1.89800000 4.11000000
S             -1 -3.28800000 -1.89800000 1.05400000
S             0 -2.19200000 0.00000000 4.11000000
Cu             0 1.09600000 -1.89800000 6.46000000
Cu             0 1.09600000 -1.89800000 1.76000000
S             -1 0.00000000 -3.79700000 1.05400000
S             0 1.09600000 -1.89800000 4.11000000
Cu             0 -2.19200000 3.79700000 1.76000000
Cu             0 -1.09600000 1.89800000 4.11000000
S             -1 -3.28800000 1.89800000 1.05400000
S             0 -2.19200000 3.79700000 4.11000000
S             -1 -3.28800000 5.69500000 1.05400000
H             0 0.94000000 0.94000000 9.74400000
H             0 7.98600000 0.00000000 1.05400000
H             0 -1.87800000 -5.69500000 1.05400000
H             0 -1.87800000 5.69500000 1.05400000

作者
Author: wzkchem5 时间: 2021-10-10 23:53
CuS的结构和CdS不一样，挺复杂的，还有二硫键，参见https://en.wikipedia.org/wiki/Copper_monosulfide
此外，B3LYP算过渡金属一般不如PBE0可靠。重新建模以后建议改用PBE0算

作者
Author: sobereva 时间: 2021-10-11 04:17
又不考虑背景电荷又不考虑capped ECP或赝氢之类处理，直接这么算，而且团簇尺寸又不大，合理性有限。建议至少对比一下团簇取再大一圈时候的结果，以证明当前大小团簇能说明问题。

写int=(ultrafine,acc2e=12) 完全多余，本身B3LYP对积分格点要求又不算高，基组又没带弥散函数

作者
Author: fffff 时间: 2021-10-11 06:36

wzkchem5 发表于 2021-10-10 23:53
CuS的结构和CdS不一样，挺复杂的，还有二硫键，参见https://en.wikipedia.org/wiki/Copper_monosulfide
此 ...

老师，请问一下这么多自旋中心不用考虑多参考方法吗？还是说几何优化用单参考多参考影响不大呢？
我看Grimmr提出FOD的原文里有两个Cu，静态相关就比较强了

作者
Author: biogon 时间: 2021-10-11 10:40

fffff 发表于 2021-10-11 06:36
老师，请问一下这么多自旋中心不用考虑多参考方法吗？还是说几何优化用单参考多参考影响不大呢？
我看Gr ...

要考虑啊，orca的casscf练习里面就有一个和这个很相似的例子

作者
Author: luffys 时间: 2021-10-11 11:12

wzkchem5 发表于 2021-10-10 23:53
CuS的结构和CdS不一样，挺复杂的，还有二硫键，参见https://en.wikipedia.org/wiki/Copper_monosulfide
此 ...

谢谢老师，我去看一看。

作者
Author: luffys 时间: 2021-10-11 11:20

sobereva 发表于 2021-10-11 04:17
又不考虑背景电荷又不考虑capped ECP或赝氢之类处理，直接这么算，而且团簇尺寸又不大，合理性有限。建议至 ...

谢谢社长，结构优化第一轮收敛就失败了，两个不同尺寸的团簇都没撑过优化的第一轮，我正在尝试您那篇解决SCF不收敛的博文里提到的方法。。。

作者
Author: sobereva 时间: 2021-10-11 14:23

fffff 发表于 2021-10-11 06:36
老师，请问一下这么多自旋中心不用考虑多参考方法吗？还是说几何优化用单参考多参考影响不大呢？
我看Gr ...

这都得用多参考的话，第一性原理领域DFT算的一堆体系都根本没法做了
没这么夸张

作者
Author: wzkchem5 时间: 2021-10-11 15:09

fffff 发表于 2021-10-10 23:36
老师，请问一下这么多自旋中心不用考虑多参考方法吗？还是说几何优化用单参考多参考影响不大呢？
我看Gr ...

能跑得动的话需要，但是这么多Cu肯定跑不动，所以只能用DFT
最多只能是用双核的团簇做一个多参考计算，再做一系列DFT计算，看哪个泛函合适，再用这个泛函算实际的团簇

作者
Author: wzkchem5 时间: 2021-10-11 15:09

luffys 发表于 2021-10-11 04:20
谢谢社长，结构优化第一轮收敛就失败了，两个不同尺寸的团簇都没撑过优化的第一轮，我正在尝试您那篇解决 ...

先解决结构问题，你结构是错的。

作者
Author: fffff 时间: 2021-10-11 15:56

wzkchem5 发表于 2021-10-11 15:09
能跑得动的话需要，但是这么多Cu肯定跑不动，所以只能用DFT
最多只能是用双核的团簇做一个多参考计算， ...

噢噢，是以多参考为标准做一个benchamrk。谢谢老师~

作者
Author: fffff 时间: 2021-10-11 17:31

sobereva 发表于 2021-10-11 14:23
这都得用多参考的话，第一性原理领域DFT算的一堆体系都根本没法做了
没这么夸张

好的，谢谢sob老师，我不太了解多参考方法

作者
Author: luffys 时间: 2021-10-11 23:46

wzkchem5 发表于 2021-10-11 15:09
先解决结构问题，你结构是错的。

嗯嗯好的，一开始没往复杂处想，觉得先优化一个小团簇看看结构，想得太简单了，看来这个真的是结构问题，我去调整一下结构再优化。

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