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凑合,那基本只是个玩具 |
Zhangchj 发表于 2023-10-18 20:18 你离合理的计算还有不少距离。需要多向内行人请教。内存和核数都不指定,那计算效率十分低下。 “用CCSD算偶极矩”这句话并不等于 用CCSD做结构优化+频率分析+弛豫密度计算。结构优化和频率分析往往是用DFT算的,得到平衡键长的结构后,用CCSD算单点和弛豫密度。 如果硬要追求结构优化和频率分析也要在CCSD下算(例如导师逼迫),那么对于双原子分子一般是构造一系列结构,每个结构有不同的键长,然后算CCSD单点,得到势能曲线,拟合得到极小点结构和频率,然后只对极小点这一帧结构做弛豫密度计算(从而获得偶极矩)。 |
参与人数Participants 1 | eV +4 | 收起 理由Reason |
---|---|---|
Zhangchj | + 4 | 谢谢 |
zjxitcc 发表于 2023-10-18 13:12 好的 感谢指导,刚上手这个软件。用内置基组算了一个很快,还不太了解这个基组,没考虑这个基组差别这么大。输入文件按照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 想着只是算两个原子应该不是太慢 核和内存就用默认值了 |
wjc404 发表于 2023-10-18 13:59 感谢指导,刚上手这个软件。还没尝试过分配内存,所以就默认分配了。学到了 |
本帖最后由 wjc404 于 2023-10-18 14:02 编辑 把gaussian的rwf文件放到高性能的固态硬盘上,然后内存给多一些(从打印的"call FoFJK"来看是没给够的),这样可以让l1111.exe跑得快些。 |
本帖最后由 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.