计算化学公社

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

[综合交流] 结构优化后到底要不要做频率分析确认?

[复制链接 Copy URL]

219

帖子

0

威望

3614

eV
积分
3833

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
          以前一直以为结构优化后必须要在相同水平下做频率分析确认为结构为势能面全局极小点或局部极小点。文献中也经常见到如下类似表达:The structures were optimized using the aug-cc-pVTZ  Dunning’s correlated consistent basis set29 and the B3LYP functional.30 Frequency calculations were carried out to check that all the geometries correspond to energy minima.

最近看了一些弱相互作用的文件,出现两种情况:
一、结构优化并做频率分析确认
如文献:
DOI: 10.1021/jp512020q J. Phys. Chem. C 2015, 119, 7095−7108
dx.doi.org/10.1021/jp506277n | J. Phys. Chem. C 2014, 118, 26003−26016

二、只做结构优化,没有频率分析确认
如文献:
A.A. Peyghan et al. / Superlattices and Microstructures 59 (2013) 115–122
dx.doi.org/10.1021/jp503126q | J. Phys. Chem. C 2014, 118, 16165−16174
S. M. Seyed-Talebi and M. Neek-Amal J. Appl. Phys. 116, 153507 (2014)
Wang, Zhang, and Wang J. Chem. Phys. 140, 094302 (2014)
DOI 10.1007/s00894-015-2768-0
dx.doi.org/10.1021/jp411979x | J. Phys. Chem. C 2014, 118, 2650−2659
Journal of Computational Chemistry 2015, 36, 1763–1771

我有两个问题需要确认一下:第一个问题对于弱相互作用的复合物结构优化后到底要不要做频率分析确认?第二个问题,用半经验方法计算复合物结构,然后用DFT方法计算其它性质,是否可行?

谢谢!

评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
翟yingxin123 + 3 不明觉厉
一条君 + 1 好物!

查看全部评分 View all ratings

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2016-1-19 01:09:49 | 只看该作者 Only view this author
做振动分析很耗时,如果算得动得话,优化完了做一下振动分析确认无虚频是最保险的。但是当体系较大,比如大的弱相互作用体系,光是优化就挺费时间了,再做振动分析耗时太高了,而且很多情况会由于内存、硬盘占用过大而根本没法做。对于有经验的人,对于熟悉的体系从直觉上一般就可以直接判断出是否优化到了极小点,所以振动分析不是必须要做的。

问题2:完全可行。
北京科音自然科学研究中心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

219

帖子

0

威望

3614

eV
积分
3833

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2016-1-19 08:38:25 | 只看该作者 Only view this author
非常感谢社长指导。

294

帖子

2

威望

3018

eV
积分
3352

Level 5 (御坂)

4#
发表于 Post on 2016-1-19 17:01:20 | 只看该作者 Only view this author
sobereva 发表于 2016-1-19 01:09
做振动分析很耗时,如果算得动得话,优化完了做一下振动分析确认无虚频是最保险的。但是当体系较大,比如大 ...

苏忠民教授跟我说过他们有时会优化一系列结构,然后取其中一个算频率,确认没有虚频,然后就说其他的结构也很可能是极小值点的结构。

另外,优化时加int=ultrafine应该可以较大程度上减少出现虚频的概率吧?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

5#
发表于 Post on 2016-1-19 17:11:20 | 只看该作者 Only view this author
虎王 发表于 2016-1-19 17:01
苏忠民教授跟我说过他们有时会优化一系列结构,然后取其中一个算频率,确认没有虚频,然后就说其他的结构 ...

是的。ultrafine和tight一定要结合使用,否则效果都不好
一般对于大分子弱相互作用体系都可以不用振动分析来判断,越是小体系、有对称性的情况,越可能和直觉判断有异,值得振动分析判断一下
北京科音自然科学研究中心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

294

帖子

2

威望

3018

eV
积分
3352

Level 5 (御坂)

6#
发表于 Post on 2016-1-19 22:06:57 | 只看该作者 Only view this author
sobereva 发表于 2016-1-19 17:11
是的。ultrafine和tight一定要结合使用,否则效果都不好
一般对于大分子弱相互作用体系都可以不用振动分 ...

那在Gaussian 09里能不能先用
opt(loose) int=ultrafine
优化,
再用
opt  int=ultrafine
优化?
还是说这样不如直接用
opt  int=ultrafine
优化?
至于opt=tight这个应该是Gaussian 09默认的,所以不用写。

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

7#
发表于 Post on 2016-1-20 02:32:13 | 只看该作者 Only view this author
虎王 发表于 2016-1-19 22:06
那在Gaussian 09里能不能先用
opt(loose) int=ultrafine
优化,

g09默认的是scf=tight,并不是默认的opt=tight,所以要写

直接opt=tight int=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

545

帖子

0

威望

5972

eV
积分
6517

Level 6 (一方通行)

8#
发表于 Post on 2016-1-20 08:51:51 | 只看该作者 Only view this author
sobereva 发表于 2016-1-20 02:32
g09默认的是scf=tight,并不是默认的opt=tight,所以要写

直接opt=tight int=ultrafine就可以

请问g09的int=ultrafine相当于orca的grid多少?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

9#
发表于 Post on 2016-1-20 09:52:33 | 只看该作者 Only view this author
dingniu2 发表于 2016-1-20 08:51
请问g09的int=ultrafine相当于orca的grid多少?

默认的fine格点就是75,302,ultrafine是99,590
对照手册,ultrafine差不多是grid6
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 17:47 , Processed in 0.271567 second(s), 21 queries , Gzip On.

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