24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1440  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 087100初试311求调剂 +3 任雅琴 2026-04-09 3/150 2026-04-09 22:42 by lbsjt
[考研] 266求调剂 +25 阳阳哇塞 2026-04-07 25/1250 2026-04-09 20:32 by only周
[考研] 计算机11408,286分求调剂 +9 木子念晞 2026-04-05 9/450 2026-04-09 20:04 by vgtyfty
[考研] 085400 328分 求调剂 +8 喂你一个大橙子 2026-04-09 12/600 2026-04-09 19:20 by zl8213662
[考研] 269求调剂 +7 跪求收留。 2026-04-04 7/350 2026-04-09 19:06 by 探123
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 302分求调剂 +4 凡语祈愿 2026-04-08 5/250 2026-04-08 22:03 by 土木硕士招生
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +20 努力奋斗112 2026-04-07 21/1050 2026-04-08 14:54 by screening
[考研] 301求调剂 +10 细胞相关蛋白 2026-04-03 10/500 2026-04-08 10:36 by tjzhao
[考研] 304求调剂(085602,过四级,一志愿985) +25 化工人999 2026-04-04 26/1300 2026-04-07 22:06 by hemengdong
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 081700学硕,323分,一志愿中国海洋大学求调剂学校 +19 披星河 2026-04-04 19/950 2026-04-07 15:00 by 上岸快快
[考研] 考研调剂 +7 15615482637 2026-04-04 7/350 2026-04-06 22:56 by chenzhimin
[考研] 生物学调剂 可调剂到生物与医药 +3 李政莹 2026-04-06 3/150 2026-04-06 19:02 by macy2011
[考研] 一志愿南航,数一英一学硕317求调剂!! +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 344材料与化工调剂 +9 调剂上岸玘 2026-04-03 9/450 2026-04-04 23:10 by happyddm
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 085601,一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
信息提示
请填处理意见