Multiwfn查找开壳层分子pi轨道问题

treeustcer

我在使用AICD计算一价pentacene分子的pi轨道磁感应电流时，先用Multiwfn的LMO功能找出了pi轨道的序号，再添加到gjf文件的末尾。但是Gaussian在计算时总是报错，表明找不到β轨道对应的轨道序号，不知道这是哪里出了问题？
我用LMO找出的pi轨道序号如下：
54
58
60
65
67
68
69
70
71
72
73
74
645
648
651
655
657
658
659
660
661
662
663

两位数的是α轨道，可以正常读出，三位数的是β轨道，一直会在这里报错。

granvia

一种可能是，如果用g16结合NBO做自旋非限制性计算，那么输出的fchk文件里将只有alpha轨道信息，这是g16的bug

treeustcer

granvia 发表于 2020-6-20 12:21
一种可能是，如果用g16结合NBO做自旋非限制性计算，那么输出的fchk文件里将只有alpha轨道信息，这是g16的bu ...

你的意思是即使我在后面指定要输出β轨道的数据，g16也输不出来吗？因为我不需要从fchk文件里导出来，用iop(10/93=2)就直接把磁感应电流导出到txt文件里去了，而且我用g09 D.01版算得时候也是报了同样的错误，之前我还一度以为是轨道序号两边不匹配呢。

granvia

treeustcer 发表于 2020-6-20 12:35
你的意思是即使我在后面指定要输出β轨道的数据，g16也输不出来吗？因为我不需要从fchk文件里导出来，用i ...

我指的是chk或fchk文件，而且只要不用NBO就没事

sobereva

那不叫“LMO功能”。Multiwfn的轨道定域化功能里并没有给出pi轨道序号的功能。自动判断pi轨道序号的是主功能100的子功能22。

对于非限制性的情况，Multiwfn程序中beta轨道序号在alpha序号后面，仔细看此文的说明：
使用Multiwfn观看分子轨道
http://sobereva.com/269（http://bbs.keinsci.com/thread-462-1-1.html）

我并不清楚IOp(10/93=2)怎么读取beta轨道的序号，也许是需要在alpha序号后面空一行再写beta的序号，如同guess=alter的情况。测试不通过的话可以问问ACID的作者。

当前问题和NBO没有直接关系。