24小时热门版块排行榜    

查看: 2933  |  回复: 8

ldzh0531

木虫 (小有名气)

[求助] qst3寻找过渡态时出现不需要的虚频

如题,用qst3寻找过渡态时,在非反应方向上出现了虚频怎么消除。有些反应的反应物和产物的结构除了在反应方向上不一样外,其他地方也有不一样的。理论上,这些不一样的地方都可能出现虚频,可是那些非反应方向的虚频是不想要的,g09算出一个虚频就终止了,根本不对。怎么有办法在输入文件的时候就大致规定一下反应的方向呢?
回复此楼

» 收录本帖的淘帖专辑推荐

催化与分子模拟

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

我做过一个计算,计算出来的虚频就只是单键的旋转,根本不对!也不知道该怎么改
2楼2012-03-02 16:33:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

哪位大神知道一下
3楼2012-03-04 21:41:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenjinmiao

银虫 (初入文坛)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
ldzh0531(金币+3): ★★★很有帮助 2012-03-05 10:52:35
消除虚频的主要方法是改变构型。其次在计算上还可以尝试:nosymm;加大循环次数;提高收敛度;iop(1/8=1)等。通常对于势能面很平,较小的虚频很难消除。当出现两个以上虚频时,可以先看一下振动模式。如果多余的虚频,和想要的振动模式发生在不同原子上,可以先固定过渡态模式所对应的原子,然后做opt到minimal,一般可以把多余的虚频去掉,再放开做全优化。或许会有所收获。
4楼2012-03-05 10:13:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990(金币+2): 2012-03-06 12:21:02
ldzh0531(金币+3): ★★★很有帮助 2012-03-09 12:53:42
半经验的程序组
如 mopac 和 ampac 有一个叫EF的方法 (eigenvalue following)
可指定跟踪某个虚频 (详见ampac manual::EF)
可以先用mopac或ampac优化
得到结果去算gaussian
更痛苦的是梦没醒路已经不能走了
5楼2012-03-06 01:43:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ludeng8710

新虫 (小有名气)

引用回帖:
5楼: Originally posted by mchen10 at 2012-03-06 01:43:42:
半经验的程序组
如 mopac 和 ampac 有一个叫EF的方法 (eigenvalue following)
可指定跟踪某个虚频 (详见ampac manual::EF)
可以先用mopac或ampac优化
得到结果去算gaussian

请问您有 MOPAC的教程 么?“EF的方法 (eigenvalue following)
可指定跟踪某个虚频”具体要怎么做从哪下手呀? 我一直用GAUSSIAN,这个不太会用,谢谢~
6楼2012-03-06 11:18:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)

引用回帖:
6楼: Originally posted by ludeng8710 at 2012-03-05 15:18:43:
请问您有 MOPAC的教程 么?“EF的方法 (eigenvalue following)
可指定跟踪某个虚频”具体要怎么做从哪下手呀? 我一直用GAUSSIAN,这个不太会用,谢谢~

http://openmopac.net/manual/ef.html

用其中的mode=n 来指定某个eigenvalue
更痛苦的是梦没醒路已经不能走了
7楼2012-03-06 11:45:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

引用回帖:
7楼: Originally posted by mchen10 at 2012-03-06 11:45:41:
http://openmopac.net/manual/ef.html

用其中的mode=n 来指定某个eigenvalue

高斯的opt里也有ef关键字,我想问下高斯里的ef是不是跟这个一样。如果是的话该怎么改?
8楼2012-03-09 20:41:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jessie46

木虫 (正式写手)

引用回帖:
4楼: Originally posted by wenjinmiao at 2012-03-05 10:13:47
消除虚频的主要方法是改变构型。其次在计算上还可以尝试:nosymm;加大循环次数;提高收敛度;iop(1/8=1)等。通常对于势能面很平,较小的虚频很难消除。当出现两个以上虚频时,可以先看一下振动模式。如果多余的虚频 ...

请问如何固定过渡态模式所对应的原子,再做opt到minimal呢?
9楼2015-11-25 09:32:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ldzh0531 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 04:47 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见