24小时热门版块排行榜    

查看: 1715  |  回复: 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 wangmingqi 2026-08-28 17/850 2026-08-29 08:43 by appleren
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +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
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+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
信息提示
请填处理意见