24小时热门版块排行榜

>论坛更新日志 (3697)
>文献求助 (324)
>导师招生 (314)
>虫友互识 (303)
>休闲灌水 (97)
>论文投稿 (81)
>招聘信息布告栏 (72)
>考博 (69)
>硕博家园 (68)
>博后之家 (67)
>公派出国 (52)
>催化 (51)
>考研 (43)
>基金申请 (41)
>论文道贺祈福 (40)
>教师之家 (34)

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chemzhh

至尊木虫 (著名写手)

应助: 0 (幼儿园)
金币: 19616.2
帖子: 1110
在线: 180.6小时
虫号: 94115
注册: 2005-09-26
性别: GG

[交流] 【求助】请问计算稳态结构时出现虚频是什么原因

我用gaussian优化了MePhPhMe(就是两个甲苯对位联起来的结构)的结构后，计算正常结束(normal termination),然后计算其红外频率，发现有一个虚频-51.07，强度为0.0061，这是为什么，频率计算也是正常结束，没有报错。谢谢了。

回复此楼

» 猜你喜欢

1楼 2009-05-17 22:58:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

neweroica

木虫 (著名写手)

应助: 53 (初中生)
金币: 5042.5
散金: 301
红花: 31
帖子: 2480
在线: 216.9小时
虫号: 36775
注册: 2004-02-11
专业: 理论和计算化学

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 5-20 10:08

引用回帖:

Originally posted by chemzhh at 2009-5-19 23:33:
谢谢楼上的详尽解释。
既然通常优化的结构只是一个极小值，或者仅仅只是势能面上的一个驻点，那么根据这个结构计算出的能量有意义吗？还是每个结构必须找到全局最优，然后计算出的能量才行，尽管找到全局最优很难

也许我前面解释得可能有些误导，补充一下：

1. 并非只有全局极小值点才最重要。就好比不同的同分异构体，固然最稳定的那个同分异构体更容易得到，但其他相对能量稍高的异构体同样也可能存在，甚至更有特殊的意义。稳定性还得看不同极小点之间的能垒（决定异构化转换的难易，即动力学稳定性）。如果虽然局部极小点的能量较高，但由它转变成邻近极小点的能垒很高，那么它依然可能稳定存在（即亚稳）。

2. 通常优化的结构仅仅只是势能面上的一个驻点。所以，得进行振动频率的计算来看看有没有虚频率，如果没有，才能真正令人信服地确定是真正的极小值点（但这也取决于计算采用的方法，我就遇到过某个结构在MP2方法下是极小指点，而在B3LYP方法下却存在虚频，相反的情况也遇到过。）

3. 但实际上在很多情况下我们往往并不需要进行频率计算，因为频率计算的代价较大。比如，让你优化苯分子，结果肯定是极小值点，没必要进行频率分析。为什么？因为，我们没有必要说明都依赖于计算化学——因为我们还有化学的基本常识。所以，准备计算之前，以及得到计算结果之后，都要充分运用自己的化学基础知识，不能为计算而计算，不能把程序当成黑盒子。有时候，即使你运用常识，也还是拿不定某个结构究竟是不是极小值点，到底有没有虚频率（特别对于柔性分子，往往会存在一些比较小的虚频率）。我们可以用比较廉价的方法进行优化和频率计算。比如，用MP2进行频率计算所需要的代价要比相同基组下用B3LYP进行频率计算所需要的代价大得很多。这时，你就可以用B3LYP进行频率计算，来确认一下是否是极小值点（当然，这仍然不能保证其在MP2下就一定也是极小值点，但有总比没有好：Algo es algo.）