24小时热门版块排行榜    

查看: 1419  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +15 tcxiaoxx 2026-03-20 17/850 2026-03-26 13:27 by 小小麦片
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +5 丹青奶盖 2026-03-26 5/250 2026-03-26 12:55 by 0906ljy
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 材料求调剂 +4 .m.. 2026-03-25 4/200 2026-03-25 21:30 by peike
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 求调剂一志愿武汉理工大学材料工程(085601) +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
信息提示
请填处理意见