24小时热门版块排行榜    

查看: 1420  |  回复: 11
当前主题已经存档。

chemzhh

至尊木虫 (著名写手)

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

我用gaussian优化了MePhPhMe(就是两个甲苯对位联起来的结构)的结构后,计算正常结束(normal termination),然后计算其红外频率,发现有一个虚频-51.07,强度为0.0061,这是为什么,频率计算也是正常结束,没有报错。谢谢了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

neweroica

木虫 (著名写手)

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 5-20 10:08
引用回帖:
Originally posted by chemzhh at 2009-5-19 20:54:
能否说频率计算无虚频,那么对应的结构就是全局最优呢?谢谢

当然不是。无虚频只是保证在当前的计算方法下,得到的结构对应势能面上的极小值点(即势能面的各个方向上都是最小值)。但完整的势能面可能有很多极小值点(称为局部极小点),当然还可能有鞍点(在某些个方向上是最小值,但在其他方向上则是极大值。如果是一级鞍点,那就是过渡态)。

在完整的势能面上的所有可能的极小值中,能量最低的那个极小值点就是全局极小值点,对应该分子的热力学最稳定的结构(当然有时还应该考虑零点能)。要找全局极小值点,是个非常非常非常复杂的事情,需要对整个势能面进行全面的扫描。
Simpler, stronger, more tolerant
10楼2009-05-19 22:43:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
lei0736(金币+2,VIP+0):谢谢 5-18 20:02
表明结构还需要进一步优化,不是最稳定结构……
2楼2009-05-18 07:21:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suntao1982

木虫 (著名写手)

小木虫

★ ★
lei0736(金币+2,VIP+0):谢谢 5-18 20:02
还有一个办法,就是找到有虚频的那一列,把原子坐标加上去,再从新计算!
做中国人的化学!!!!
3楼2009-05-18 07:59:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chemzhh

至尊木虫 (著名写手)

还是不太明白

谢谢二楼和三楼的回答,但我还是有疑问。
我用的基组是6-31+g*,为什么优化已经得到四个yes后,结构还要优化。再次优化时用opt=tight是否有助于得到正确结果。
另外三楼说的,如何找到有虚频那一列的坐标。
谢谢了
4楼2009-05-18 09:51:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
lei0736(金币+2,VIP+0):谢谢 5-18 20:03
优化很多时候与初始结构有关。

另外,优化大部分不能获得全局最优,只是局部最优,而这个局部最优也可以满足收敛条件,也可以得到 4 个 yes 的。
5楼2009-05-18 11:54:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chemzhh

至尊木虫 (著名写手)

能否说频率计算无虚频,那么对应的结构就是全局最优呢?谢谢
6楼2009-05-19 20:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★
lei0736(金币+2,VIP+0):谢谢 5-20 10:07
引用回帖:
Originally posted by snoopyzhao at 2009-5-18 11:54:
优化很多时候与初始结构有关。

另外,优化大部分不能获得全局最优,只是局部最优,而这个局部最优也可以满足收敛条件,也可以得到 4 个 yes 的。

非也。 4 个 yes表明是势能面上的驻点,而不一定是极小值点。
Simpler, stronger, more tolerant
7楼2009-05-19 22:35:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★
lei0736(金币+2,VIP+0):谢谢 5-20 10:08
引用回帖:
Originally posted by suntao1982 at 2009-5-18 07:59:
还有一个办法,就是找到有虚频的那一列,把原子坐标加上去,再从新计算!

正解。不过gaussview4已经带有这个功能了,可以自动导出坐标,不再需要我们手工去做了。
Simpler, stronger, more tolerant
8楼2009-05-19 22:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★
lei0736(金币+2,VIP+0):谢谢 5-20 10:08
引用回帖:
Originally posted by chemzhh at 2009-5-17 22:58:
我用gaussian优化了MePhPhMe(就是两个甲苯对位联起来的结构)的结构后,计算正常结束(normal termination),然后计算其红外频率,发现有一个虚频-51.07,强度为0.0061,这是为什么,频率计算也是正常结束,没有报错 ...

怀疑你是不是把两个苯环放在同一个平面内了。因为,按照有机化学常识,联苯是有很大的扭转角,因为氢原子间的Pauli排斥。
Simpler, stronger, more tolerant
9楼2009-05-19 22:37:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chemzhh 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602化学工程求调剂。 +3 平乐乐乐 2026-03-26 3/150 2026-03-26 16:37 by 3Strings
[考研] 0856求调剂 +8 zhn03 2026-03-25 9/450 2026-03-26 13:42 by zzll406
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +6 哈哈157349 2026-03-21 6/300 2026-03-26 10:39 by 醉在风里
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 300求调剂,材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 生物学一志愿985,分数349求调剂 +6 zxts12 2026-03-21 9/450 2026-03-23 18:37 by macy2011
[考研] 315分,诚求调剂,材料与化工085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
信息提示
请填处理意见