24小时热门版块排行榜    

查看: 985  |  回复: 11

tjuZidan

木虫 (小有名气)

[交流] 【求助】关于频率计算中Warning的校正已有2人参与

优化之后计算频率,结果中存在:
Warning -- explicit consideration of  33 degrees of freedom as
           vibrations may cause significant error
我在OPT的时候加了opt=Verytight,int=ultrafine,所以结构应当是局域极小值;
在频率计算的时候加了freq=(readfc,hinderedrotor)
这种情况下要么依然存在warning,要么没有计算频率和热力学,我该怎么办?
感谢大家
回复此楼
别以后想起又是遗憾
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lb1586551

铜虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
你可以通过Gaussianview看啊,内转动是容易看出来的,我感觉频率一般在几十吧,内转动的
10楼2010-08-23 20:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

tjuZidan

木虫 (小有名气)

我用的是B3lyp/genecp
对C,H,O,N用的是6-311++g**
对过渡金属用的是Lanl2dz+f

[ Last edited by tjuZidan on 2010-8-21 at 09:21 ]
别以后想起又是遗憾
2楼2010-08-20 17:17:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lb1586551

铜虫 (小有名气)

★ ★
tjuZidan(金币+1):谢谢参与
heyo_123(金币+1):鼓励交流!! 2010-08-22 16:24:33
这个问题,主要是存在内转动时,简单的看是简谐振动可能引起错误,
3楼2010-08-20 19:17:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjuZidan

木虫 (小有名气)

是有这个可能,所以我做了hinderedrotor analysis,结果中仍然有warning,而且数目没有发生变化!
别以后想起又是遗憾
4楼2010-08-21 09:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyr3365754

木虫 (著名写手)

★ ★
tjuZidan(金币+1):谢谢参与
heyo_123(金币+1):鼓励交流!! 2010-08-22 16:24:42
引用回帖:
Originally posted by tjuZidan at 2010-08-20 17:17:51:
我用的是B3lyp/genecp
对C,H,O,N用的是6-311++g**
对过渡金属用的是Lanl2dz+f

[ Last edited by tjuZidan on 2010-8-21 at 09:21 ]

请问Lanl2Dz+f 是怎么加的,谢谢!我问了好久了,请指导
5楼2010-08-21 10:03:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjuZidan

木虫 (小有名气)

引用回帖:
Originally posted by zyr3365754 at 2010-08-21 10:03:33:



请问Lanl2Dz+f 是怎么加的,谢谢!我问了好久了,请指导

我是先在这儿https://bse.pnl.gov/bse/portal找到需要的基组,然后查文献找到所用基组的f函数,然后把f函数加到d函数之后,我用的是Pd:

****
Pd     0
S   3   1.00
      2.7870000             -1.6102393        
      1.9650000              1.8489842        
      0.6243000              0.6037492        
S   4   1.00
      2.7870000              1.3540775        
      1.9650000             -1.6780848        
      0.6243000             -0.8559381        
      0.1496000              1.0200299        
S   1   1.00
      0.0436000              1.0000000        
P   3   1.00
      5.9990000             -0.1034910        
      1.4430000              0.7456952        
      0.5264000              0.3656494        
P   2   1.00
      0.7368000              0.0763285        
      0.0899000              0.9740065        
P   1   1.00
      0.0262000              1.0000000        
D   3   1.00
      6.0910000              0.0376146        
      1.7190000              0.5200479        
      0.6056000              0.5706071        
D   1   1.00
      0.1883000              1.0000000        
****

PD     0
PD-ECP     3     28
f potential
  5
0    598.3336444             -0.0563177        
1    162.4298290            -20.1288036        
2     51.5714771           -105.8197923        
2     16.4888260            -42.5733345        
2      5.8287656             -3.6165086        
s-f potential
  5
0     73.3806304              3.0003651        
1     14.7550438             32.4350093        
2     17.8350204            459.0830383        
2     12.7111477           -868.0629029        
2      9.3292063            514.4726098        
p-f potential
  5
0     55.6689877              4.9593099        
1     64.2337771             21.1711029        
2     17.6254952            605.0560092        
2     11.9058155           -726.9641846        
2      8.5100832            396.3274883        
d-f potential
  5
0     49.9994728              3.0508745        
1     39.7477547             22.2506580        
2     11.4321366            674.8357698        
2      9.1790080          -1040.8554048        
2      7.5624429            505.9375147

然后把
F   1   1.00
      1.4720000              1.0000000
加到D函数之后,就行了
别以后想起又是遗憾
6楼2010-08-21 12:07:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lb1586551

铜虫 (小有名气)


heyo_123(金币+1):鼓励交流!! 2010-08-22 16:24:10
tjuZidan(金币+3):我找到那篇文章了,非常感谢 2010-08-23 20:07:00
简单的加hinderedrotor analysis, 是不解决问题的,好像在J. chem Phys 1998上有一篇文章,有介绍怎么处理,你也可以收一下内转动的文章
7楼2010-08-21 19:48:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyr3365754

木虫 (著名写手)

引用回帖:
Originally posted by tjuZidan at 2010-08-21 12:07:56:

我是先在这儿https://bse.pnl.gov/bse/portal找到需要的基组,然后查文献找到所用基组的f函数,然后把f函数加到d函数之后,我用的是Pd:

****
Pd     0
S   3   1.00
      2.7870000         ...

谢谢你拉,嘿嘿,你感觉加了极化和没加极化的时候,对于Pd来说算出来的差距大不大???
8楼2010-08-22 09:16:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjuZidan

木虫 (小有名气)

那么请问,如果不对小于625的频率进行校正,热力学数据还可靠吗?
别以后想起又是遗憾
9楼2010-08-23 08:32:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tjuZidan 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见