24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1441  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856专硕求调剂 希望是a区院校 +18 好好休息好不好 2026-04-09 21/1050 2026-04-10 10:35 by l_paradox
[考研] 材料工程日语考生求调剂 +3 0856?调剂 2026-04-10 3/150 2026-04-10 09:45 by yutian743
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 调剂 +3 电气300求调剂不 2026-04-08 6/300 2026-04-08 09:39 by 电气300求调剂不
[考研] 283分求调剂 +14 试试看呗 2026-04-04 14/700 2026-04-08 07:03 by lijunpoly
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 085602调剂 初试总分335 +3 19123253302 2026-04-06 3/150 2026-04-07 18:00 by jp9609
[考研] 材料求调剂 +18 一样YWY 2026-04-05 18/900 2026-04-07 15:49 by dxlg
[考研] 求调剂 +5 小沢 2026-04-03 5/250 2026-04-06 22:45 by 875465
[考研] 295求调剂 +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by 蓝云思雨
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 272求调剂 +4 电气李 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
信息提示
请填处理意见