24小时热门版块排行榜    

查看: 1475  |  回复: 7

crosschannel

金虫 (正式写手)

[求助] 势能面最高点,算出来的振动频率是什么样的?

我知道每个最低点算出红外频率全是正值,因为是坑;鞍点是正负值都有……那么最高点呢?我推测是全是负值?

找到一个点,算出来频率全是正值,有人跟我说:这个就是最高点!
纳尼?最高点也全是正值是吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
crosschannel: 金币+10, ★★★★★最佳答案 2014-07-17 23:34:05
我觉得讨论这个问题是无意义的,因为在全方向上的能量极大点我认为有两种情况:分子键全断,或者所有原子挤一块。这两种情况对于研究分子体系来说有用么?
但是从LZ的表述来看LZ对频率怎么得出的不甚了解,所以我仅以自己的理解来讨论一下LZ的问题吧,欢迎大神们批评指正。
频率是根据此能量上对原子的方向进行二阶求导之后的二阶偏导数矩阵的计算而来。这个矩阵叫Hessien矩阵。一个频率对应一个Hessien矩阵
所以这个问题简单来说是一个线代里的概念问题了,即H矩阵的值是正是负的问题。
在其一阶导数为0或近似于0的情况下:
H矩阵负值,反应在频率上应该为虚频,结合鞍点的意义,可以理解为什么过渡态有一个虚频。
H矩阵正值,还有两种情况:正定和负定。
正定,极小值
负定,极大值
所以以我的理解来说,极大值的频率也是全为正的。
也可能高斯中对于频率的处理进行了简化,这样的情况下全为负也不无可能。
2楼2014-07-17 10:24:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bittersmile

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
crosschannel: 金币+15, ★★★★★最佳答案 2014-07-17 23:43:02
楼上说得不对吧。
1)一个频率不是对应一个Hessien矩阵,而是对应矩阵的一个本征值。
2)因此极大值肯定频率都是负的。
3)pes上的极大值点,并不是完全没研究意义。例如一个二维的pes,二维都是旋转角度的话,这种pes上的极大点,并不涉及任何键断裂,只是键角的旋转而已,完全可以有研究价值。
3楼2014-07-17 23:26:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crosschannel

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2014-07-17 10:24:40
我觉得讨论这个问题是无意义的,因为在全方向上的能量极大点我认为有两种情况:分子键全断,或者所有原子挤一块。这两种情况对于研究分子体系来说有用么?
但是从LZ的表述来看LZ对频率怎么得出的不甚了解,所以我仅 ...

哦,懂了,hessian matrix正值可以是负定,这个我没想到
不过大侠,研究local maximum真不是没意义的,起码现在对我的研究很有影响……
4楼2014-07-17 23:33:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crosschannel

金虫 (正式写手)

引用回帖:
3楼: Originally posted by bittersmile at 2014-07-17 23:26:35
楼上说得不对吧。
1)一个频率不是对应一个Hessien矩阵,而是对应矩阵的一个本征值。
2)因此极大值肯定频率都是负的。
3)pes上的极大值点,并不是完全没研究意义。例如一个二维的pes,二维都是旋转角度的话,这 ...

哎呦我怎么才看见你的回复
5楼2014-07-17 23:43:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crosschannel

金虫 (正式写手)

引用回帖:
3楼: Originally posted by bittersmile at 2014-07-17 23:26:35
楼上说得不对吧。
1)一个频率不是对应一个Hessien矩阵,而是对应矩阵的一个本征值。
2)因此极大值肯定频率都是负的。
3)pes上的极大值点,并不是完全没研究意义。例如一个二维的pes,二维都是旋转角度的话,这 ...

谢谢这位大侠,我就想搞一个不断键的最高点呢,是一大坨原子的旋转
6楼2014-07-17 23:46:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crosschannel

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2014-07-17 10:24:40
我觉得讨论这个问题是无意义的,因为在全方向上的能量极大点我认为有两种情况:分子键全断,或者所有原子挤一块。这两种情况对于研究分子体系来说有用么?
但是从LZ的表述来看LZ对频率怎么得出的不甚了解,所以我仅 ...

这位大侠,你回复下面真有人批评指正你呢,你再来发表一下意见吧
我个菜鸟喜欢看不同观点
7楼2014-07-17 23:48:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

引用回帖:
7楼: Originally posted by crosschannel at 2014-07-17 23:48:06
这位大侠,你回复下面真有人批评指正你呢,你再来发表一下意见吧
我个菜鸟喜欢看不同观点...

关于第一个问题,我看过一些讨论,确实在讲虚频的时候是说频率是hessien的本征值,所以这样一来极大点的频率全是负值。
关于这个,我也询问过我的做量子化学的博士生同学,他说一个hessien矩阵只能给出一个频率,他推荐我去看levien的quantum chemsitry里的推导,但是鄙人不才,还没有找到。。。

第二个问题,在我看来,虽然我可能有点强词夺理,那位大神说的还是在一定的限定条件下,即二维和旋转角扫描。我所说的是就全局来看,如果将分子中的所有原子完全放开了去找这个极大点,那应该是键全处于断裂的状态下的吧?但是你的问题里确实没有限定条件。。。

就这个问题,我觉得可否以H2做这样一个尝试,即用opt=modredundant扫描两个H之间的键长变化,找到极大点,再固定此时的键长opt=modredundant之后算频率?
8楼2014-07-18 08:11:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 crosschannel 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-15 7/350 2026-03-16 09:43 by 闲人终南山
[考研] 304求调剂 +4 素年祭语 2026-03-15 4/200 2026-03-16 09:42 by 闲人终南山
[考研] 321求调剂 +4 大米饭! 2026-03-15 4/200 2026-03-16 08:41 by Linda Hu
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 材料专硕274一志愿陕西师范大学求调剂 +4 薛云鹏 2026-03-13 4/200 2026-03-13 10:40 by 学员8dgXkO
信息提示
请填处理意见