计算化学公社

 找回密码 Forget password
 注册 Register
Views: 16717|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 关于opt=calcall 的频率分析出现虚频

[复制链接 Copy URL]

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

本帖最后由 冰释之川 于 2017-6-12 20:23 编辑

用opt=calcall 优化了一个花菁分子,但频率分析发现有2个甲基旋转的虚频。查阅了《Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法http://sobereva.com/278 。其中PS2有一句写道“几何优化的时候用了GDIIS(默认的GEDIIS也有一定GDIIS的成份),由于这种方法预测下一步位移的时候还会参考之前步的信息,因此会和freq时基于牛顿法判断出的位移有所不同。”那么结合本例,出现两个虚频是否还能接受呢? 另外,有没有可能是对称性过高导致的呢?

T11.com (3.4 KB, 下载次数 Times of downloads: 18)

T11.log (559.28 KB, 下载次数 Times of downloads: 48)



Stand on the shoulders of giants

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2017-6-12 19:54:10 | 只看该作者 Only view this author
就是对称性太高。旋转一下甲基,破坏Cs对称性在优化即可
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

3#
发表于 Post on 2017-6-28 17:15:12 | 只看该作者 Only view this author
本帖最后由 mutron 于 2017-6-28 17:18 编辑

(1)试了一下选用别的基组或者使用更好的积分格点会得到C2v对称的ground state,对于有对称性又不太小的分子是否需要考虑使用更好的积分格点以排除数值问题导致的小虚频?(参见此贴:密度泛函计算出虚频,MP2没有虚频
T11-C2v
631G(d)
CBSB7
GTBas3
cc-pvtz
wb97xD-fine
IF


wb97xD-ultrafine
G
IF
IF

wb97xD-superfine
G
G
G
G
(IF:虚频,G:基态)

丙酮分子也碰到过这个问题,参考一下:
Acetone-C2v
631G(d)
CBSB7
GTBas3
cc-pvtz
wb97xD-fine
IF
IF
IF
IF
wb97xD-ultrafine
IF
IF
G
G
wb97xD-superfine
G
IF
G
G
CCSD


G
B2GPPLYP-ultrafine


G
DSDPBEP86-D2-ultrafine


G
(丙酮基态应该是C2v对称,参考NIST-acetonefar infrared torsional spectra of acetoneacetone in excited states



(2)如果一个分子几个不同对称性都没有虚频,计算自由能的时候该怎么选?选对称性最高的那个?
T11-wb97xD-ultrafine
C1
C2
Cs
C2v
Gibbs correction
0.391535
0.392657
0.391542
0.391365
kcal/mol
0.11
0.81
0.11
0.00
不同对称性算出来的Gibbs校正还是有点差别(几个不同对称的几何结构几乎看不出来变化,SCF能量没差别)

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2017-6-28 18:47:40 | 只看该作者 Only view this author
mutron 发表于 2017-6-28 17:15
(1)试了一下选用别的基组或者使用更好的积分格点会得到C2v对称的ground state,对于有对称性又不太小的分 ...

wB97XD/6-31G* int=fine确实有虚频,ultrafine时候倒是没发现虚频。
当前问题应当是仅限于wB97XD或个别泛函,测试了B3LYP、PBE0、M062X、B97D,在fine级别皆无此问题
由于wB97XD比起M062X也没什么优点,索性别用它

M062X在某些情况我也发现会由于格点精度问题导致破坏实际对称性的虚频。比如H2O--HBr,在fine格点下居然不是平面的。

一些泛函格点敏感性问题主要是在涉及到弱相互作用的时候容易凸显,上面的二聚体是这样,丙酮的甲基间的微弱耦合也是这样。

对于wB97XD、明尼苏达系列,用ultrafine比较保险。G16倒也是已经默认ultrafine了。

热校正量还是应当在没有虚频、能得到真实点群的格点设定下计算
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

5#
发表于 Post on 2017-6-28 19:42:31 | 只看该作者 Only view this author
本帖最后由 mutron 于 2017-6-28 20:19 编辑
sobereva 发表于 2017-6-28 18:47
wB97XD/6-31G* int=fine确实有虚频,ultrafine时候倒是没发现虚频。
当前问题应当是仅限于wB97XD或个别 ...

感谢sober老师的解答

想再请教一下,如果暂时没找到分子真实点群的资料的话,有没有什么约定俗成的规定使用哪个点群(都没有虚频,碰到过的几个有C2点群的分子和完全没有对称C1计算出来的热校正量都差了0.4-0.8kcal/mol)

PS:其实丙酮的B3LYP也会在某些基组+fine格点下会有这个问题

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2017-6-29 03:19:38 | 只看该作者 Only view this author
mutron 发表于 2017-6-28 19:42
感谢sober老师的解答

想再请教一下,如果暂时没找到分子真实点群的资料的话,有没有什么约定俗成的规 ...

DFT的话用superfine绝对保险。要么用后HF方法,没有积分格点的影响。
如果得到了不同点群下的极小点,比较能量高低判断该用哪个。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

7#
发表于 Post on 2017-6-29 16:18:08 | 只看该作者 Only view this author
sobereva 发表于 2017-6-29 03:19
DFT的话用superfine绝对保险。要么用后HF方法,没有积分格点的影响。
如果得到了不同点群下的极小点,比 ...

“比较能量高低判断该用哪个”   是不是然后选取Gibbs能量最低的那个点群来计算(并且其他点群就不再考虑)?

因为题主给出的T11在wB97XD/6-31G(d) int=ultrafine条件下C1、C2、Cs、C2v点群都是极小点,而且空间结构只有细微坐标差别,电子能量没区别,但是算出来热校正量有不小区别。这种仅仅因为点群识别不一样导致结果有差别让我有点困惑。而且按道理四个点群都是极小点的话,应该都需要统计进玻尔兹曼分布(把空间结构完全一样只是点群不一样的分子看做异构体)?这看上去又好像不太对。。。

碰到过C2点群算出来的Gibbs自由能比完全不考虑对称的C1点群高,这个是不是只能选取没有对称的C1点群进行后续计算(以前算过的两个分子都是如此,包括题主的T11 C2点群Gibbs校正也是比C1的大,不知道是不是普遍都是这样)?以前碰到分子有对称的我默认都用gview识别一下全部用对称来计算,如果正确做法是需要判断Gibbs自由能高低来选取点群恐怕我的一些结果得重算了   T.T


6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

8#
发表于 Post on 2017-6-29 21:49:25 | 只看该作者 Only view this author
mutron 发表于 2017-6-29 16:18
“比较能量高低判断该用哪个”   是不是然后选取Gibbs能量最低的那个点群来计算(并且其他点群就不再考虑 ...

取Gibbs能量最低的那个点群

"点群识别不一样导致结果有差别"这是因为计算配分函数时涉及转动对称数,转动对称数是和点群有关系的。要让程序判断对实际点群,否则转动对热力学数据的贡献是不准的。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

9#
发表于 Post on 2017-6-30 17:07:03 | 只看该作者 Only view this author
sobereva 发表于 2017-6-29 21:49
取Gibbs能量最低的那个点群

"点群识别不一样导致结果有差别"这是因为计算配分函数时涉及转动对称数, ...

如果是取Gibbs能量最低的那个点群的话,可能会在(某些?)C2、C2h点群选择的时候碰到问题

比如:1,2-Dichloroethane

1,2-Dichloroethane(anti)
C1
C2h
total electronic energy
-997.8598602
-997.8598603
Gibbs correction
0.030038
0.030692
Gibbs free energy
-997.8298222
-997.8291683
(MP2/cc-pVTZ)

可以看到电子能量没区别,但是最后Gibbs自由能C1反而更低(~0.4kcal/mol),按照给出条件需要选择C1点群。
但是1,2-Dichloroethane(anti构象)不少都归到C2h点群,gauche构象归到C2点群(参考:NISTInfrared and Raman Spectra of 1,2-DichloroethaneStructure of Molecules and Internal Rotation, Page 32

所以是不是优先考虑最高点群会更合适(物质结构还没有对应实验资料)?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

10#
发表于 Post on 2017-7-1 01:42:05 | 只看该作者 Only view this author
mutron 发表于 2017-6-30 17:07
如果是取Gibbs能量最低的那个点群的话,可能会在(某些?)C2、C2h点群选择的时候碰到问题
...


如果优化出来的结构,只差一丁点(肉眼都几乎难以分辨)就能归到更高对称上,这往往是因为数值精度等原因造成的虚假的对称性下降(积分格点精度不够、用了弥散函数等容易导致),这时候应该gview里对称化成高对称后优化。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

88

帖子

0

威望

1955

eV
积分
2043

Level 5 (御坂)

11#
发表于 Post on 2017-7-1 12:38:11 | 只看该作者 Only view this author
sobereva 发表于 2017-7-1 01:42
如果优化出来的结构,只差一丁点(肉眼都几乎难以分辨)就能归到更高对称上,这往往是因为数值精度等原 ...

嗯,谢谢老师耐心的解答  :D

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 05:29 , Processed in 0.209434 second(s), 23 queries , Gzip On.

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