24小时热门版块排行榜    

查看: 2659  |  回复: 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的回帖

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的回帖
查看全部 11 个回答

绍敏郡主

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +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的回帖

绍敏郡主

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见