24小时热门版块排行榜    

查看: 1687  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +12 爱看书的可乐 2026-08-26 14/700 2026-08-26 15:48 by coolysnow
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 我不理解! +8 Edward_pc 2026-08-26 13/650 2026-08-26 15:37 by Tide man
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 系统查不到 +9 董八千 2026-08-26 9/450 2026-08-26 12:30 by yanmofang
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了,明天出结果吗 +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
[基金申请] 明天应该可查了!? +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
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见