24小时热门版块排行榜    

查看: 1402  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 8/400 2026-03-16 09:45 by xujiaoszu
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
信息提示
请填处理意见