计算化学公社

标题: 想算carbonmonoxide偶极矩，静电势。CCSD-cc-PV5Z基组已经一天了 [打印本页]

作者
Author: Zhangchj 时间: 2023-10-18 08:57
标题: 想算carbonmonoxide偶极矩，静电势。CCSD-cc-PV5Z基组已经一天了
本帖最后由 Zhangchj 于 2023-10-18 09:00 编辑

看到其他人的文章中采用这个基组，所以选了相同基组。基组来源于BSE网站，内置的基组大概也就十几秒就算完了，

作者
Author: zjxitcc 时间: 2023-10-18 13:12
本帖最后由 zjxitcc 于 2023-10-18 13:25 编辑

这个计算量本来就很大，你得对5Z基组有点概念。不过，我还是建议先看看输入文件，有可能你关键词没写对，白算了。另外，你用多少核并行，多少内存？

作者
Author: wjc404 时间: 2023-10-18 13:59
本帖最后由 wjc404 于 2023-10-18 14:02 编辑

把gaussian的rwf文件放到高性能的固态硬盘上，然后内存给多一些（从打印的"call FoFJK"来看是没给够的），这样可以让l1111.exe跑得快些。

作者
Author: Zhangchj 时间: 2023-10-18 20:11

wjc404 发表于 2023-10-18 13:59
把gaussian的rwf文件放到高性能的固态硬盘上，然后内存给多一些（从打印的"call FoFJK"来看是没给够的）， ...

感谢指导，刚上手这个软件。还没尝试过分配内存，所以就默认分配了。学到了

作者
Author: Zhangchj 时间: 2023-10-18 20:18

zjxitcc 发表于 2023-10-18 13:12
这个计算量本来就很大，你得对5Z基组有点概念。不过，我还是建议先看看输入文件，有可能你关键词没写对，白 ...

好的感谢指导，刚上手这个软件。用内置基组算了一个很快，还不太了解这个基组，没考虑这个基组差别这么大。输入文件按照sobevera老师的方法写的：

%chk=C:\Users\Zhangchj\Desktop\gaussian_CO\CO\ma-cc-PV5Z.chk
# opt freq ccsd/gen geom=connectivity

ma-cc-PV5Z

0 1
C -3.22012573 1.25157231 0.00000000
O -2.08520335 1.26484625 0.00000000

1 2 2.0
2

@C:\ma-cc-PV5Z.txt

想着只是算两个原子应该不是太慢

核和内存就用默认值了

作者
Author: zjxitcc 时间: 2023-10-18 21:24

Zhangchj 发表于 2023-10-18 20:18
好的感谢指导，刚上手这个软件。用内置基组算了一个很快，还不太了解这个基组，没考虑这个基组差别这么 ...

你离合理的计算还有不少距离。需要多向内行人请教。内存和核数都不指定，那计算效率十分低下。

“用CCSD算偶极矩”这句话并不等于用CCSD做结构优化+频率分析+弛豫密度计算。结构优化和频率分析往往是用DFT算的，得到平衡键长的结构后，用CCSD算单点和弛豫密度。

如果硬要追求结构优化和频率分析也要在CCSD下算（例如导师逼迫），那么对于双原子分子一般是构造一系列结构，每个结构有不同的键长，然后算CCSD单点，得到势能曲线，拟合得到极小点结构和频率，然后只对极小点这一帧结构做弛豫密度计算（从而获得偶极矩）。

作者
Author: sobereva 时间: 2023-10-19 06:09
ma-cc-PV5Z.txt是什么指代不明
如果你是按下文方式加的弥散，完全不适合。把偶极矩算好光靠s、p角动量的弥散根本不够
给def2以ma-方式加弥散函数的Gaussian格式的基组定义文件（含所有def2支持的元素）
http://sobereva.com/509（http://bbs.keinsci.com/thread-14370-1-1.html）

此体系高精度偶极矩、静电势计算用aug-cc-pVQZ就完了，毫无必要考虑cc-pV5Z。另外，倘若你真要极高精度计算偶极矩，应该用CCSD(T)的能量对外电场做有限差分得到，(T)的贡献不可完全忽视

耦合簇结合大基组的计算找个像样的linux服务器，甭拿windows版Gaussian凑合，那基本只是个玩具

作者
Author: Zhangchj 时间: 2023-10-26 22:06

sobereva 发表于 2023-10-19 06:09
ma-cc-PV5Z.txt是什么指代不明
如果你是按下文方式加的弥散，完全不适合。把偶极矩算好光靠s、p角动量的弥 ...

感谢！

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