用高斯优化硫化铜团簇的设置求助

查看数: 5554 | 评论数: 12 | 收藏 Add to favorites 1

关灯 | 提示：支持键盘翻页<-左右->

发布时间: 2021-10-10 23:35

各位老师好，我想计算CuS的激发态性质和实验上对比，感觉用第一性原理软件计算周期结构的计算量非常大，所以想先研究CuS小团簇。我先用MS建了个CuS周期性结构，扩超胞后挖了一个团簇，再用氢饱和。结构图我贴在下 ...

回复 Reply

luffys 发表于 Post on 2021-10-11 23:46:03

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

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

fffff 发表于 Post on 2021-10-11 17:31:25

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

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

fffff 发表于 Post on 2021-10-11 15:56:47

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

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

wzkchem5 发表于 Post on 2021-10-11 15:09:56

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

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

wzkchem5 发表于 Post on 2021-10-11 15:09:34

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

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

sobereva 发表于 Post on 2021-10-11 14:23:06

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

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

luffys 发表于 Post on 2021-10-11 11:20:30

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

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

luffys 发表于 Post on 2021-10-11 11:12:10

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

谢谢老师，我去看一看。

biogon 发表于 Post on 2021-10-11 10:40:36

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

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

fffff 发表于 Post on 2021-10-11 06:36:42

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

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

SharedScreenshot.jpg (24.95 KB, 下载次数 Times of downloads: 3)

sobereva 发表于 Post on 2021-10-11 04:17:38

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

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

wzkchem5 发表于 Post on 2021-10-10 23:53:41

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

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