24小时热门版块排行榜    

查看: 1748  |  回复: 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的回帖

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的回帖

春华秋实2020

金虫 (正式写手)

引用回帖:
4楼: Originally posted by xieshen0321 at 2014-11-11 22:22:47
你是只算NBO 还是 opt freq nbo一块算的啊??...

我开始是想一起计算,后来发现一起算NBO出错,就先优化及计算频率了,
5楼2014-11-11 22:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

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

春华秋实2020

金虫 (正式写手)

我从网上看到有计算Fe的例子,我试了一下,可以计算NBO,难道就只是对Cs不能算,不应该这样吧?
7楼2014-11-12 14:10:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

%chk=test.chk
%nproc=8
%mem=1850MB
# b3lyp/genecp 6D 10F GFinput IOP(6/7=3) pop=nbo freq=NoRaman scf

AlFe2

0 8
Al            0.000000    0.000000    1.803136
Fe           0.000000    1.201089   -0.450784
Fe           0.000000   -1.201089   -0.450784

Al 0
6-311++G(d)
****
Fe 0
lanl2dz
****

Fe 0
lanl2dz
8楼2014-11-12 14:11:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春华秋实2020

金虫 (正式写手)

我把AlFe2的输入文件中Fe替换成其他元素,都试着算了一下,有的元素可以计算,有的就出错,能算的有Ca,Ru,Co,不能算的有K,Rb,Cs。注意了一下,能算的这些元素在计算中都没有用到l607.exe,不能算的都要使用l607.exe。应该是607的问题,但是应该怎么解决呢?网上似乎有将NBO5.0替换的,不知道是不是只有这种方法。
9楼2014-11-12 17:54:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 春华秋实2020 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:02 by 5BDX0d0WFp7t
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:26 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:01 by mibUvS8DDCwf
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 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 萨仁其其格
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见