计算化学公社

标题: 结构优化后到底要不要做频率分析确认？ [打印本页]

作者
Author: hzfish 时间: 2016-1-18 23:21
标题: 结构优化后到底要不要做频率分析确认？
以前一直以为结构优化后必须要在相同水平下做频率分析确认为结构为势能面全局极小点或局部极小点。文献中也经常见到如下类似表达：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方法计算其它性质，是否可行？

谢谢！

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

问题2：完全可行。

作者
Author: hzfish 时间: 2016-1-19 08:38
非常感谢社长指导。

作者
Author: 虎王 时间: 2016-1-19 17:01

sobereva 发表于 2016-1-19 01:09
做振动分析很耗时，如果算得动得话，优化完了做一下振动分析确认无虚频是最保险的。但是当体系较大，比如大 ...

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

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

作者
Author: sobereva 时间: 2016-1-19 17:11

虎王发表于 2016-1-19 17:01
苏忠民教授跟我说过他们有时会优化一系列结构，然后取其中一个算频率，确认没有虚频，然后就说其他的结构 ...

是的。ultrafine和tight一定要结合使用，否则效果都不好
一般对于大分子弱相互作用体系都可以不用振动分析来判断，越是小体系、有对称性的情况，越可能和直觉判断有异，值得振动分析判断一下

作者
Author: 虎王 时间: 2016-1-19 22:06

sobereva 发表于 2016-1-19 17:11
是的。ultrafine和tight一定要结合使用，否则效果都不好
一般对于大分子弱相互作用体系都可以不用振动分 ...

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

作者
Author: sobereva 时间: 2016-1-20 02:32

虎王发表于 2016-1-19 22:06
那在Gaussian 09里能不能先用
opt(loose) int=ultrafine
优化，

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

直接opt=tight int=ultrafine就可以

作者
Author: dingniu2 时间: 2016-1-20 08:51

sobereva 发表于 2016-1-20 02:32
g09默认的是scf=tight，并不是默认的opt=tight，所以要写

直接opt=tight int=ultrafine就可以

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

作者
Author: sobereva 时间: 2016-1-20 09:52

dingniu2 发表于 2016-1-20 08:51
请问g09的int=ultrafine相当于orca的grid多少？

默认的fine格点就是75,302，ultrafine是99,590
对照手册，ultrafine差不多是grid6

欢迎光临计算化学公社 (http://bbs.keinsci.com/)