计算化学公社

 找回密码 Forget password
 注册 Register

想算carbonmonoxide偶极矩,静电势。CCSD-cc-PV5Z基组已经一天了

查看数: 626 | 评论数: 7 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-10-18 08:57

正文摘要:

本帖最后由 Zhangchj 于 2023-10-18 09:00 编辑 看到其他人的文章中采用这个基组,所以选了相同基组。基组来源于BSE网站,内置的基组大概也就十几秒就算完了,

回复 Reply

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

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

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

耦合簇结合大基组的计算找个像样的linux服务器,甭拿windows版Gaussian凑合,那基本只是个玩具
zjxitcc 发表于 Post on 2023-10-18 21:24:38
Zhangchj 发表于 2023-10-18 20:18
好的 感谢指导,刚上手这个软件。用内置基组算了一个很快,还不太了解这个基组,没考虑这个基组差别这么 ...

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

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

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

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
Zhangchj + 4 谢谢

查看全部评分 View all ratings

Zhangchj 发表于 Post on 2023-10-18 20:18:34
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

想着只是算两个原子应该不是太慢 核和内存就用默认值了
Zhangchj 发表于 Post on 2023-10-18 20:11:31
wjc404 发表于 2023-10-18 13:59
把gaussian的rwf文件放到高性能的固态硬盘上,然后内存给多一些(从打印的"call FoFJK"来看是没给够的), ...

感谢指导,刚上手这个软件。还没尝试过分配内存,所以就默认分配了。学到了
wjc404 发表于 Post on 2023-10-18 13:59:22
本帖最后由 wjc404 于 2023-10-18 14:02 编辑

把gaussian的rwf文件放到高性能的固态硬盘上,然后内存给多一些(从打印的"call FoFJK"来看是没给够的),这样可以让l1111.exe跑得快些。
zjxitcc 发表于 Post on 2023-10-18 13:12:59
本帖最后由 zjxitcc 于 2023-10-18 13:25 编辑

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

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 07:52 , Processed in 0.180250 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list