计算化学公社

标题: 密度泛函计算出虚频,MP2没有虚频 [打印本页]

作者
Author:
chemjay    时间: 2016-3-30 16:43
标题: 密度泛函计算出虚频,MP2没有虚频
遇到一个体系,密度泛函用高HF的泛函如m062x计算得到线性结构,但是有虚频,B3lyp计算得到弯曲结构,也有虚频,
ccsd(t)和mp2得到线性结构,没有虚频。请问是不是在这个情况中密度泛函就不适用了??或者有没有什么解决的办法?

作者
Author:
冰释之川    时间: 2016-3-30 19:50
B3LYP下优化不行,不代表其它泛函优化不行呀
作者
Author:
liyuanhe211    时间: 2016-3-31 08:54
什么体系?
CCSD(t)优化太浪费,CCSD足够了
作者
Author:
chemjay    时间: 2016-3-31 09:30
冰释之川 发表于 2016-3-30 19:50
B3LYP下优化不行,不代表其它泛函优化不行呀

M06等高hf的泛函优化得到线形,但是有虚频。原帖中的图一个是m06一个是ccsd(t)

作者
Author:
chemjay    时间: 2016-3-31 09:31
liyuanhe211 发表于 2016-3-31 08:54
什么体系?
CCSD(t)优化太浪费,CCSD足够了

问题关键不在于计算浪费资源问题,在于DFT不能得到合理的结果。这个问题应该怎么解决
作者
Author:
chemjay    时间: 2016-3-31 09:34
补充说明一下,帖子中虚频的是M062x计算的结果,没有虚频的是CCSD(T)的结果,两个计算的结构相似,但是mm062x的频率计算貌似
有问题, B3lyp等低hf的泛函计算的结构和CCSD(T)都相差很远,一个是线形一个是弯曲的
作者
Author:
冰释之川    时间: 2016-3-31 10:58
chemjay 发表于 2016-3-31 09:34
补充说明一下,帖子中虚频的是M062x计算的结果,没有虚频的是CCSD(T)的结果,两个计算的结构相似,但是mm06 ...

http://www.sobereva.com/278    这个试过了吗
作者
Author:
liyuanhe211    时间: 2016-3-31 11:36
chemjay 发表于 2016-3-31 09:31
问题关键不在于计算浪费资源问题,在于DFT不能得到合理的结果。这个问题应该怎么解决

(, 下载次数 Times of downloads: 28)



作者
Author:
beefly    时间: 2016-3-31 14:16
有些体系用DFT算不好。例如线形H-He-CCH,CCSD(T)和MP2有虚频,各种DFT都没有虚频
作者
Author:
chemjay    时间: 2016-3-31 16:21
beefly 发表于 2016-3-31 14:16
有些体系用DFT算不好。例如线形H-He-CCH,CCSD(T)和MP2有虚频,各种DFT都没有虚频

对这个体系奇怪的是,不杂化HF交换相关的时候得到的结构和CCSDt差别很大, 含高杂化高HF交换相关的泛函
计算的结构和CCSD(t)接近,但频率有虚频。像这种情况是不是可不可以通过调整所用的密度泛函来得到正确的结果
作者
Author:
sobereva    时间: 2016-4-3 12:06
你先确认是不是因为数值问题所致的,观看虚频振动模式考察对应什么模式。M062X计算时用int=ultrafine和tight收敛限再试。数值问题很容易导致小虚频出现。
另外可以再试双杂化泛函。





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