计算化学公社

标题: 求助 NBO7 $DEL与g16A IOp [打印本页]

作者
Author: lqhenwunai 时间: 2020-5-28 15:00
标题: 求助 NBO7 $DEL与g16A IOp
大家好。我现在想使用NBO7的$DEL的功能，并通过G16A获得电子结构。特别的，我想通过$DEL删除所有非LEWIS轨道。手册上说仅仅通过
$nbo  $end
$del  lewis  $end
这样的关键字即可。

但是因为我的高斯版本是G16A，所以不能直接使用手册上的关键字。
”!!!!!!!!!!!!!!!!!!!!!!!!!!
Alert -- for G16 A.03 only
!!!!!!!!!!!!!!!!!!!!!!!!!!

NBO deletions (pop=nbo6del) fail with G16 A.03 unless you use a
non-standard route.  This is fixed in G16 B.01, and users are
encouraged to upgrade to B.01 if possible“

因此我必须用非标准的IOp方法来控制输入。

我看nbo7给的例子中有
#nonstd
1/38=1,163=2,172=1/1;
2/12=2,17=6,18=5,40=1/2;
3/5=4,6=6,7=1111,11=1,25=1,30=1,116=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1,40=12,113=1,114=1,124=2103/1,12(1);
99/5=1,9=1/99;
5/5=2,7=1,13=1,38=1,48=100000/2;
6/7=2,8=2,9=2,10=2,28=1,40=3,113=1,114=1,124=2103/12(-1);
99/5=1,9=1/99;

Butadiyne...rhf/6-311++g**

0  1
H    0.000000 0.000000 2.933869
C    0.000000 0.000000 1.878441
C    0.000000 0.000000 0.693074
C    0.000000 0.000000 -0.693074
C    0.000000 0.000000 -1.878441
H    0.000000 0.000000 -2.933869

$nbo $end
$del lewis $end

但是我不清楚应该怎么通过IOp来设置。其他的我还可以在官网上查到，但是6/124=2103这里，2103是什么?还有后面括号里的正1和负1是什么？

谢谢。

作者
Author: hebrewsnabla 时间: 2020-5-28 15:10
本帖最后由 hebrewsnabla 于 2020-5-28 15:21 编辑

那你直接运行这个例子能成功么

6/124定义的是高斯和External程序之间交互的方式，诸如产生什么格式的中间文件之类。既然手册给了例子，照他给的写就是。

作者
Author: lqhenwunai 时间: 2020-5-28 15:20

hebrewsnabla 发表于 2020-5-28 15:10
那你直接运行这个例子能成功么

嗯。例子是可以正常结束的。但是似乎没有相关说明，或者是我没找对地方。

#nonstd 1/38=1,163=2,172=1/1; 2/12=2,17=6,18=5,40=1/2; 3/5=4,6=6,7=111
1,11=1,25=1,30=1,116=1/1,2,3; 4//1; 5/5=2,38=5/2; 6/7=2,8=2,9=2,10=2,2
8=1,40=12,113=1,114=1,124=2103/1,12(1); 99/5=1,9=1/99; 5/5=2,7=1,13=1,
38=1,48=100000/2; 6/7=2,8=2,9=2,10=2,28=1,40=3,113=1,114=1,124=2103/12
(-1); 99/5=1,9=1/99;
----------------------------------------------------------------------
1/38=1,163=2,172=1/1;
2/12=2,17=6,18=5,40=1/2;
3/5=4,6=6,7=1111,11=1,25=1,30=1,116=1/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1,40=12,113=1,114=1,124=2103/1,12(1);
99/5=1,9=1/99;
5/5=2,7=1,13=1,38=1,48=100000/2;
6/7=2,8=2,9=2,10=2,28=1,40=3,113=1,114=1,124=2103/12(-1);
99/5=1,9=1/99;
-----------------

作者
Author: hebrewsnabla 时间: 2020-5-28 15:24
本帖最后由 hebrewsnabla 于 2020-5-28 15:32 编辑

lqhenwunai 发表于 2020-5-28 15:20
嗯。例子是可以正常结束的。但是似乎没有相关说明，或者是我没找对地方。

#nonstd 1/38=1,163=2,172= ...

那还有什么问题

你是说输出文件没有del lewis的信息？

作者
Author: lqhenwunai 时间: 2020-5-28 16:10

hebrewsnabla 发表于 2020-5-28 15:24
那还有什么问题

你是说输出文件没有del lewis的信息？

我看不懂IOP 6/124=2103是干什么的。

作者
Author: hebrewsnabla 时间: 2020-5-28 16:24
本帖最后由 hebrewsnabla 于 2020-5-28 20:23 编辑

lqhenwunai 发表于 2020-5-28 16:10
我看不懂IOP 6/124=2103是干什么的。

2L写了。要更具体的只有源码里有，应该是subroutine RunExt里面的IFUse。但是知道这个有什么用处呢？你又不用改这个，能算不就行了。

作者
Author: biogon 时间: 2020-5-28 19:49

hebrewsnabla 发表于 2020-5-28 16:24
2L写了。要更具体的只有源码里有。

IOP 6/124在官网上描述也就一句话，没任何详细的东西

作者
Author: hebrewsnabla 时间: 2020-5-29 10:58

biogon 发表于 2020-5-28 19:49
IOP 6/124在官网上描述也就一句话，没任何详细的东西

不是所有功能都会开放给用户的

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