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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » qst3寻找过渡态时出现不需要的虚频
91/1返回列表
查看: 2261  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭