24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1373  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

春华秋实2020

金虫 (正式写手)

[求助] gauss 03 lanl2dz 计算Cs时,能计算NBO电荷吗?

我尝试计算文献的结果,采用lanl2dz计算Cs2-NO3-3结构,想知道NBO电荷,但是计算时只要有pop=NBO就出错,
“Subroutine NAOANL could not find a d-type valence orbital on atom Cs    9。”
输入文件也改了很多次,有写赝势名称的,也有直接写赝势的,无论那种方法都是同样的问题,不计算NBO就可以正常计算。所以想问问,是不是gauss03使用赝势时,不能计算NBO?
1、赝势名称
# opt  freq=raman b3lyp gen Pseudo=read geom=connectivity  iop(6/76=8) pop=nbo
分子结构
******
   Cs   0
   lanl2dz
********
(空行)
Cs  0
lanl2dz

至少空两行            
2、直接写赝势
%chk=Cs-2-NO3-3-anion-singp.chk
%nproc=8
%mem=1850MB
# b3lyp/gen  Pseudo=read pop=NBO geom=connectivity

Cs-2-NO3-3-anion

-1 1
N                 -0.00027700    1.96184600   -0.09393800
O                 -1.08709000    2.57155100   -0.18683700
O                 -0.00132000    0.68957800    0.00826800
O                  1.08752100    2.57687100   -0.10008300
N                 -6.39449300   -1.39394400    0.12925500
O                 -5.35059300   -2.04603900   -0.18674100
O                 -6.28248600   -0.15158000    0.37146600
O                 -7.49440200   -1.95563200    0.19922600
Cs                -3.23421900    0.22605100   -0.07914100
N                  6.40089600   -1.38339000   -0.01109000
O                  5.44935400   -1.78427100   -0.75165500
O                  6.19264700   -0.39683000    0.76241700
O                  7.50490200   -1.94050900   -0.04230000
Cs                 3.23074400    0.23219000    0.06532900

1 2 1.5 3 1.5 4 1.5
2 9 1.0
3
4 14 1.0
5 6 1.5 7 1.5 8 1.5
6 9 1.0
7 9 1.0
8
9
10 11 1.5 12 1.5 13 1.5
11 14 1.0
12 14 1.0
13
14

N O 0         
6-311++g(2d)   
****                                                      
Cs     0                                          
S   3   1.00                                       
      0.8709000             -1.1604798            
      0.5393000              1.4198515        
      0.1724000              0.5872905        
S   4   1.00                                 
      0.8709000              0.8237727        
      0.5393000             -1.0709181            
      0.1724000             -0.9085122        
      0.0393000              1.0774049        
S   1   1.00                                 
      0.0151000              1.0000000        
P   3   1.00                                 
      1.4680000             -0.1321033        
      0.4134000              0.6572410        
      0.1536000              0.4756456        
P   2   1.00                                 
      0.1458000             -0.1680556        
      0.0279000              1.0665746        
P   1   1.00                                 
      0.0113000              1.0000000   
****

Cs  0
CS-ECP     3     46                                
f-ul potential                                 
  6                                            
0    636.7840130             -0.0805109        
1    147.8912446            -33.4311448        
2     42.5724911           -193.9270887        
2     11.2427883            -56.8301433        
2      3.5689707            -18.6504975        
2      1.3107462             -2.5411399        
s-ul potential                                 
  3                                            
0     40.8648070              3.0499934        
1     10.9329329             37.9677514        
2      2.2294903             45.2340673        
p-ul potential                                 
  4                                            
0    118.8598171              4.6406951        
1     30.6847958             62.0313533        
2     10.1762999            129.5909481        
2      1.8272850             31.5186825        
d-ul potential                                 
  4                                            
0     16.0081704              7.3054473        
1      2.6488291             33.2785760        
2     14.5555672             53.1296948        
2      0.6703604             16.2118878        
   
   
   
                                               

                                     Sample Text
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

我从网上搜了一下,说要换成NBO5.9,我使用gauss03,能换吗?
6楼2014-11-12 09:33:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

xieshen0321

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
春华秋实2020: 金币+2 2014-11-11 22:00:03
关键字:b3lyp/genecp pop=nbo
你只写了Cs的基组,其他原子的基组呢??
2楼2014-11-11 20:58:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

我给省略了,完整的是这样的
%chk=Cs-2-NO3-3-anion.chk
%nproc=8
%mem=1850MB
# opt  freq=raman b3lyp gen Pseudo=read geom=connectivity

Cs-2-NO3-3-anion

-1 1
N                  3.46570000    0.51170000   -0.69720000
O                  2.59060000    1.12460000   -1.34940000
O                  3.14400000   -0.29600000    0.22720000
O                  4.66050000    0.72720000   -0.99250000
N                 -3.41700000    0.04800000   -0.03680000
O                 -2.86640000   -0.78110000    0.75280000
O                 -2.86530000    0.87050000   -0.83580000
O                 -4.65690000    0.04170000   -0.03870000
Cs                -0.03790000    0.04310000   -0.03710000
N                  9.22190000   -2.21910000    1.86930000
O                  8.26320000   -2.83540000    2.42830000
O                  9.17800000   -1.29100000    1.00370000
O                 10.35480000   -2.59020000    2.21240000
Cs                 6.12850000   -1.22400000    0.96130000

1 2 1.0 3 1.0 4 1.0
2 9 1.0
3 9 1.0
4
5 6 1.0 7 1.0 8 1.0
6 9 1.0
7 9 1.0
8
9
10 11 1.0 12 1.0 13 1.0
11 14 1.0
12 14 1.0
13
14

N O 0
6-311++g(2d)
****
Cs 0
lanl2dz
****

Cs 0
lanl2dz
3楼2014-11-11 21:59:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xieshen0321

木虫 (正式写手)

引用回帖:
3楼: Originally posted by 春华秋实2020 at 2014-11-11 21:59:40
我给省略了,完整的是这样的
%chk=Cs-2-NO3-3-anion.chk
%nproc=8
%mem=1850MB
# opt  freq=raman b3lyp gen Pseudo=read geom=connectivity

Cs-2-NO3-3-anion

-1 1
N                  3.46570000    ...

你是只算NBO 还是 opt freq nbo一块算的啊??
4楼2014-11-11 22:22:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见