应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » qst3寻找过渡态时出现不需要的虚频
91/1返回列表
查看: 2742  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ldzh0531

木虫 (小有名气)

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

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

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

催化与分子模拟

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-03-01 19:56:16
已阅   回复此楼   关注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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？ +5 13108017953 2026-09-21 5/250 2026-09-27 16:46 by Leogzhya
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +3 zju2000 2026-09-26 3/150 2026-09-27 16:26 by lala1008
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:26 by 4IW0sJvtEqX8
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:21 by 4IW0sJvtEqX8
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:15 by 4IW0sJvtEqX8
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ViqFlLSxHDGt 2026-09-26 3/150 2026-09-27 10:03 by 4IW0sJvtEqX8
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 GxESXJmptOtQ 2026-09-25 3/150 2026-09-27 06:32 by XgCC8uTcwILl
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
信息提示
关闭
请填处理意见
关闭