应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » gauss 03 lanl2dz 计算Cs时，能计算NBO电荷吗？
91/1返回列表
查看: 1375  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

春华秋实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
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-11-11 20:44:59
已阅   回复此楼   关注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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭