关于gaussian NBO3.1 的bug

nkallwar

最近用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

sobereva

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

nkallwar

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

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

nkallwar

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

测试结果：

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

多谢大神