24小时热门版块排行榜    

查看: 2112  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 311求调剂 +13 冬十三 2026-03-15 14/700 2026-03-21 22:10 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 0703化学调剂 +11 妮妮ninicgb 2026-03-15 15/750 2026-03-21 19:15 by ColorlessPI
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 材料专硕英一数二306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 5/250 2026-03-20 22:11 by 云游重阳
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
信息提示
请填处理意见