24小时热门版块排行榜    

查看: 2716  |  回复: 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一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:49 by Aj1rhIDL5ixY
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:25 by Aj1rhIDL5ixY
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:12 by Aj1rhIDL5ixY
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 16:00 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:57 by Aj1rhIDL5ixY
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:00 by Aj1rhIDL5ixY
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[有机交流] 科研求助 10+4 我吃一只鱼 2026-09-07 4/200 2026-09-10 08:32 by 0940628
信息提示
请填处理意见