24小时热门版块排行榜    

查看: 2643  |  回复: 10

cczjkd

木虫 (正式写手)

[求助] 高斯计算负频如何修改

我用高斯计算,计算方式是opt+freq在计算结尾是 Normal termination of Gaussian 03 结束的。在查看频率时出现负频。
Frequencies --   -18.6466                37.8502                83.8667
Frequencies --    85.1228                86.5643                89.2846
Frequencies --   101.3354               119.4033               120.7554
Frequencies --   125.8665               126.1519               132.6895
Frequencies --   144.7401               159.4350               167.4659
Frequencies --   176.2744               182.7707               192.3313
Frequencies --   194.9447               199.1771               241.0531
Frequencies --   288.7987               295.7388               371.3347
Frequencies --   432.3451               435.5862              1841.7032
请问如何修改呢?
回复此楼

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

经验

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

绍敏郡主

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gkf高: 金币+3 2013-11-28 22:24:51
cczjkd: 金币+10, 有帮助 2013-12-01 16:06:16
引用回帖:
4楼: Originally posted by cczjkd at 2013-11-28 20:24:41
还有输入的坐标是:
W
W                    1    B1
W                    2    B2       1    A1
W                    1    B3       2    A2       3    D1       0
W                    3    B4    ...

在输出文件里有每个振动所对应的坐标。
用gaussview打开输出文件,进入“results”,“vibrations”,就能看到各个振动的列表,选中虚频,勾选最下面一行manual displacement,把那个数改成0.2,然后右边“save structure”就行了。如图。
高斯计算负频如何修改
屏幕快照 2013-11-28 下午9.33.54.png

5楼2013-11-28 20:38:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

绍敏郡主

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gkf高: 金币+2 2013-11-28 22:24:38
将虚频对应的坐标乘上一个系数比如0.2,然后与优化后的分子坐标相加,以此作为初始坐标,再进行优化
2楼2013-11-28 20:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cczjkd

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 绍敏郡主 at 2013-11-28 20:13:59
将虚频对应的坐标乘上一个系数比如0.2,然后与优化后的分子坐标相加,以此作为初始坐标,再进行优化

你好,谢谢你的回答。但是如何知道负频和哪个坐标对应呢?
3楼2013-11-28 20:21:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cczjkd

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 绍敏郡主 at 2013-11-28 20:13:59
将虚频对应的坐标乘上一个系数比如0.2,然后与优化后的分子坐标相加,以此作为初始坐标,再进行优化

还有输入的坐标是:
W
W                    1    B1
W                    2    B2       1    A1
W                    1    B3       2    A2       3    D1       0
W                    3    B4       2    A3       1    D2       0
W                    5    B5       3    A4       2    D3       0
W                    6    B6       5    A5       3    D4       0
W                    7    B7       6    A6       5    D5       0
W                    8    B8       7    A7       6    D6       0
C                    3    B9       2    A8       1    D7       0
O                    10   B10      3    A9       2    D8       0
       Variables:
  B1                    2.55099                  
  B2                    2.68336                  
  B3                    2.5502                  
  B4                    2.64943                  
  B5                    2.6605                  
  B6                    2.64198                  
  B7                    2.51139                  
  B8                    2.49995                  
  B9                    2.59                     
  B10                   1.19                     
  A1                   62.48564                  
  A2                   60.02585                  
  A3                   61.12406                  
  A4                   60.36224                  
  A5                   88.36781                  
  A6                   62.05267                  
  A7                   59.97786                  
  A8                  109.42695                  
  A9                  178.49146                  
  D1                  109.59536                  
  D2                 -109.25265                  
  D3                  107.54165                  
  D4                   55.60575                  
  D5                   54.27173                  
  D6                 -111.17867                  
  D7                  103.93658                  
  D8                 -168.02495                  
而输出的坐标是:
Standard orientation:                        
---------------------------------------------------------------------
Center     Atomic     Atomic              Coordinates (Angstroms)
Number     Number      Type              X           Y           Z
---------------------------------------------------------------------
    1         74             0        2.114332    0.564364    1.301887
    2         74             0        2.115539    0.628731   -1.268488
    3         74             0       -0.000938    1.509963    0.035475
    4         74             0        2.244623   -1.460891   -0.035160
    5         74             0        0.000698   -0.849630   -1.326705
    6         74             0        0.000794   -0.915975    1.289214
    7         74             0       -2.114356    0.561610    1.302208
    8         74             0       -2.242400   -1.463527   -0.035154
    9         74             0       -2.115693    0.626058   -1.268780
   10          6             0       -0.010442    3.541844    0.030611
   11          8             0       -0.016202    4.737118    0.027939
---------------------------------------------------------------------
这样话怎么对应相加呢?
4楼2013-11-28 20:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marson

银虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
gkf高: 金币+2 2013-11-28 22:24:55
gmy1990: 金币+2 2013-11-28 22:28:10
将虚频对应的坐标乘上一个系数比如0.2,然后与优化后的分子坐标相加,以此作为初始坐标,再进行优化
这只是方法一,另外也可以用GV打开含有虚频的log文件然后看虚频对应的振动方向,朝着振动方向微调构型保存重新opt即可~
思想有多远,我们就能走多远~
6楼2013-11-28 20:42:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cczjkd

木虫 (正式写手)

引用回帖:
5楼: Originally posted by 绍敏郡主 at 2013-11-28 20:38:49
在输出文件里有每个振动所对应的坐标。
用gaussview打开输出文件,进入“results”,“vibrations”,就能看到各个振动的列表,选中虚频,勾选最下面一行manual displacement,把那个数改成0.2,然后右边“save s ...

没有那功能
高斯计算负频如何修改-1
9`AG7D)V](USEN3D6U%AY[U.jpg

7楼2013-11-28 22:19:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绍敏郡主

木虫 (著名写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by cczjkd at 2013-11-28 22:19:57
没有那功能

9`AG7D)V](USEN3D6U%AY[U.jpg
...

http://muchong.com/html/201104/3035135.html
8楼2013-11-28 22:29:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
如果你是多个虚频的话,如下几种方法可尝试,建议先把你的输入文件改成直角坐标(save as,下面的cartesian打钩)
1、增加收敛精度:scf=verytight;
2、降低对称性,我一般是用nosymm打破;
3、找到最大虚频,将原始分子坐标或优化后的分子坐标移动过去,我一般乘0.5,怎么移动上面都讲了,如果你体系不大,可以用excel来做,注意坐标格式修改成输入文件那样的。
然后再优化,这几个方法可以同时使用,我一般1,2,3连用。
9楼2013-11-29 07:43:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cczjkd

木虫 (正式写手)

引用回帖:
9楼: Originally posted by 枪下游魂 at 2013-11-29 07:43:13
如果你是多个虚频的话,如下几种方法可尝试,建议先把你的输入文件改成直角坐标(save as,下面的cartesian打钩)
1、增加收敛精度:scf=verytight;
2、降低对称性,我一般是用nosymm打破;
3、找到最大虚频,将 ...

直角坐标(save as,下面的cartesian打钩)这个是在GView里面吗?找不到呢。
还有请问如何用GView提交任务啊?
10楼2013-11-29 09:53:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cczjkd 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 8/400 2026-08-25 05:12 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见