计算化学公社

标题: 关于gaussian NBO3.1 的bug [打印本页]

作者
Author: nkallwar 时间: 2014-12-9 20:15
标题: 关于gaussian NBO3.1 的bug

最近用gaussian 09 做一个NBO分析，出现报错：Highest occupied NBOs are not at the beginning of the NBO list;

在网上查询了一下，说是NBO3.1的一个bug，需要将轨道占据数最高的原子放在第一个

查看了一下输出作业的占据数，首先把CR型轨道占据数最高的O放在第一位，依然出现同样的错误。然后又把BD轨道占据数最高的N放在第一个原子，还是出现同样的错误

目前输出作业中轨道占据数 BD和CR占据数最高的情况如下：

1. (0.99150) BD ( 1) N 1 - C 2
            ( 60.79%) 0.7797* N 1 s( 38.60%)p 1.59( 61.35%)d 0.00(  0.05%)
                                       -0.0001  0.6211  0.0163  0.4820  0.0223
                                          0.4374  0.0181  0.4345  0.0162  0.0153
                                          0.0091  0.0094 -0.0044 -0.0059
            ( 39.21%) 0.6262* C 2 s( 30.54%)p 2.27( 69.34%)d 0.00(  0.12%)
                                          0.0000  0.5520  0.0262 -0.5728 -0.0242
                                       -0.4873 -0.0192 -0.3557 -0.0190  0.0243
                                          0.0177  0.0157  0.0022 -0.0075

91. (0.99995) CR ( 1) O 6          s(100.00%)p 0.00(  0.00%)
                                          1.0000  0.0000  0.0000  0.0002  0.0000
                                          0.0004  0.0000 -0.0002  0.0000  0.0000
                                          0.0000  0.0000  0.0000  0.0000

报错信息：
NATURAL LOCALIZED MOLECULAR ORBITAL (NLMO) ANALYSIS:
Highest occupied NBOs are not at the beginning of the NBO list;
The NLMO program is not currently set up to handle this.
Illegal range for read:  IUnit=  7 Address=    -256 Length=          -256 File length=       9393664.
Error termination in NtrErr:
NtrErr called from NtrIO.

除了去买NBO 5.0，还有其他的解决方案么？还望高人指点

输入文件：

%nprocshared=12
%mem=800MW
%chk=nbo-Cr-asym-S3-simple-O2.chk
#p opt ub3lyp/genecp nosymm pop=(savemixed,nboread) geom=connectivity
dfstructure
0 3
N                1.35611100 2.71222500 -0.18191800
C                2.22439700 3.50406200 0.45391700
N                -2.81927900 -1.15835900 1.14442600
N                -1.45921600 -2.48917900 -1.03847500
N                2.84862500 0.65413100 1.39382600
O                0.30149300 -1.08231600 1.05365700
......
H                -1.91409500 -1.28572500 -3.39659000
1 2 1.5 52 1.0
2 41 1.5 53 1.0
3 17 1.0 18 1.5
......

C H O N 0
6-31G(d)
****
Cr 0
LANL2DZ
****

Cr 0
LANl2DZ

$nbo bndidx $end

作者
Author: sobereva 时间: 2014-12-9 21:08
把savemixed去掉，不产生NLMO也就不会报错了。
如果非要产生NLMO，就用archive关键词输出.47文件，然后再用NBO5/6来做，NBO5这里有：http://bbs.keinsci.com/forum.php ... &extra=page%3D1

作者
Author: nkallwar 时间: 2014-12-9 21:15

sobereva 发表于 2014-12-9 21:08
把savemixed去掉，不产生NLMO也就不会报错了。
如果非要产生NLMO，就用archive关键词输出.47文件，然后再 ...

好的，谢谢sob大神，我先去掉savemixed试试

作者
Author: nkallwar 时间: 2014-12-11 15:08

sobereva 发表于 2014-12-9 21:08
把savemixed去掉，不产生NLMO也就不会报错了。
如果非要产生NLMO，就用archive关键词输出.47文件，然后再 ...

测试结果：

去掉了savemixed后没有出现报错。

多谢大神

欢迎光临计算化学公社 (http://ccc.keinsci.com/)