24小时热门版块排行榜    

查看: 2462  |  回复: 10

vigaryang

木虫 (正式写手)

[交流] 【讨论】消除虚频时是采用scf=tight还是opt=tight?

优化两个分子的相互作用,出现一个虚频(-7.2)。根据本论坛的FAQ,似乎应该是以上述带虚频的构型为初始构型、命令行中加入SCF=tight关键字后继续优化。不过我在网上看到有其他人说应该加入OPT=tight命令继续优化。不知这两种情况孰是孰非,抑或两者皆可?如果采用opt=tight的话,是否应如gaussian 03手册中所说一定要再加上Int=UltraFine关键词?谢谢。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

userhung

禁虫 (文学泰斗)

木虫博士

10楼2010-03-27 08:27:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

c111999

铜虫 (小有名气)

★ ★ ★ ★
vigaryang(金币+1):谢谢参与
yjcmwgk(金币+3): 2010-03-12 13:24
引用回帖:
Originally posted by vigaryang at 2010-03-10 22:13:09:
优化两个分子的相互作用,出现一个虚频(-7.2)。根据本论坛的FAQ,似乎应该是以上述带虚频的构型为初始构型、命令行中加入SCF=tight关键字后继续优化。不过我在网上看到有其他人说应该加入OPT=tight命令继续优化 ...

我觉得应该opt=tgigjht
你也可以在加上scf=tighth
这个虚品多半可以忽略
加ultrafine如果你用了dft
2楼2010-03-10 22:22:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

引用回帖:
Originally posted by c111999 at 2010-03-10 22:22:10:

我觉得应该opt=tgigjht
你也可以在加上scf=tighth
这个虚品多半可以忽略
加ultrafine如果你用了dft

如果忽略的话,对相互作用能会有多大的影响呢?虚频在什么样的情况下可以考虑忽略呢?谢谢。
3楼2010-03-10 22:41:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c111999

铜虫 (小有名气)

★ ★
yjcmwgk(金币+2): 2010-03-12 13:24
小虚频对于相互作用的影响,我不知道多大,估计不大。

过渡态需要的虚频
1应该是在要断键/成键的方向上的震动
2应该比较大~几百波数
3几十~100,争取减少,如果无法消减——如果不在1内,则可以忽略。
4平动转动也会产生虚频,很小就是了,看log文件中的“Low frequencies”

[ Last edited by c111999 on 2010-3-10 at 22:57 ]
4楼2010-03-10 22:48:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

我的log文件如下:
Low frequencies ---  -11.9736  -10.4981   -6.7810   -0.0081    0.0114    0.0118
Low frequencies ---    8.6975   24.4838   34.1247
虽然这里有好几负数,但随后的频率输出结果里只有一个:
                    1                      2                      3
                     A                      A                      A
Frequencies --    -7.1160                24.0276                34.0197
Red. masses --     4.1740                 4.8112                 4.9130
Frc consts  --     0.0001                 0.0016                 0.0034
IR Inten    --     1.0249                 1.4418                 2.4991
这说明什么问题呢?
5楼2010-03-10 23:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c111999

铜虫 (小有名气)

vigaryang(金币+1): 2010-03-10 23:42
平动+转动Low frequencies ---  -11.9736  -10.4981   -6.7810   -0.0081    0.0114    0.0118
最低震动Low frequencies ---    8.6975   24.4838   34.1247

-7 ~ 平动转动的频率,因此可以忽略

奇怪的是Frequencies --    -7.1160   =!=  8.6975   不知是如何算的?

[ Last edited by c111999 on 2010-3-10 at 23:25 ]
6楼2010-03-10 23:22:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

引用回帖:
Originally posted by c111999 at 2010-03-10 23:22:08:
平动+转动Low frequencies ---  -11.9736  -10.4981   -6.7810   -0.0081    0.0114    0.0118
最低震动Low frequencies ---    8.6975   24.4838   34.1247

-7 ~ 平动转动的频率,因此可以忽略

奇怪的是F ...

后面的两个最终频率也并不是和最低震动频率完全吻合。看来最后给出的编了号的频率应该是“平动+转动”和“震动”两者的一个函数?
7楼2010-03-10 23:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

467266420

铁杆木虫 (著名写手)


vigaryang(金币+1):谢谢参与
引用回帖:
Originally posted by vigaryang at 2010-03-10 22:13:09:
优化两个分子的相互作用,出现一个虚频(-7.2)。根据本论坛的FAQ,似乎应该是以上述带虚频的构型为初始构型、命令行中加入SCF=tight关键字后继续优化。不过我在网上看到有其他人说应该加入OPT=tight命令继续优化 ...

无论什么方法能消掉虚频才是王道,提供一个消虚频的小软件http://muchong.com/bbs/viewthread.php?tid=1874636  类似于小卒大师在本版GaussianFAQ提到的坐标中合法。

[ Last edited by yjcmwgk on 2010-3-14 at 12:14 ]
8楼2010-03-12 19:21:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

引用回帖:
Originally posted by 467266420 at 2010-03-12 19:21:02:

无论什么方法能消掉虚频才是王道,提供一个消虚频的小软件http://muchong.com/bbs/viewthread.php?tid=1874636&pid=18892365&page=1#pid18892365  类似于小卒大师在本版GaussianFAQ提到的坐标 ...

非常感谢!
9楼2010-03-12 22:48:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vigaryang 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 CyTWnZiM0cS6 2026-08-02 5/250 2026-08-02 20:43 by 3VGOm6paSgqv
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:32 by 3VGOm6paSgqv
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:31 by 3VGOm6paSgqv
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 AS6fX1tr1on4 2026-08-02 5/250 2026-08-02 20:31 by 3VGOm6paSgqv
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:24 by 3VGOm6paSgqv
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 20:14 by 3VGOm6paSgqv
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:11 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 AS6fX1tr1on4 2026-08-02 14/700 2026-08-02 20:11 by 3VGOm6paSgqv
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 10/500 2026-08-02 20:03 by 3VGOm6paSgqv
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cBGhfJFuNMca 2026-08-02 17/850 2026-08-02 19:54 by 3VGOm6paSgqv
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 19:54 by 3VGOm6paSgqv
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 cBGhfJFuNMca 2026-08-02 21/1050 2026-08-02 19:52 by 3VGOm6paSgqv
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 18:34 by ca3DVAx4OB19
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[基金申请] 准备明年的基金了 +7 Tide man 2026-07-27 7/350 2026-08-02 16:43 by gltch
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
信息提示
请填处理意见