计算化学公社

标题: 关于CDA分析原子顺序不一致的问题 [打印本页]

作者
Author:
小李子赛    时间: 2020-4-13 12:37
标题: 关于CDA分析原子顺序不一致的问题
本帖最后由 小李子赛 于 2020-4-23 16:00 编辑

老师好,分析CDA,我在单点计算复合物时,发现其原子顺序发生了变化,两个片段的原子顺序未发生变化,不知道是哪里出了问题?我加了弥散函数复合物:

输入文件关键词及坐标
#p MP2/genecp nosymm pop=full iop=(3/33=1,3/32=2) 5d 7f











作者
Author:
snljty    时间: 2020-4-13 12:55
本帖最后由 snljty 于 2020-4-13 13:23 编辑

顺便一提,用MP2通常已经不是什么好主意了,用MP2做CDA分析还不用Density=Current(还得改Population=Full,限制型方法改Population=NO,非限制型方法改Population=NOAB),用的是HF级别的波函数做的分析,是很浪费计算力而且低效率的方法。还有你的弥散函数...刚看到GenECP没注意
作者
Author:
sobereva    时间: 2020-4-13 13:01
加了nosymm的情况下,不太可能出现这种情况。上传你的输入输出文件,我看看。

如果没有特殊理由,不建议用MP2,用DFT可以直接用fch,还省得写pop=full iop=(3/33=1,3/32=2)。

如我博文(http://sobereva.com/166)、Multiwfn手册、GCDA论文(http://www.hanspub.org/journal/PaperInformation.aspx?paperID=16370)所述,做CDA强烈不建议带弥散函数,可能导致结果根本没意义或者有极大的误导性。

作者
Author:
小李子赛    时间: 2020-4-13 17:28
sobereva 发表于 2020-4-13 13:01
加了nosymm的情况下,不太可能出现这种情况。上传你的输入输出文件,我看看。

如果没有特殊理由,不建议 ...

老师,由于文件太大上传不了,我把输入输出文件放网盘里了,麻烦老师帮我看一下哪里出了问题,链接:https://pan.baidu.com/s/1AKVuUmMQJB8FdlR3CNsYSQ
提取码:gmxm

作者
Author:
小李子赛    时间: 2020-4-13 17:33
snljty 发表于 2020-4-13 12:55
顺便一提,用MP2通常已经不是什么好主意了,用MP2做CDA分析还不用Density=Current(还得改Population=Full, ...

我是看前几届的师兄师姐用的这种计算方法,类比算的,老师,如果在算CDA时与之前的优化等其它计算不在同一基组以及计算水平下,对计算结果有影响吗?
作者
Author:
snljty    时间: 2020-4-14 19:32
小李子赛 发表于 2020-4-13 17:33
我是看前几届的师兄师姐用的这种计算方法,类比算的,老师,如果在算CDA时与之前的优化等其它计算不在同 ...

问题是你写MP2不写Density,结果和用HF没区别啊...
作者
Author:
小李子赛    时间: 2020-4-14 20:06
snljty 发表于 2020-4-14 19:32
问题是你写MP2不写Density,结果和用HF没区别啊...

哦哦,我明白了。非常感谢

作者
Author:
sobereva    时间: 2020-4-15 00:14
小李子赛 发表于 2020-4-13 17:33
我是看前几届的师兄师姐用的这种计算方法,类比算的,老师,如果在算CDA时与之前的优化等其它计算不在同 ...

显然没影响

不要盲目效仿师兄师姐。实验室祖传的以讹传讹的错误东西就是这么传下来的。
作为CDA代码的开发者(以及GCDA方法的提出者),我认为我在这方面比你师兄师姐更有发言权。

不仅用MP2对于CDA没有任何好处,单纯从理论计算的常识来说,用MP2通常也是非常糟糕、很容易被审稿人骂的做法,仔细看下文中的相关部分
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html


作者
Author:
sobereva    时间: 2020-4-15 00:15
小李子赛 发表于 2020-4-13 17:28
老师,由于文件太大上传不了,我把输入输出文件放网盘里了,麻烦老师帮我看一下哪里出了问题,链接:http ...

网盘失效了
甭用pop=full,文件能小得多得多


作者
Author:
小李子赛    时间: 2020-4-15 10:50
sobereva 发表于 2020-4-15 00:15
网盘失效了
甭用pop=full,文件能小得多得多

链接:https://pan.baidu.com/s/1--T-uPkhUp4b6Z3ySX_UPw
提取码:euwh
不好意思老师,我没有注意看有效时间,这个是输入输出文件。
作者
Author:
小李子赛    时间: 2020-4-15 11:03
sobereva 发表于 2020-4-15 00:15
网盘失效了
甭用pop=full,文件能小得多得多

老师,同样的关键词,我使用B3LYP计算方法,得到的坐标还是这个样子。
作者
Author:
sobereva    时间: 2020-4-16 13:31
小李子赛 发表于 2020-4-15 11:03
老师,同样的关键词,我使用B3LYP计算方法,得到的坐标还是这个样子。

你的输出文件里出现了两次nosymm,完全是不合理的




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3